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From experimental absorption cross sections for ex- 
citation of nuclei by x-rays, radiation widths, i.e., emission 
probabilities are deduced. For the lowest activation levels of 
In'® and Pb widths of the order of a few millivolts are ob 
tained. These widths are of the same order as those of the 
low levels of the natural radioactive nuclei. The quad- 
rupole radiation widths obtained are within the upper 
limits set by a sum formula. Formulae for the transition 
probabilities are given depending both upon energy and 
spin. A numerical constant entering into the formulae, 
which is usually calculated from assumed models, is ob- 


I. INTRODUCTION 


NTIL recently, direct information about the 

absolute values of nuclear transition prob- 
abilities and the level systems of the nuclei 
involved was obtained only for a few elements. 
Three groups of experiments have been used for 
this purpose. (1) The transition probabilities 
between low levels of natural radioactive nuclei 
were estimated by comparing the intensity of 
the long range a-particle groups from the excited 
levels of these nuclei with that of the y-rays from 
the same level. I'-widths of about a few millivolts 
were obtained. (2) The transition probabilities 
between high levels of heavy nuclei were esti- 
mated from experiments on the capture of slow 
neutrons. I'-widths of about a few tenths of a 
volt were obtained. (3) The transition probability 
between a continuous and a discrete state of 


tained from experiment. A simple expression is given, 
valid in the neighborhood of the short wave limit, for the 
intensity of x-rays produced by fast electrons. At the short 
wave limit the intensity is finite and changes but little 
both with energy and frequency in the neighborhood of the 
limit. Therefrom an excitation curve with many flat 
plateaus follows, each plateau indicating another nuclear 
level. These flat plateaus reveal directly the shape of 
isochromats, thus providing a qualitative check on the 
theory of production of x-rays by fast electrons. 


light nuclei was estimated from experiments on 
the capture of protons. I'-widths from a tenth 
to about a few volts were obtained. 

Recently, however, Waldman and Collins 
excited low metastable states of In'® and Pb 
by x-rays and electrons.'! This method is capable 
of extension to higher states and probably to 
other elements. The present note shows how the 
experimental excitation curves may be under- 
stood quantitatively and how the radiation 
widths, i.e., the emission probability, may be 
calculated from the experimental absorption 
cross sections. The nuclear levels are given 
directly by the distances between discontinuities 

! These data were furnished me prior to their publication 
by Drs. Waldman and Collins. Cf. also B. Waldman, 
G. B. Collins, E. M. Stubblefield and M. Goldhaber, Phys. 
Rev. 55, 1129 (1939); G. B. Collins and B. Waldman, 


Phys. Rev. 59, 109A (1940); B. Waldman and G. B 
Collins, Phys. Rev. 57, 338 (1940). 


325 


on. 

yur 

j108 

ng 

of 

the 

ch. 

ory 

| 
| 
| 


326 EUGENE GUTH 


in the excitation curves. The radiation widths 
obtained for In" and Pb are of the order of a 
few millivolts, i.e., of the same order as those 
of the low states of the natural radioactive nuclei. 
Incidentally, the excitation curves revealing the 
shape of the isochromats provide a qualitative 
check of the theory of the production of x-rays 
by fast electrons, particularly in the neighbor- 
hood of the short wave limit. 

For the interpretation of the experimental 
data, theoretical formulae had to be introduced 
for: (a) the intensity of the x-rays produced by 
electrons near the short wave limit; (b) the 
dependency of the multipole matrix elements 
involved upon both energy and spin. 

The importance of such semi-empirical infor- 
mation about heavy nuclei is borne out by the 
well-known fact that there are no reliable models 
for these nuclei. Calculations of the transition 
probabilities for an a-pargigle model or the 
liquid droplet model may jand were, in fact, 
carried out. But the results are not trustworthy. 

The excitation of a nucleus by photons may be 
(1) a line absorption, (2) a Raman effect, and 
(3) a photo-effect. We consider here the first 
process only.? For In" it may be concluded 
directly that a line absorption must take place. 
In the experiments of Goldhaber, Hill, and 
Szilard, y-rays from } g of Ra irradiated indium 
without producing any activity.’ It is now 
known that the energy of these y-rays is suf- 
ficient to excite In™ and the absence of activity 
is attributed to the fact that none of the known 
y-ray energies are equal to the energy of the 
activation level. The probability of excitation is 
therefore small. Thus, only irradiation by a 
continuous x-ray spectrum is likely to lead to a 
line absorption. 


II. CONNECTION BETWEEN ABSORPTION CROsS 
SECTION AND EMISSION TRANSITION 
PROBABILITY 


A general formula may be derived, connecting 
the cross section o, for the excitation in line 


* Estimates show that the cross section for Raman pro- 
cesses on nuclei is rather small and probably not observable. 
An estimate of widths from photo-effect and photo-fission 
was Carried out by V. F. Weisskopf, to be published soon. 

3M. Goldhaber, R. D. Hill, and L. Szilard, Phys. Rev. 
55, 57 (1939). Also, no effect was observed for the y-rays of 
ThC, according to a personal communication of Dr. 
Goldhaber. 


absorption of a level with the emission transition 
probability. 

Let us define o(/v) as the absorption cross 
section for excitation by one photon. Then the 
dispersion formula holds: 


2jat+1 rr, 
2(2jg+1) h?(v—vo)?+T?/4 


T is the total gamma-width of a level @ with the 
energy hy, above the ground level g. I, is the 
partial width for the resonance transition. At the 
peak of the line, o(/v) reduces to the expression 


(1) 


2(2j,+1) 


ja and j, denote the spins of the states a and g, 
respectively. Formula (1) is valid not only for 
dipole, but also for electric and magnetic 2'- and 
2'_pole radiation, if and T,=T,™ are 
taken then as widths for the corresponding 
multipole radiation. 

Let n(hv)d(hv) designate the number of x-ray 
quanta produced by one incident electron of the 
energy £ with frequencies between v and v+dy. 
hv-n(hv, E)d(hv) is then the intensity. Then the 
absorption cross section o, for production of 
nuclei in the excited state with energy hv, above 
the ground state is given by 


= f n(hva, 


over the line 
2jat+1 
=n(hv,, E)—-———-: (2a) 
2(2j,+ 1) 
or conversely 
2(2j,+1) 


= . (2b) 
Qjat+1 E) 


4 


III. CONNECTION OF THE CROSS SECTION 
WITH THE EXPERIMENTAL CROsS 
SECTION Gexp 


In the actual experiments, the cross section o. 
as defined above cannot be observed. Rather, the 
internally converted electrons are observed di- 
rectly because of the decay of a metastable state 
which lies between the ground state g and the 
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activation state a, giving an experimental cross 
section @exp. The cross section o, refers to the 
production of this activation level. From o-x, the 
cross section for production of stable nuclei in 
the metastable state ¢, may be obtained by divi- 
sion by the coefficient of internal conversion a: 


Tm =Cexp' (3) 


The next step is to connect om with og. In 
order to do that, we must know the relative 
transition probabilities, activation—ground and 
activation—metastable. The ratio, f, of the 
transition probabilities 


activation—ground 


activation—metastable 


consists of the product of two factors: f= fs: fr. 
The factor fs (statistical factor) is simply equal 
to the ratio of the statistical weights of the 
ground and metastable levels. The factor fe 
(energy factor) arises from the dependency of 
the transition probabilities upon the energy. 

Knowing this ratio, f, one has only to multiply 
it by om, the cross section for production of 
nuclei in the metastable state in order to obtain 
the excitation cross section oa: 


Ca=from; f=fs- fr. (4) 
This holds if the two transitions to be compared 


have the same multipole order as is the case for 
the level scheme of Fig. 1. However, if the 
transitions have different multipole orders, as is 
the case for the level scheme of Fig. 2, f will 
contain a third factor f; where ] denotes the 
difference in spins of the levels considered. f; 
may be taken from formula (10) of Appendix 1. 
One then has: 


fer fr (4a) 


IV. DEPENDENCY OF MULTIPOLE TRANSITION 
PROBABILITIES UPON THE ENERGY 


The factor fx consists of the ratio of two 
multipole transition probabilities which depend 
upon the energy £. 

For an electric 2'-pole or a magnetic 2'~!-pole 
we have, according to definition 


(1), 
|2 } (5) 


where | Ma Zand | M |? designate the square 
of the absolute values of the coordinate electric 
2'- and magnetic 2''-pole matrix elements, re- 
spectively. There are good arguments which 
make it very plausible that the coordinate matrix 
elements ||? are to a good approximation, in- 
dependent of the energy. Rough estimates and 
more accurate simple model calculations indicate 
this, and also a more general argument put for- 
ward by Weisskopf and Ewing‘ yields the same 
conclusion. 

We have then: 


(1) 
| Cc 


where the C and C’ depend upon the actual 
values of the matrix elements involved and are 
best taken, whenever possible, from experiment. 
Model calculations are unreliable in this respect 
and generally give too high values for the con- 
stants. For a more complete formula which, 
however, is based on more special assumptions, 
cf. Appendix 1. 


‘V. F. Weisskopf and D. H. Ewing, Phys. Rev. 57, 472 
(1940); cf. also Weisskopf, reference 2. 
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V. Assumptions ABOUT THE LEVEL SCHEME 
NEEDED FOR THE EVALUATION OF THE 
Factors f 


The values of the factors fs and fz, introduced 
in formula (4), depend upon the level scheme 
assumed, namely, upon the values of the spin 
assigned to the metastable and the activation 
level and also whether another level between the 
metastable and the activation level is assumed 
to exist or not. 

For In'® the simplest plausible scheme is 
shown in Fig. 1. On the right side of the scheme 
we inserted the spins attributed to the levels. 
All three levels may have the same parity. The 
height of the activation level corresponds to the 
threshold of the x-ray excitation and the height 
of the metastable level is the energy of the elec- 
trons arising from internal conversion of the 
y-rays emitted by the decay of the metastable 
state. Both have been observed directly. 

This system does not contain the level at 
approximately 800 kev postulated by Lawson 
and Cork,® since the evidence for it is meager. 
Contrary to the usual conclusion that the 
relatively large lifetime of In'® necessitates the 
assumption of a change in parity between the 
ground and metastable states in addition to a 
change of the spin by two, no change in parity 
is postulated.® The need for introducing a change 
in parity may be overcome either by the intro- 
duction of a second level between the metastable 
level and the activation level, or by assuming 
that the matrix elements involved are much 
smaller than is generally thought.’ 

It should be pointed out that the transition 
from the ground to the activation level may 
turn out eventually to be a magnetic or even 
an electric dipole transition. In either case, it 
would be necessary to assume the existence of 


5J. L. Lawson and J. M. Cork, Phys. Rev. 57, 982 
(1940), Fig. 12, p. 993. 

6 My attention was called by Drs. Serber and Dancoff to 
the interesting fact that the assumption of a change of 
parity between m and g fits the lifetime of m; K/L con- 
version ratio and K conversion. However, in the estimate 
of the lifetime model calculations were used. Furthermore, 
the conversion coefficients for the magnetic 2*-pole radi- 
ation were taken from nonrelativistic calculations of M. H. 
Hebb and E. Nelson, Phys. Rev. 58, 486 (1940) which 
do not hold well for In (Z =49) and conversion electrons of 
0. ne energy 
ass H. Hebb and G. E. Uhlenbeck, Physica 5, 605 


an intermediate level between the activation 
level and the metastable level. A probable level 
scheme is shown in Fig. 2 (cf. ref. 5). The indices 
e and o designate even and odd. According to 
well-known selection rules, Fig. 2 means that the 
transition g—a is assumed to be a magnetic 
dipole transition. For the height of the inter- 
mediate level i the value given by Lawson and 
Cork may be used, better evidence being lacking. 

From data on excitation by x-rays alone it is 
hardly possible to arrive at an unambiguous 
decision about the multipole character of the 
activation transition. However, additional data 
about excitation of nuclear levels by electrons 
and heavy charged particles may help in this 
decision, as will be shown at another place (cf. 
also end of Appendix 1). 

In the case of lead, it is not known to which 
of its four known isotopes, namely, Pb**, Pb®°, 
and the observed activity should 
be attributed. For the even isotopes with zero 
spins a transition having a 0-0 spin change 
between two states of different parity is also 
possible,’ in addition to transitions with changes 
of the spin. For a spin change 0-0, the internally 
converted electrons will not be monochromatic. 
Furthermore, the conversion coefficient will be 
much smaller than for a transition with change 
of the spin. Neither of these conclusions has as 
yet been experimentally tested, because of the 
short life of the metastable lead level. 

Another possibility is a 0-0 transition between 
states of the same parity by direct interaction 
between the nucleus the atomic electrons, which 
penetrate, in this case, into the nucleus. Such a 
nonradiative transition takes place, according to 
Fowler,’ from the metastable level of 1.41 Mev of 
RaC’ to its ground level. One obtains, however, 
on all plausible assumptions much too short life- 
times (r<10-% sec.) to account for the observed 
periods of lead or any other known isomers.’ 

For simplicity, it is assumed here that in the 
actual transition, a change in the spin takes 
place just as for In". A possible level scheme 
is then shown in Fig. 3. The first column on the 


8 R. G. Sachs, Phys. Rev. 57, 194 (1940). 

9R.H. Fowler, Proc. Roy. Soc. A129, 1 (1930). It may 
be, however, that the observed transition is connected with 
a spin change of two or more instead of being a 0-0 
transition. 
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Fic. 4(a). Schematic nuclear level system and (b) the 
corresponding excitation curve. The theoretically discon- 
tinuous corners of the plateaus are rounded off. 


right-hand side of the scheme refers to the even 
active isotope Pb** or Pb®* or Pb*’, the second 
to Pb®”’. 


VI. NUMBER oF X-RAY QUANTA PRODUCED BY 
Fast ELECTRONS NEAR THE SHORT 
WAVE LIMIT 


In the experiments, excitation curves were 
taken as a function of the energy of the electrons 
incident on a gold or lead target. The x-rays 
produced this way were used to excite In" or 
Pb. We need, therefore, n(hv, E), i.e., the 
number of quanta produced by one electron at 
a fixed iv, namely, hv=hy,, the value of the 
threshold (energy of activation level) as a func- 
tion of energy. I(hy.-E)=hva-n(hv.-E), as a 
function of EF at fixed v, is a so-called isochromat. 
In checking the theory of the continuous x-ray 
radiation, it is, from the experimental point of 
view, best to take isochromats. It is interesting 
to point out that the shape of an excitation curve 
in passing an activation threshold reveals exactly 
the shape of the corresponding isochromat. For in 
formula (2), u(hv., EZ) is the only quantity 
varying with E. 

According to the Born approximation the num- 
ber of quanta at a fixed v produced by one elec- 
tron would increase with increasing energy. In 
fact, the Bethe-Heitler formula gives n(hve, hva) 
=0, i.e., zero number of quanta at the short 
wave limit. However, Born’s approximation does 
not hold near the short wave limit. The exact 
theory is known only in the nonrelativistic region. 
Both Sommerfeld’s original approximation and 
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the exact formulae of Sommerfeld and Maue” 
show that the number of quanta is finite at the 
short wave limit and changes but little both with 
v and E in the immediate neighborhood of the 
threshold. As a matter of fact, the isochromats 
may even decrease with increasing energy. This 
is true both for the number of quanta integrated 
over all angles or in a certain specified direction. 
Quantitatively, one may prove that the exact 
formulae of Sommerfeld and Maue may be 
obtained from the Born approximation (which 
holds rigorously at the long wave end) by multi- 
plication with a factor f(Z; v, v’) where v and v’ 
are the velocities of the electron before and after 
collision and Z is the atomic number." 


1 —e72raZe/v 


f(Z; 2, v’)= (7a) 


1— v 


It is now a plausible assumption that the same 
factor will convert the Bethe-Heitler formula into 
a formula valid with good approximation at and 
around the threshold; i.e., 


n(hv, E)=n(hv, E)p_n -f(Z; 2,0’). (7b) 


Of course, the exact n(hvy, E) may contain 
additive terms which go to zero in the nonrela- 
tivistic limit. Calculations in which approximate 
wave functions of Sommerfeld and Maue are 
used point in this direction. However, comparison 
of n(hv, E), according to formula (7b), and of 
exact numerical calculations of Jaeger’ shows 
that such additive terms are not too important. 
In fact, formula (7b) gives for hy,~1 Mev and 
E~1 Mev an n(hy, E) differing by about 15 
percent from Jaeger’s value. In the immediate 
neighborhood of the threshold, (7b) may be 
written in the form: 


n(hv, E)=A- -—-+N-dx; 
137 hv 
A23.44 for E=1.1 Mev, (7b’) 


A=4.89 for E=20.65 Mev, 


10 A. Sommerfeld and A. W. Maue, Ann. d. Physik 23, 


589 (1935). 
2raZc/v 2naZc/v’ 


11 The factor -= converts the Born 


approximation into the old Sommerfeld formula (Sommer- 
feld, Ann. d. Physik 9, 257 (1931)) which, however, is 
itself only an approximation. 

#2 J. C. Jaeger, Nature 140, 108 (1937). 
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where JN is the number of atoms per cm* and 
dx is the thickness of the (thin) target. Further 
away from the threshold, however, the full 
formula (7b) has to be used. (7b’) does not 
depend upon E at all. 

This behavior of n(hv, E) shows that an excita- 
tion curve with many flat plateaus is to be ex- 
pected, each plateau indicating another level ; cf. 
Fig. 4(a) and (b). In fact, for In"® and Pb, the 
first plateau was observed and for In" indication 
of a second level was also observed. The “‘flat- 
ness” of the first plateau is illustrated by the 
fact that, according to (7b), n(hv, E) increases 
by less than 2 percent if E is raised from 1.1 to 
1.4 Mev. The experimental curves, of course, do 
not show sharp discontinuities at each level. 
This is just as for the isochromats, because of 
the unavoidable slight inhomogeneity of the elec- 
tron beam. 


VII. CALCULATION OF THE WIDTHS 


The experimental cross sections of Waldman 
and Collins" are: 


=1X10-* for In! 
and 
= 2.5 X10-* for Pb. 


The internal conversion coefficients are™ 
a=0.5 for In™, a~1 for Pb. 
Therefore, one has according to formula (2) 


n=2X10-" cm? for In", 
om = 2.5 X10~* cm? for Pb. 


For the simple schemes of Figs. 1 and 3, the 
statistical factors are: 


fs=(2jg+1)/(2jm+1) =5 for In" 
and 
fs=} for Pb (odd isotope). 


The energy factors are, according to the energy 
values indicated in these figures: fg =6.4 for In'® 
and fg=11.4 for Pb. Inserting these factors into 


18 These cross sections are upper limits, taken from thick 
targets approaching a thin one. The cross sections are 
taken for electron energies above the threshold. They are 
sensibly independent of the electron energies. 

4 The value of a for Pb is a plausible assumption, but 
needs verification by experiment. 


formula (4), one has 


oa? =6.4X10-* for 
oa =5.7 K10-* for Pb. 


Finally, from formula (2b), using formula (7b’) 
for n(hv, E), one obtains 


=4.09 millivolts for In", 
lr, =0.28 millivolt for Pb. 


These values for the widths are smaller than an 
upper limit set by a sum formula for quadrupole 
transitions. (Cf. Appendix 2.) 

Table I gives a comparison of the widths 
lifetimes 7? =h/T,, and C=T,@/E® of 
excited states of natural radioactive nuclei and 
those of In’ and of Pb. In view of the incom- 
pleteness of the experimental data and the 
resulting ambiguity in the assumed level scheme, 
it is somewhat encouraging to find the same 
orders of magnitude. 

In order to calculate the f-factors for the level 
scheme of Fig. 2, one has to use formula (10) 
of Appendix 1. On account of the uncertainty 
about the intermediate level and about the 
validity of formula (10) we defer, however, the 
calculation of the width. 


VIII. DoprLer EFrrect AND SELF- 
ABSORPTION 


In connection with the above evaluation of 
the widths of excited states, it is important to 
point out that both for In" and Pb, the Doppler 
widths are large compared with the radiation 
widths. The Doppler width of this level of In™ 
is about 0.7 ev. This explains the smallness of 
self-absorption of the line in question in In" 
which was observed by Waldman and Collins. 


TABLE I. Widths of the excited states of RaC’, ThC’, Pb, In". 


T,@ 
E (MILLI- tr) ¢ 
NucLteus Levet (MEv) VOLTS) (SEc.) (Mev) 
RaC’ 1 0.61 2. 2.41078 
3* 1.41 0.01 
ThC’ 1 0.73 1.5 4X1078 
2 1.80 1.0 6 X1078 6.0 X107" 
Pb 2 0.65 0.28 2.34 2.40 
2 1.10 4.09 1.55 X107-" 2.59 X10 


* Only conversion electrons observed. 


Cf. H. A. Bethe, Rev. Mod. Phys. 9, 229 (1937), 
Table XXXXIV. 
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For the Doppler width, one has 


= (2kT/ Mc*)!; 
AJ: mass of nucleus. (8) 


At low temperatures, kT in this formula has to 
be replaced'’® by k(T+30), where © designates 
the Debye temperature of the nucleus. It would 
be possible to measure the width of the activation 
level more accurately by direct absorption 
measurements at low temperatures, were it not 
for the fact that the value for the Debye tem- 
perature for In is still too high for an effective 
reduction of the Doppler width. All other ele- 
ments with known metastable levels also have 
Debye temperatures too high to allow an 
appreciable reduction of the Doppler width, e.g., 
lead (Op = 88) Kr, Xe, Ag. The @p-values for In, 
Kr, Xe were calculated by means of the Linde- 
mann formula: 


0 = 123(T,,/Av!)}, (9) 


Tm, melting point; A, atomic weight; v, atomic 
volume. No data on the specific heats at low 
temperatures are known for these three elements. 

It is a pleasure for the author to thank Drs. 
B. Waldman and G. B. Collins for discussion of 
their experiments and Drs. H. A. Bethe, W. E. 
Lamb, Jr., J. A. Wheeler and especially Dr. V. F. 
Weisskopf for helpful discussions and suggestions. 


APPENDIX 1 
More complete formula for multipole matrix elements 
A semi-empirical formula for the multipole transition 
probabilities like (2): 

T(E) =c, (2) 
gives the dependency upon E but the dependency upon / 
is not contained completely. Now if one wants to relate 
the transition probabilities for transitions of different 
multipole order, one needs a formula giving the dependency 


1 Cf. W. E. Lamb, Jr., Phys. Rev. 55, 190 (1939). The 
above formula was recommended to the author by Dr. H. 
A. Bethe. According to Dr. Lamb (personal communica- 
tion) it is a fair approximation at the two limits To. 


Lamb replaces T by é(7) and (Fig. 1, p. 198) gives a plot of 
eT) against 7/0. For the above-mentioned limits, one 
obtains (O is the well-known “order of” symbol): 


TKO. 


Therefore, the replacement of T by (7+ 0) instead of 
(T+ 30) is somewhat better. 
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upon / more accurately. It may be shown that such a more 
complete formula will have for electric 2'- or magnetic 
2'-1.pole transitions the general form: 
po (4 ‘) 

(10) 
Here j; and jx are the spins of lower and upper state with 
energies E; and E,, respectively; /= | ji—je|; E= — 
measured in units mc*?; R: nuclear radius in units of 
h/mc; B is a constant, to be taken from experiment." 
This formula, however, is based on less certain assumptions 
than (2) and may be used only with caution. 

Using (10) and the value of the width in Table I for 
fixing B, one obtains for the lifetime of the metastable 
state of Pb, assuming again the level scheme of Fig. 3: 
7“4)~67 sec. in amazing agreement with the observation 
(100 sec.). This agreement has, however, to be taken only 
with great caution, because of the uncertainties in formula 
(10), in the assumed level scheme, and in the experimental 
data, inasmuch as for In" the agreement is less good 
using either the level scheme of Fig. 1 or of Fig. 2. We do 
not want to discuss this matter further here, but defer it 
to another paper, where some more applications to isomers 
will also be treated. 


APPENDIX 2 


Sum formula for quadrupole radiation 
An upper limit for the dipole matrix element 


2 
may be derived, as is well-known, from the sum formula: 


(11) 


k 
n 

Ex. =n; namely: | Dix| 
n designates the number of radiating particles (protons) 
with mass mp. 

Similarly, an upper limit for the quadrupole matrix 

element 


On= Sr 
may be derived from the sum formula 


k 
namely: 
2h 1 
| Qual (12) 


In contrast to the dipole formula, this last formu!a contains 
itself an average value. Assuming, however, | gx«|*210° ** 
and remembering that 
4r 2? 
13 
Po = (13) 


one sees that the values of Table I satisfy our inequality. 


7 For an a-particle moving in a potential well, a formula 
of this type (with A~'B(2j;+ was given 
by Hebb and Uhlenbeck, reference 7. 
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Neutron bursts were produced by modulation of the arc 
source of a small cyclotron. Neutrons subsequently slowed 
down in paraffin were detected by a BF; ionization chamber 
and linear amplifier. Modulation of an additional amplifier 
allowed the counting of all detected neutrons as well as 
those arriving in a specified time interval. Time of flight 
experiments were carried out at 3.0 meters with ‘‘on” 
times of 50 or 100 usec. (microseconds) out of a period of 
2500 usec. The modulation of the source and detector was 
accurate to about 5 usec. The absorption of thin Cd 
(0.052 g/cm?) shows a resolved resonance absorption with 
maximum at 0.14 ev. The thick Cd (0.45 g/cm*) absorp- 
tion was also examined and it was found that the absorp- 
tion edge was not sharp, the transmission being 0.5 at 0.47 
ev. Experiments with thick absorbers of Rh and In 
showed definite evidence of resonance absorption at 1.0 


ev for both elements. Experiments were also carried out 
to study the delayed emission of thermal neutrons from 
the paraffin surrounding the source. An effective mean life 
of 170 ywsec. was found. The effect of the resolution of the 
apparatus was examined and found to be small except for 
the resonances at 1.0 ev. The effect of the delayed emission 
of thermal neutrons on the time of flight measurements 
was also examined. From the thin Cd absorption curve I, 
the total width of the resonance level, was found to be 
0.12 ev and from the cross section at resonance 


where T,, is the neutron width and 7 is the angular mo- 
mentum of the absorbing nucleus, it was found to be 
5X10~ ev. In this calculation 13 percent abundance of 
the effective isotope was assumed. 


T is now well known that the absorption of 
slow neutrons by various nuclei changes 
rapidly with the kinetic energy of the neutron. 
All known nuclei show an absorption of slow 
neutrons which increases with decreasing energy 
and in addition many show the property of 
selective absorption for slow neutrons which are 
Lelieved to lie in a narrow kinetic energy band. 
These absorption bands are believed to represent 
states of a compound nucleus and their width is 
believed to be roughly 10-§ of the excitation 
energy. Through the use of slow neutrons, one is 
permitted to get a tremendously magnified view 
of a small section of the energy spectrum in the 
region of 8X10° ev, the binding energy of a 
neutron. In principle it is possible to get a 
detailed picture of these absorption bands and 
thus to obtain useful information for the study of 
nuclear structure. There are, however, no mono- 
energetic sources of slow neutrons. Neutrons of 
the proper energy for study of these resonance 
effects are produced from the fast neutrons 
generated in nuclear disintegrations, by subjec- 
tion of the neutrons to repeated energy losses. In 
this fashion a continuous neutron energy spec- 
trum is produced and in order to investigate 


nuclear properties as a function of neutron 


* This work was reported previously in brief. Baker and 
Bacher, Phys. Rev. 57, 351 fio40) ; 57, 1076 (1940). 


energy, the energy of neutrons for the particular 
process under investigation must be determined 
either by some direct method or by some indirect 
method. 

A number of indirect methods of ascertaining 
the energy of slow neutrons have been developed, 
of which the most useful is probably the boron 
method which is founded on the very reasonable 
assumption that the absorption by boron is 
inversely proportional to the neutron velocity. 
The first direct method was the mechanical 
velocity selector of Dunning, Pegram, Fink, 
Mitchell and Segré.' Two sets of cadmium- 
covered sector disks, separated by a given dis- 
tance, absorbed neutrons in a given velocity 
interval which was determined by the speed of 
rotation of the disks. This selector was used in the 
neighborhood of thermal energy (0.025 ev) and in 
order to obtain much higher energy neutrons, 
excessive rotational speeds would have been 
necessary. The resolution of the selector depends 
both upon the geometrical arrangements and 
upon the available neutron intensity which, with 
the radon sources used, was not very large. 

With the introduction of large artificial sources 
of neutrons, direct control of the production of 


1 Dunning, Pegram, Fink, Mitchell and Segré, Phys. 
Rev. 48, 704 (1935); Phys. Rev. 49, 103 (1936); Fink, 
Phys. Rev. 50, 738 (1936). 
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neutrons was made possible. With such a modu- 
lated source, methods of measuring the neutron 
time of flight for a given distance were developed 
by Alvarez? and by Fertel, Gibbs, Moon, 
Thomson, and Wynn-Williams*® (referred to as 
FGMTW). Alvarez used the 37” California 
cyclotron and produced the desired modulation 
by control of the plate voltage on the oscillator 
tubes. Neutrons were detected by a BF; ioniza- 
tion chamber located at a distance and the 
amplifier output was short-circuited except during 
a certain interval after the burst. Modulation of 
the oscillator required large amounts of power 
and the bursts which were produced were of 
rather long duration and somewhat inaccurately 
defined, but the experiments which were carried 
out clearly demonstrated the power of the 
method. Alvarez found it possible to work with 
neutrons which had energies up to about 0.025 ev 
although the resolution was rather low at this 
energy. FGM TW produced bursts of neutrons by 
modulating the ion source of a high voltage tube 
in which deuterons were accelerated and subse- 
quently used to produce neutrons. While Alvarez 
counted only those neutrons which arrived 
during a certain interval after the burst, F@MTW 
obtained a continuous record of neutron arrival 
by use of an oscillograph. Primarily because of 
the fact that the modulation of the source was 
much easier to control, they were able to extend 
the time of flight measurements to energies of 
about 0.1 ev. 

The present work is based on the use of the arc 
type ion source‘ in the cyclotron and neutron 
bursts were produced by modulation of this 
source. The method of modulation therefore was 
in principle similar to that used by FGMTW. It 
seemed difficult, however, to get the desired 
accuracy in determination of the time of flight 
from an oscillograph record and hence a modu- 
lated amplifier in principle similar to that em- 
ployed by Alvarez was used. .\ schematic dia- 
gram of the apparatus used in the present work is 
shown in Fig. 1. Neutrons produced by the 
cyclotron (1.4-Mev deuterons on Be) were first 
slowed down in paraffin and then detected by a 

?L. W. Alvarez, Phys. Rev. 54, 609 (1938). 

’Fertel, Gibbs, Moon, Thomson and Wynn-Williams, 
Proc. Roy. Soc. 175, 316 (1940); Nature 142, 829 (1938). 


a one Holloway and Baker, Rev. Sci. Inst. 10, 63 


BF; ionization chamber connected to a linear 
amplifier which operated continuously. The input 
thyratron of the discriminating amplifier flashed 
for each neutron detected and the total number of 
flashes counted through the scale-of-16 counter. 
The pulses arriving in the desired time interval 
were selected by the discriminating amplifier 
which was controlled by the amplifier modulator. 
The design and operation of the various com- 
ponents of the apparatus will be discussed below. 


COUNTER 
SCALE 


OF 16 


LINEAR DISCRIMINATING 
AMPLIFIER AMPLIFIER 
BF, CHAMBER 
CYCLOTRON 
AMPLIFIER 
ARC 
CONTROL 


FREQUENCY OSGLLOSCOPE 400 CYCLE 
METER STABILIZATION SUPPLY 


Fic. 1. Schematic diagram of apparatus. 


APPARATUS 
Source 


From previous work on the development of the 
arc type ion source it was known that this source 
would operate in the cyclotron from an alter- 
nating supply. Preliminary experiments on source 
modulation were therefore undertaken with the 
are operating from a 400-cycle generator. In 
order to find out how well this operation con- 
trolled the deuteron beam, the target was 
grounded through a resistance of 10‘ or 105 ohms 
and connected to an oscilloscope. With the 
oscilloscope sweep synchronized to the a.c. sup- 
ply, square pulses which rose and fell in 10 
microseconds or less, were observed and at this 
stage it was quite possible that the observed 
width of rise or fall was due to the amplifier which 
was used with the oscilloscope tube. 

It is necessary to be able to regulate the dura- 
tion of the deuteron burst. This was first done by 
connecting the arc in series with a thyratron used 
in a conventional grid-controlled rectifier circuit. 
The same alternating potential was applied to the 
grid and plate of this thyratron except that the 
phase of the potential applied to the former was 
shifted so that the grid was held negative over 
part of the operating cycle. The arc was thus 
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Fic. 2. Arccontrol circuit. R; = 50 adjustable ; Re = 30,000 
variable; R;=30,000; C,=1; C2=0.01. Resistance in 
ohms, capacity in uf. 


started when the thyratron fired and extinguished 
when the applied potential became too low. This 
method has the disadvantage that the critical 
extinguishing potential and hence the time of 
extinguishing, depends upon the gas pressure and 
filament condition in the arc and these are not 
stable over long periods of time. 

It appeared necessary both to start and stop 
the arc by some thyratron control. The circuit 
finally adopted is shown in Fig. 2. The upper 
thyratron starts the arc discharge at some time in 
the cycle determined by the phase of the applied 
grid potential and the lower thyratron short- 
circuits the arc and thereby stops the discharge. 
The timing of these two events is controlled by 
the phase-shifting network*7,R2C2. Both thyra- 
trons are extinguished when the supply voltage 
changes sign. The resistor R; is used to prevent 
excessive currents when both tubes are con- 
ducting. For most stable operation the grid 
potential should vary rapidly with time in the 
critical region so that changes in the supply 
voltage or the firing voltage of the thyratron will 
have little effect upon the time of firing. The 
large grid swing (400 volts) used to produce this 
rapid variation requires the use of the grid re- 
sistors R; in order to limit the grid current and 
prevent damage to the tubes. None of the circuit 
constants was critical except, of course, the phase 
shifting components which had to be adjusted for 
the desired duration of the burst. With the 400- 
cycle modulation used throughout this work, the 
pulse applied to the arc could be varied from 50 
microseconds (usec.) to somewhat more than 
500 usec. The voltage pulse controlling the arc 
rose or fell within 5 usec. 

The modulation of the deuteron beam controls 
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the production of fast neutrons but delays may 
be introduced in the slowing down process. The 
slowing down took place in a paratfin block in 
which the Be target was imbedded as shown in 
Fig. 3. The detecting chamber for the time of 
flight measurements was placed on a line ap- 
proximately perpendicular to the face of the 
paraffin block. Since the average energy of the 
neutrons is reduced to E/e" by n collisions with 
the hydrogen in the paraffin, about 16 collisions 
are necessary to reduce the energy from 5 Mev to 
1 ev. The variation of the mean free path with 
neutron energy somewhat complicates the calcu- 
lation of the time of slowing down. Assuming a 
mean free path of 5 cm for E>5000 ev and 1 cm 
for 5000 ev>E>0.1 ev, it takes 0.13 usec. for a 
neutron to go from 5 Mev to 5000 ev, 1.5 usec. to 
go to 1.4 ev and 4.0 usec. to go from 5 Mev to 
0.19 ev. In general the slowing down time will be 
of the order of the time for a few mean free paths 
at the final velocity or a few cm and this distance 
is small compared to the distance of source to 
detector. 

Somewhere in the region of 0.25 ev, the as- 
sumptions made above break down because the 
energy loss per collision is affected by the pres- 
ence of the hydrogen chemical bond. In addition, 
the neutrons at still lower energy are in thermal 
equilibrium with the paraffin and the delay in 
emission will depend upon their effective mean 
life in the paraffin. This latter point was in- 
vestigated and will be discussed later. 


CYCLOTRON 
CHAMBER 


PARAFFIN BLOCK 
20cm 


Fic. 3. Target arrangement. Observations were made on 
a line perpendicular to the face of the paraffin indicated at 
the bottom of figure. 
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BE, IONIZATION 
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B,C SHIELD CADMIUM PLATED 


ABSORBER POSITION 


Fic. 4. Ionization chamber and boron carbide collimator. 


Detector 


The slow neutrons were detected with a BF; 
ionization chamber and linear pulse amplifier. 
The ionization chamber was cylindrical with 
inside diameter of 5 cm and 15 cm long and was 
filled with BFs gas at atmospheric pressure. A 
collecting potential of 800 volts was used to 
collect negative charges. The chamber was 
shielded with cadmium and boron carbide and 
has a collimating tube as shown in Fig. 4. The 
boron carbide of the collimating tube was about 
1 g/cm? and about half as much for the chamber 
shield. Because of the geometrical arrangement 
and the presence of the first amplifier tube, there 
was no boron carbide shielding at the back of the 
ionization chamber although there was a cadmium 
shield not indicated in Fig. 4. The extension of 
the boron carbide shield then allowed cadmium 
penetrating neutrons to come from only a small 
solid angle. The general position of the ionization 
chamber with respect to the walls of the room and 
the cyclotron is shown in Fig. 5. The ionization 
chamber was placed so that in back it could see 
only an adjacent room, connection to the ampli- 
fier being made by a 1.5-m flexible cable. The 
collimation in front was such that at 3 m from 
the source, the view of the ionization chamber 
was restricted almost entirely to the face of the 
paraffin block shown in Fig. 3. 

Certain special precautions must be observed 
in the use of the ionization chamber and linear 
amplifier in time of flight measurements if the 
time of recording a neutron pulse at the output 
of the amplifier is to be an accurate measure of 
the time of arrival of the neutron in the ionization 
chamber. The collection of negative ions in a BF; 
chamber had been found to give rise to much 


sharper pulses, as might be expected if the nega- 
tive ions are electrons. This collection time for 
positive ions is about 250 usec. and for electrons 
about 6 usec. Since the collection time is de- 
pendent on the position of the disintegration in 
the chamber, it is highly desirable to make this 
time as short as possible. 

Certain modifications were made in the linear 
amplifier in order to obtain speed of response and 
avoid excessive smearing of the sharp pulses 
from the ionization chamber. To do this it is 
necessary to increase the frequency band passed 
by the amplifier. The time constants of the 
interstage coupling elements were reduced by 
decreasing the size of the coupling condensers to 
0.001 uf. This did not give sufficient high fre- 
quency response, so it was necessary to reduce 
the plate and grid resistors to 50,000 ohms. In 
order not to reduce the gain per stage, tubes with 
a high amplification factor were substituted for 
the conventional type 6C6. These tubes were 
type 1851, with an amplification factor of 6000, 
and were used in the second and third stages. The 
shielded wire in the cable between the first and 
second tubes had a capacity of 3 wuf per cm, and 
it was found necessary to substitute coaxial cable 
which had a capacity of a tenth this amount. The 
wider response of the amplifier made the signal to 
noise ratio somewhat less than that in con- 
ventional amplifiers. This is unavoidable as the 
noise is roughly proportional to the band width 
of the amplifier. 

Even with these steps taken to increase the 
amplifier speed a definite time lag was found. A 
sketch of a neutron pulse as it appears on an 
oscilloscope screen is shown in Fig. 6. While the 
pulse seems to start at a time 7), it is not re- 
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corded until 72 when the potential rises to the 
flash point of the thyratron. The time lag 
72—7\, was estimated to be about 25 usec. In 
the course of the experiments it was found pos- 
sible to make a measurement of the over-all time 
lag and this is considered to be more accurate. 
The output of the linear amplifier was fed into 
a separate unit called the discriminating amplifier 
which is shown in Fig. 7. The input thyratron 
(type 885) flashed for all neutrons detected. The 
bias of this tube was adjusted in the usual way to 
eliminate the noise background and was not 
changed in the course of an experiment. The 
thyratron output was fed into two channels. The 
first of these was connected to a scale-of-16 
counter and recorded all neutrons detected. The 
total number of neutrons detected was used as a 
measure of the integrated beam intensity and has 
several advantages over the usual method of 
beam current integration. Changes in the sensi- 
tivity of the detector or efficiency of the target do 
not, to a first approximation, affect the results if 
the number of timed neutrons is always referred 
to the total number. The 56 tube in this channel 
was used principally as a buffer tube to isolate 
the output circuit of the thyratron from possible 
interaction with the scale counter or external 
disturbances. The second channel was the timing 
channel, and led to the grid of the type 57 tube. 
The bias on this tube was adjusted so that the 
amplification was a critical function of the screen 
potential. The screen potential was supplied by a 
square wave generator (amplifier modulator) 
whose output was adjustable in width and phase 
so that only those pulses which occurred in a 
predetermined time interval were amplified. The 
following 56 was used to amplify and invert the 


WEIGHT OF CELING 32M 
HEIGHT OF TARGET 11M 


CYCLOTRON - 
DETECTOR 
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| | 
} | 


Fic. 5. Room plan. 


timed pulses and the final 885 operated the 
mechanical counter which recorded the number of 
timed neutrons. The circuit components of the 
output of the first 885 were chosen so that the 
pulses into the discriminating tube were sharp 
and of uniform height, which simplified the 
modulation problem. 

The square wave generator or amplifier modu- 
lator is shown in Fig. 8. The firing of the first 885 
starts the pulse which controls the discrimination 
and the firing of the second tube stops it by 
short-circuiting, in a similar manner to the arc 
control. The time of starting and stopping is 
determined by the phase-shifting network 7,C2R, 
and since the control of each tube is roughly 
independent of the other, any desired phase and 
width of pulse may be obtained. The power 
supply indicated as a battery must float, but a 
transformer rectifier system was finally used. 
The time taken for the pulse to rise from zero to 
its full value, or to fall, is less than 5 micro- 
seconds. The same considerations involving sta- 
bility that apply to the arc control apply here, so 
the grid swing used was large (300 to 500 volts). 


Calibration 


All time measurements were made on the 
screen of a 5-inch oscilloscope which was synchro- 


NEUTRON PULSE 
THYRATRON BIAS 


T 


Fic. 6. Sketch of oscilloscope trace showing neutron pulse 
which starts at 7; and fires the thyratron at 7». 


nized to the 400-cycle supply. For accurate 
measurements the oscilloscope sweep could be 
expanded in any region. With a 5-inch sweep the 
time scale was normally 500 usec./inch and this 
was frequently increased to 100 uwsec./inch. A 
central tap switch permitted one to examine the 
are control, amplifier modulator, oscillator or 
various outputs with the oscilloscope. Time 
settings with the arc control and amplifier 
modulator could be made to better than 0.05 
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TO SCALE OF 16 


TO MODULATOR 


57 
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+200 TO BIAS 


Fic. 7. Discriminating amplifier circuit. R; =40,000; R:=200,000; 100,000; R4= 4000; 
R; = 30,000; R,=20,000; R;=10,000 adjustable; C;=0.00025; C,=0.0001; C,;=0.005; Cy=1. 


Resistance in ohms, capacity in pf. 


inch on the oscilloscope screen which corre- 
sponded to 5 usec. 

Since a strictly linear sweep is difficult to ob- 
tain and the sweep frequently used in an ex- 
panded, nonlinear condition, all time measure- 
ments were referred to a 10-kilocycle oscillator 
synchronized to the 400-cycle supply. This was 
used to calibrate the oscilloscope trace. For 
greater accuracy, time measurements were made 
in 50-usec. units, corresponding to one loop on 
the oscillator wave ; i.e., the source time, detector 
time and interval between were set to be a 
certain number of 50-ysec. loops. The oscillator 
used was the so-called “‘transitron’’®> and the 
synchronization was introduced on the control 
grid of the 57 tube normally connected to the 
cathode. 

The frequency of the 400-cycle motor generator 
set was checked against a separate stable oscil- 
lator of known frequency. It was found that the 
generator frequency was constant to about 1 
percent. 

It is important to know the time lag of the 
entire time measuring apparatus. In addition to 
the time lag in the amplifier, there is a lag of the 
order of 10 usec. due to the acceleration of the 
deuterons. This lag depends on the cyclotron 
“D” potential, the 10-ysec. estimate being based 
on the assumption of 100 revolutions of ions at a 


5 C. Brunetti, Rev. Sci. Inst. 10, 85 (1939). 


frequency of 10 megacycles. The sum of all the 
time lags was measured directly by finding the 
time interval between production of the are 
pulse and the detection of the first neutrons at 
the detector. With no time lags, this interval 
should be approximately zero for the very fast 
neutrons which first reach the detector in the 
same position used for time of flight measure- 
ments. In order to increase the number of fast 
neutrons the paraffin was removed from around 
the Be target. The detector time was then 
adjusted to be just before the source time and the 
detector time changed until it overlapped the 
source time. Counts were taken as a function of 
the detector time and the results are shown in 
Fig. 9. With 30 usec. overlap of the source and 
detector the observed count was apparently just 
background, but for greater overlap the number 


Fic. 8. Amplifier modulator circuit. R,=10,000; 
R,=50,000; R;=500,000; R4g=20,000; C,;=0.1; C.=0.02; 
C;=0.005. Resistance in ohms, capacity in yf. 
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Fic. 9. Determination of over-all time lag. Neutrons counted 
vs. overlap of “‘on”’ times of source and detector. 


of counts increased approximately linearly. The 
linear increase was extrapolated to the back- 
ground and an effective time lag of 35 usec. 
obtained. This time lag must be subtracted from 
all time measurements. It is gratifying but 
possibly fortuitous that this figure agrees with 
the sum of the amplifier lag and estimated lag in 
the cyclotron. 


EXPERIMENTAL PROCEDURE 


The apparatus was put into operation at first 
with the cyclotron ion source unmodulated and 
the amplifier was then adjusted for the best ratio 
of counting rate of neutrons to background. The 
ion source was next operated from the arc control 
which was adjusted for the length of burst 
desired. When the ionization chamber was placed 
at the proper distance and the amplifier con- 
nected to the oscilloscope the neutron velocity 
spectrum could be seen, the group at thermal 
energies being particularly striking. The square 
topped wave which modulated the amplifier was 
then adjusted for the desired time of flight. 
Counting was started with a three-pole switch 
which simultaneously started the scale-of-16 
counter, the timed counter, and an electric clock. 
The observations were timed in order to be able 
to correct for noise and contamination back- 
ground. Counts were taken until the desired 
number of timed neutrons were recorded, the 
number being determined by the statistical 
accuracy desired and convenient to obtain, 
usually 10 percent. For convenience in recording, 


observations were made for an integral number of 
revolutions of the mechanical counter dial at the 
output of the scale-of-16 counter, or in units of 
960 recorded neutrons. The number of timed 
neutrons was referred to the total number of 
neutrons recorded in the same time interval. The 
time settings were made by switching the con- 
trolling pulse and the synchronized oscillator 
alternately onto the oscilloscope screen so that a 
direct comparison was made. At the end of an 
observation the time settings were checked before 
a new adjustment was made, although it seldom 
happened that any unexpected shift occurred. 
For absorption experiments observations were 
made with and without the absorber in place, 
with nochange in the timing so that if the settings 
were stable any error in the initial adjustment of 
the width of the setting had little effect on the 
measured absorption. 

For absolute absorption measurements it was 
necessary to compare the untimed neutrons with 
and without the absorber. This was necessary 
because the timed neutrons were referred to the 
untimed and the latter represented different 
numbers of incident neutrons for the runs with 
and without absorber. In order to avoid errors 
due to change in intensity of the source, this 
comparison was made by a succession of 30- 
second runs with and without absorber over a 
period of 10 minutes with the cyclotron operating 
under steady conditions. 

It was found that with a thick boron carbide 
absorber (2 g/cm*) in front of the ionization 
chamber, there was an appreciable background 
which decreased rapidly (Fig. 10), being ap- 
proximately 5 percent at 300 ywsec. after the end 
of the burst. This background is probably due to 
neutrons which are scattered around the room 
and which might penetrate the shielding and 
possibly be absorbed in the chamber. Such 
neutrons would, of course, take much longer to 
lose energy because of the long distance between 
collisions. The data were always corrected for this 
background. 

Since it is known that the neutrons have a 
considerable life in the paraffin as_ thermal 
neutrons it is to be expected that thermal 
neutrons will be emitted after the neutron burst. 
This delayed emission has a considerable effect in 
these experiments and was therefore measured 
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directly with the ionization chamber in contact 
with the paraffin. The results and their effect 
upon the resolution of the apparatus will be 
discussed. With the ionization chamber in the 
same position, the neutrons from the block were 
also investigated with a thick Cd absorber (0.45 
g/cm*) interposed. It was found that very few 
neutrons were detected after the production had 
stopped. If the gas pressure in the cyclotron was 
deliberately raised some neutrons were detected 
after the arc had stopped. With increasing gas 
pressure the number of these neutrons increased 
and the time after the extinguishing of the arc 
during which they were seen, increased until 
finally little etfect of the arc modulation could be 
detected. When the gas pressure was too high a 
bright glow could be seen between the ‘‘D’s” 
indicating that the unmodulated neutrons came 
from ions originating outside the arc. In these 
experiments, care was taken to operate the 
cyclotron chamber at sufficiently low pressure to 
avoid this difficulty. 


RESULTS 


The data were first corrected for the noise and 
contamination background taking into account 
the fact that for the timed neutrons the amplifier 
was sensitive for only a small fraction of the 
time. The timed neutron counts were then 
corrected for the fast neutron background as 
shown in Fig. 10. In the experiments carried out 
at 3 meters, this introduced an appreciable cor- 
rection only for those neutrons with energy above 
0.5 ev. For the data without the absorber the 
number of fast neutrons per unit total was 
subtracted from the number of counts per unit 
total to obtain the final result. The data with 
absorber, after being corrected for noise, were 
multiplied by the ratio of the untimed totals so 
that the numbers correspond to the same number 
of incident neutrons from the cyclotron. The 
number of fast neutrons per unit total already 
computed was then subtracted directly to obtain 
the corrected number with absorber. The ratio of 


the corrected numbers with and without absorber _ 


then gives the transmission directly and with the 
known number of atoms per cm? in the absorber, 
the cross section is obtained. The above pro- 
cedure for the fast neutron background assumes 
that this background is the same with and with- 


out the absorber. Since it was rather laborious to 
do, this was not checked experimentally, but it 
seems very reasonable. 

The time of flight associated with each point is 
the interval between the middle of the “on” 
time of the arc and the ‘‘on”’ time of the amplifier. 
These times were measured from the position of 
the controlling pulses on the oscilloscope screen 
and the 35-usec. time lag was subtracted. 

The results for a thin Cd absorber (0.052 ¢ /em) 
are shown in Fig. 11. The points indicated as 
circles were taken with a source-detector distance 
of 3.0 meters and an ‘‘on” time of source and 
detector of 100 usec. The points indicated by 
crosses were taken with a source-detector dis- 
tance of 1.5 meters and a source and detector 
“on” time of 50 usec. so the resolution is the 
same for both sets. The times of flight have in 
each case been divided by the distance of flight so 
the axis represents the time of flight for one 
meter. The absorption maximum has_ been 
clearly resolved, the maximum cross section ap- 
pearing for a one-meter time of flight of 190 usec. 
or an energy of 0.14 ev. There is a definite indi- 
cation of an ill-defined minimum in the neighbor- 
hood of 0.03 ev and below this energy the cross 
section seems to rise slowly. The statistical 
accuracy of each point is about 10 percent in the 
measurement of the transmission. No correction 
has been applied for the possible effect of the 
resolution of the apparatus or the effect of 
delayed emission of thermal neutrons, both of 
which are discussed later. 

The absorption of thick cadmium (0.45 g¢/cm*) 
was initially undertaken in order to get some 


‘ experience at higher energy (see Fig. 12). 


Knowledge of the absorption is also useful be- 
cause such a thickness of cadmium is often used 


FRACTION OF FAST NEUTRONS 


° 100 200 300 400 500 
TIME AFTER BURST {/SEC. 


Fic. 10. Fast neutron background. 
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Fic. 11. Absorption of thin cadmium (0.052 g/cm?). Cross section for all cadmium isotopes vs. time of flight for 1 meter 
with corresponding energy points shown below. © : 3 meters source detector distance; X : 1.5 meters. 


as a screen in slow neutron experiments. These 
data were taken at a distance of 3.0 meters 
between source and detector and the “‘on”’ times 
were 100 usec. The absorption at low energy is 
practically complete and extends to much higher 
energy than the thin cadmium absorption, being 
0.5 at 0.47 ev. The absorption edge is not 
sharply defined and the effect of the resolution of 
the apparatus on this edge is not large as will be 
shown later. The statistical accuracy above 0.4 ev 
is 10 percent but at lower energy it was decreased 
because of the long time involved in taking accu- 
rate data where the transmission is low. At 0.21 
ev the error is 15 percent and the lowest energy 
point has a statistical error of 25 percent. No 
correction for delayed emission of thermal 
neutrons need be considered for these data as the 


energy range is above that where the thermal . 


distribution plays an important role. 

The results of the measurements of the trans- 
mission of indium and rhodium are shown in 
Figs. 13 and 14. These data were taken at a 
distance of 3 meters between source and detector 
and with source and detector ‘‘on’’ times of 
50 wsec. which gave the highest resolution used in 
this work. Since these measurements were of an 
exploratory nature, thick absorbers (0.39 g/cm? 


for In and 0.60 g/cm? for Rh) were used. The 


transmission curves clearly show the resolution 
of the resonances which are at approximately 1.0 
ev for each absorber. The probable errors of 
about 10 percent for In and from 10 to 20 


02! 
ENERGY ev 


percent for Rh, were governed largely by the long 
time it took to get data at high energy. For a 
single point on the Rh curve at 1.0 ev, counts 
were taken for about 40 minutes without inter- 
ruption. It will appear that for high energy, a 
considerable correction for the resolution of the 
apparatus is necessary if a reasonable determi- 
nation of the cross section as a function of energy 
is wanted. The present data are not sufficiently 
accurate to warrant these corrections and no 
attempt was made to make the measurements 
absolute. 

The effective mean life or the delayed emission 
of thermal neutrons from the paratfin block was 
measured by placing the face of the ionization 
directly against the paraffin block. A burst of 
deuterons was produced for 50 wsec. and neutrons 
were counted in 50-ysec. intervals starting 
15 usec. after production had stopped. The re- 
sults of this experiment are shown in Fig. 15 
where the timed neutrons per unit total are 
plotted on a logarithmic scale against the time 
after the end of the burst. Except for the first 
two points, they seem to lie on a straight line 
indicating an exponential decay and an effective 
mean life of 170 usec. The statistical accuracy of 
the data is better than 5 percent but in this 
measurement any error in adjusting the width of 
the detector ‘‘on’’ time enters linearly into the 
results. This error was not greater than 5 usec. or 
10 percent but it is the largest source of error in 
the results. 
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RESOLUTION OF APPARATUS 


It is quite clear that experiments with higher 
resolving power would have given greater detail 
in the results, but this was not feasible with the 
present equipment largely because of intensity. 
In the absence of such information, it is desirable 
to examine in some detail how the results are 
affected by the resolving power of the apparatus. 
No entirely satisfactory method of correcting for 
the resolving power of the apparatus has been 
found, although several methods have been used 
to attack the problem which give interesting 
results and are of help in the evaluation of the 
data. 

If there is a uniform production of neutrons for 
a time 6 and recording of neutrons for the same 
interval at a time s after the production of 
neutrons has stopped, neutrons having times of 
flight between s and s+26 will be recorded and 
the most favored time of flight will be s+6. The 
function f(s,?) expressing the distribution in 
time of flight to which the apparatus is sensitive 
starts from zero at t=s and increases linearly 
until /=s+éand decreases linearly until f= s+ 26. 
With the proportionality constant C, f(s, ¢) is 
expressed by the following relations. 


t<s, 
f\s, t) =- (t—s), s<t<(s+6), 
(s+6) <t<(s+26), 
0, (s+26) 


The distribution function f(s, 2) and its relation 
to the ‘‘on”’ times is shown in Fig. 16. 
Fundamentally the problem of finding the 
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Fic. 12. Transmission of thick cadmium (0.45 g/cm?) vs. 
time of flight for 1 meter. Data were taken at 3 meters. 
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Fic. 13. Relative transmission of indium (0.39 g/em®) vs. 
energy showing effect of the resonance at 1.0 ev. Data were 
taken at 3 meters. 
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Fic. 14. Relative transmission of rhodium (0.60 ¢/em*) 
vs. energy showing effect of resonance at 1.0 ev. Data were 
taken at 3.2 meters. 


actual number of neutrons, as a function of time 
of flight, from the observed data D(s) amounts to 
finding a solution to the integral equation 

s+25 


D(s)= cf N(Of(s, OK (dt, (1) 


where f is the function just described, N(¢) is the 
neutron distribution from the source which is the 
function that is to be found, and K(t) represents 
the variation in the efficiency of the detector 
with neutron energy. If K(t) is known (for 
instance, it would be a constant times ¢ if the 1/z 
law holds for boron) the analysis can be simplified 
by considering the product K(t)N(t) as a single 
function M(t). f(s, t) can be considered as a 
function of s—t which is zero outside the limits 
s to s+26 so the integral equation can be written 


D(s)=C f “MW fls—Ddt. (2) 


With the introduction of the Fourier transforms 
this equation can be solved for M(t). This form 
of the solution is not, however, easily adapted to 
the practical problem since a function cannot 
easily be found to fit the data and in any case the 
labor involved is considerable and the data do not 
warrant such treatment. 
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NEUTRONS PER UNIT TOTAL 


30 700 BO 
TIME 


Fic. 15. Neutron decay in paraffin. Effective mean life 
70 usec. 


The problem can be simplified by introducing 
an approximation for the function .\/(¢) in the 
form of the first three terms of the Taylor's 
expansion of J/(t) about f. Making a change of 
variable t’ = ¢— (s+), one has for the distribution 
function, 

f=0 for j|t’| >6, 
f=t+6 for —é6<?'<0, 
f=6-l for <6. 


Without the primes, D.becomes 


6 
D(to) =C f fio] (to) 


The value of the constant C is not important as 
it is fixed for any f, and the transmission which is 
the interesting datum is the ratio of two obser- 
vations of D. C can be evaluated by considering 
K(t) and N(t) to be constant and unity, the 
observation then representing a unit number of 
neutrons. 


8+25 
i=cf f(s, or C=1/8?. 


The integration gives 
(4) 


The absence of a term in M’ in this expression is 
due directly to the symmetry of f(s, f) about 
t=¢. In obtaining a transmission curve the ratios 
of the numbers with and without the absorber are 
considered, so the interesting information is 
T=D,/D;,. The analysis is the same for the two 
cases except that for D2 the corresponding func- 
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tion ./2 is not the number of neutrons from the 
source but is the number transmitted by the 
absorber or /2=A(t).M, where A is the trans- 
mission of the absorber. The measured transmis- 
sion is then expressed by 
A(t) M(t)6?+ (64/12)[A”’ (0) 


5 


M(t)82-+ (6/12) 


In powers of 6, this may be written 


T=A+[(A"+2A'M’)/M](6/12) 
+higher powers. (6) 


If the resolution width is small the true trans- 
mission is given directly by the data as might be 
expected. To a second approximation the data 
are unaffected by the slope of the curve A, if the 
function .V is large and slowly varying with 
energy. J is the product of the number of 
neutrons from the source and the efficiency of the 
detector and is therefore the effective number of 
neutrons or the number seen by the detector. 
Under these conditions only the curvature of the 
transmission curve will alter the data. In the 
vicinity of a resonance line the curvature may 
possibly be so great that the data do not give a 
good picture of the true transmission unless the 
resolution of the experiment is sufficiently 
increased. 

In order to get a clearer picture of the etlect of 
resolution on the data, the problem was ap- 
proached from the opposite direction. Instead of 
attempting a solution of Eq. (1) for M(t) in 
terms of the data, hypothetical curves for N(t) 
were plotted and the corresponding functions 
D(s) were computed for several values of 6 which 
determines the resolving power. The Breit- 
Wigner theory predicts that the cross section in 
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Fic. 16. Distribution function f(s, t) showing relation to 
“on” times of source and detector. 
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the neighborhood of a resonance absorption is 
given in terms of the resonance energy Ee and 
the total width T' by the following expression® 


o(/2) (7) 
where ga» is the cross section at E= Ex. 


1 h? 
1 ) (8) 


+t 


Ty is the neutron width and 7 the angular 
momentum of the absorbing nucleus. The trans- 
mission of an absorber which has such a cross 
section will be 

1(£) =e, 


where is the number of atoms per cm?. 7(£) was 
then plotted on an axis representing time of 
flight for 1 meter. The value of on was chosen 
such that there was almost complete absorption 
at resonance, since thick absorbers were actually 
used for indium and rhodium. The energy 
Er=1.0 ev was also picked to correspond to this 
case. Curve 0 in Figs. 17 and 18 is for 7(/) under 
these circumstances with I'=0.2 and [ =0.4 ev, 
respectively. Curves 1 and 2 represent the trans- 
mission which would be observed for the resolu- 
tion indicated by the distribution functions 
shown at the bottom of each plot. They were 
computed by approximate numerical evaluation 
of the integral in Eq. (1) for a number of values 
of s and a smooth curve drawn through the 
points. These curves correspond to ‘‘on”’ times 


46 TIME PSEC. 

Fic. 17. Dispersion curve for T=0.2 and minimum 
transmission at 1.0 ev (curve 0). Effect of resolution by 
distribution functions indicated, is shown in curves 1 and 2. 


aan A. Bethe, Rev. Mod. Phys. 9, 69 (1937). See Eq. 


TRANSMISSION 


120 160 


Fic. 18. Dispersion curve for T=0.4 and minimum 
transmission at 1.0 ev (curve 0). Effect of resolution by 
distribution functions indicated, is shown in curves 1 
and 2, 


TIME OF FLIGHT MSE 


for 3 meters of 45 and 90 usec. for Fig. 17 and of 
60 and 120 usec. for Fig. 18. Curve 1 in each case 
corresponds approximately to the resolution 
which was actually used for indium and rhodium. 
The effect of the apparatus is to make the 
minimum of the transmission curve broader and 
less deep in both cases, although curve 1 follows 
the original curve reasonably well except for the 
depth of the minimum. 

The resolution width on a time of flight axis is 
constant as this is the directly measured quan- 
tity. The resolution at a particular energy, how- 
ever, increases with decreasing energy so that for 
thin cadmium, although the resolution corre- 
sponded to an ‘‘on”’ time of 100 usec. at 3 meters, 
the transmission minimum was much more 
clearly resolved. 

The efiect of the resolution on an absorption 
edge was also studied in essentially the same 
way, in order to examine the effect of resolution 
on the thick cadmium measurements. The origi- 
nal curve in this case was obtained by fitting two 
exponential functions (see curve 0, Fig. 19). 
These two functions, A(2—e-*'), Be-*', were 
joined at 7’=0.5 and the coefficients so adjusted 
that the curve corresponded to the observed 
transmission for thick cadmium. The derived 
curves are numbered to correspond to the resolu- 
tion functions indicated below. Curve 1 which is 
indistinguishable from the original curve, was 
obtained for an ‘‘on’’ time equal to the inverse of 
the decay constant of the exponentials. Curves 2 
and 3 were obtained for ‘‘on’’ times of twice and 
three times this amount. All the curves intersect 
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TRANSMISSION 


° 50 100 150 200 
TIME OF FLIGHT {ASE C. 


Fic. 19. Resolution of double exponential function 
simulating an ‘‘absorption edge.’ Curves 1, 2 and 3 give 
the effect of the resolution of the apparatus for the distri- 
bution functions indicated. 


for a transmission value of one-half so that this 
point is given correctly regardless of the resolu- 
tion used. This will, in general, be true only if the 
original transmission curve is symmetrical about 
this point. The general shape of the transmission 
curve is not greatly altered with any of the 
resolving powers illustrated. The thick cadmium 
measurement was made with the resolution 
illustrated by curve 2. The effect of the resolution 
on a sharp edge was also investigated. This, of 
course, shows considerable effect of the resolution 
but the result was a much sharper break than 
that observed for thick cadmium. 

It was also found that the resolution of the 
apparatus had no effect on the shape of a single 
exponential function, i.e., if M(t) is an ex- 
ponential, the derived function is also an 
exponential function with the same decay con- 
stant. Consequently the delayed emission results 
do not need to be corrected for resolution. 

When the time of flight is such that there is an 
appreciable contribution to the beam from the 
thermal distribution the effect of delayed emis- 
sion must be considered. The resolving power 
will be decreased since the neutron pulse will be 
broader, and the simple resolution function f(s, ¢) 
will no longer apply. Furthermore, the energy to 
be associated with a given time of flight setting 
will be changed if the delayed neutrons make an 
important contribution. The actual importance 
of the delayed neutrons depends on two con- 
siderations: (1) the relative number of cascading 
neutrons to thermal neutrons in the energy or 


time of flight interval to which the apparatus is 
sensitive, and (2) the amount of distortion of the 
delayed emission tail of the pulse due to the 
variation in the number— energy distribution of 
thermal neutrons, i.c., neutrons emitted long 
after the pulse must have a shorter time of flight 
or higher energy to arrive at the detector at the 
same time as neutrons emitted soon after the 
pulse. Above an energy E=k7T'=0.025 ev, this 
will decrease the effect of the tail as the number 
of neutrons decreases with increasing energy. The 
amount of this decrease depends on the particular 
energy setting. 

The variation in the number-energy relation- 
ship also causes a distortion of the pulse itself, 
giving more weight (in a region above E=kT) to 
the neutrons emitted in the early part of the pulse 
because of the greater number of neutrons of 
lower energy or longer time of flight in the 
distribution. This distortion also affects the 
energy to be associated with a given setting of the 
apparatus. The amount of the distortion cor- 
rection depends on the source-detector distance 
since this distance determines the spread in 
energy which corresponds to a given spread in 
time or resolution width. 

The amount of distortion of the pulse and its 
tail can be computed with the aid of the number- 
energy relationships as given by Bethe.’ These 
are expressed as N(E)dE, but to apply to this 
problem they must be expressed not in terms of 
the number as a function of energy in an energy 
interval, but in terms of time of flight in a time 
of flight interval. When they are so expressed, Q, 
the rate of neutron production taken as unity, 
and the values of the constants introduced, 
Bethe’s Eq. (440) for the number of cascading 
neutrons becomes 


(9) 
and Eq. (461) for the thermal neutrons becomes 


N’(t)dt =4.05 X 
Xexp (— 2.01 105d?/t-)dt, (10) 


where ¢ is the time of flight in wsec. from source to 
detector and d is the source detector distance in 
meters. It is not correct to use these expressions 


7H. A. Bethe, Rev. Mod. Phys. 9, 69 (1937). See Eqs. 


(440) and (461). 
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as they stand, in order to compare the number of 
thermal and cascading neutrons since equilibrium 
of thermal neutrons with the source was assumed 
in the derivation of N’(t) above, and this con- 
dition is not satisfied in the problem under 
consideration. The cascading neutrons probably 
come into equilibrium with the source in a very 
short time (order of 10 usec.) but because of the 
long mean life (approximately 200 usec.) it takes 
much longer for the equilibrium number of 
thermal neutrons to be established. For any 
time of the order of the mean life the number of 
thermal neutrons in the paraffin increases during 
the burst. As the duration of the burst used in 
the experiments was always short compared to 
the time necessary to establish this equilibrium, 
this effect was important. If the time of slowing 
down to thermal energies may be neglected (10) 
should be multiplied by a time factor 
where ¢’ is measured from the start of the burst 
and 7 is the mean life of thermal neutrons in the 
paraffin. This time factor also causes an additional 
distortion of the pulse giving more weight to the 
end of the pulse and therefore higher energy 
neutrons. At energies greater than k7 it is in 
the opposite direction to the distribution 
distortion. 

To illustrate the effect of these corrections on 
the data a sample calculation was made of the 
effect of distortion of the pulse shape and of 
delayed emission. Simplified conditions of source- 
detector distance of 3 meters, source ‘‘on"’ time 
of 100 wsec., and an infinitesimal detector ‘‘on”’ 
time were assumed. The relative number of 
neutrons from various parts of the pulse and its 
tail which would be recorded at the detector are 
given in the last five columns of Table I for time 
of flight settings from the center of the pulse 
corresponding to 0.2, 0.15, and 0.10 ev. In the 
second column the ratio of the number of 
thermal neutrons at the pulse center to the 
number of cascading neutrons is given. It may be 
seen on the basis of these assumptions that the 
thermal distribution is not important at an 
energy of 0.2 ev and while there will be some 
distortion of the pulse shape at this energy the 
effect of the delayed emission will be in the 
opposite direction and can essentially be neg- 
lected. At 0.15 ev the effect of pulse distortion is 
not great but the effect of the tail of the pulse 


begins to be appreciable and the resolution is not 
as good. The thermal neutrons, however, still 
make only about half the contribution to the 
number of recorded neutrons so the effect on the 
data is not large. At 0.10 ev the time distortion 
begins to be more important than the distribution 
distortion and the effect of the tail is more im- 
portant. Furthermore, the beam is now pre- 
dominantly thermal neutrons so the energy to 
associate with the time setting computed from 
the time of flight from the center of the pulse is 
too low by 10 percent or more. The decrease in 
resolving power at this and all lower energies, 
begins to have an appreciable effect on the data. 
In this region, however, the time of flight is 
longer (686 usec. for 3 meters) and as the energy 
decreases the effect of the increasing spread in the 
effective ‘‘on”’ time becomes less important be- 
cause of the increasing time of flight. For energics 
less than kT the distortions will be markedly 
different, but as the data do not indicate any 
rapid variation in the transmission curves in, this 
region the resolution and distortion corrections 
were neglected as small in magnitude and 
importance. 

The effect of the chemical bond on the slowing 
down process has been neglected in the above 
calculations. If the effect is appreciable it will 
probably result in a slower loss of neutron energy 
which makes the thermal neutrons less important 
in comparison to the cascading neutrons, both 
because of the increased number of cascading 
neutrons and because of the increased delay in 
the production of thermal neutrons. These calcu- 
lations might therefore be considered as showing 
the maximum expected effect, but can only be 
considered as indicating the order of magnitude 
because of the lack of detailed knowledge of the 
slowing down process. 


TABLE I. Relative numbers of neutrons from various parts of a 
100-usec. pulse and its tail with short detector time. N’/N 
gives the ratio of thermal to cascading neutrons. 


| TIME IN wSEC. FROM THE START OF THE BuRST 
| 


50 100 125 150 
ENERGY PULSE PULSE 
EV N'/N 25 CENTER END TAIL 
0.20 0.33 | 0.91 0 0.38 0.13 0.027 


1 
0.15 1.2, 0.71 1.0 0.81 0.52 0.18 
0.10 3.9 | 0.60 1.0 1.26 0.98 0.60 
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DISCUSSION OF RESULTS 


The plot of the thin cadmium data (Fig. 11) 
clearly shows the resolution of the resonance and 
gives a maximum value of the cross section of 
5.4 10-*! cm? at 0.14 ev. Below this energy the 
cross section decreases to about 2.310-*! cm? 
and then rises slowly with decreasing energy. The 
minimum appears to be in the neighborhood of 
0.03 ev but since the minimum is not well defined 
and the points are scattered, the energy of this 
minimum cross section cannot be stated accu- 
rately. The cadmium resonance unfortunately 
comes at an energy where the corrections in- 
volving the thermal group of neutrons are most 
uncertain. As previously indicated, the effect of 
delayed emission becomes noticeable somewhere 
below 0.15 ev and is an appreciable effect at 0.10 
ev. The shape of the low energy side of the 
resonance is therefore under suspicion. At still 
lower energy the time of flight is long enough that 
the effect of delayed emission and distortion of 
the pulse shape have a less important effect. The 
relative error in the energy may still be con- 
siderable. A corrected plot of the data was made 
which included an approximate correction for the 
delayed emission but not the distortion of the 
pulse, thus making a considerable over-correction 
for the effect of the thermals. The energy of the 
maximum was shifted only 0.02 ev to 0.16 ev. 
This seems to be sufficiently small so that for the 
present it is best to omit corrections for the 
delayed emission in the thin cadmium data. Also 
the ratio of cascading to thermal neutrons is not 
well known as a function of energy in this range 
because of the uncertain effect of the hydrogen 
chemical bond. With more complete information 
about the slowing down process and with more 
accurate data such corrections will be important. 

The cross section in the region of the maximum 
agrees surprisingly well for the data taken at 3 
meters and at 1.5 meters. These data also agree 
well with the data taken at still shorter distances 
in preliminary work. This seems to preclude the 
possibility of gross accidental or systematic 
errors in the measurement since timing errors 
would be expected to be constant and not to vary 
linearly with the distance from source to detector. 
Also uncorrected errors due to delayed emission 
or pulse distortion would affect the data taken at 
various distances differently. 


A measurement of the cadmium cross section 
by the time of flight method was made by 
FGMTW.* They show a high point in the cross 
section plot at 0.10 ev. This rough agreement is 
fortuitous, since the position of the high point is 
not sufficiently well established because of the 
small number of observations and the low resolu- 
tion in their results. At energies below the high 
point they find a cross section which decreases 
continuously with energy, in disagreement with 
the present results which show a minimum and 
then an increasing cross section. The increasing 
cross section has better theoretical justification, 
for at sufficiently low energy the E-! factor in the 
Breit-Wigner formula will predominate and 
predicts a‘1/v cross section. 

A determination of the cadmium resonance was 
also made by Hoffman and Livingston® who 
arrived at their results by measuring the cross 


section in boron of those neutrons which were not — 


absorbed in cadmium. This essentially amounts 
to an inference of the constants of the cadmium 
resonance without resolving the line. They found 
Er=0.18 ev and [=0.15 ev and in view of this 
very indirect method, their results are in suffi- 
cient agreement with the present results. 

In order to determine the width of the reso- 
nance, an attempt was made to fit the data witha 
dispersion curve.’ .\ point-by-point fit of the data 
to a dispersion curve was not feasible because of 
the complexity of the effect of resolution on the 
shape in the intermediate and thermal energy 
region, and because of the statistical inaccuracy 
of the data. A simple criterion, however, is to 
compare the ratio of the cross section at reso- 
nance with the cross section at some other 
energy, i.e., to solve Eq. (7) for Tl and evaluate it 
by introducing the measured values of the cross 
sections. 


(11) 


= 


The ratio of the cross section at resonance to the 
cross section at the low energy minimum gives 
=0.12 ev. Insufficient resolution would 
make the measured cross section at resonance 
lower than it actually is, but would have less 
effect on the minimum as it is much broader on a 


8 J. G. Hoffman and M. S. Livingston, Phys. Rev. 52, 


1228 (1937). 
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time of flight axis. The observed ratio is therefore 
equal to or smaller than the true value and the 
value of [ computed from it is an upper limit. 

A check on this value of the width can be made 
by considering the ratios of the cross section at 
resonance to the cross sections at higher energies. 
With the thick cadmium measurements, the 
ratio of the cross section at resonance to the cross 
sections at 0.26, 0.36, 0.52 and 0.80 ev gives 
values of I’ of 0.12, 0.17, 0.23 and 0.31 ev, 
respectively. Whether this variation of [ com- 
puted on the basis of a single resonance is real or 
whether it is entirely due to an effect of resolution 
has not been determined. Since the resolution 
correction to the thick cadmium data is about 
equal to the probable error in the measurements, 
these data could not give a conclusive answer. 
The resolving power of the apparatus could give 
rise to a discrepancy in this direction, so that the 
smaller value of I’ which agrees with the result 
from thin cadmium is probably the best. 

From a knowledge of the absolute cross section 
at resonance and the effective width, one can 
obtain the neutron width [Ty from Eq. (8). It 
seems likely that the absorption in cadmium is by 
one of the odd isotopes whose abundances are 
12.3 and 13 percent. If one assumes 13 percent 
abundance of the effective isotope, the cross 
section at resonance becomes opr’ =4.2x10-°° 
cm*. Since i= for the odd isotopes, 'y’ = 310-4, 
or 10-* ev, depending upon the angular mo- 
mentum of the resulting nucleus. 

The data for thick cadmium are not subject to 
the uncertain corrections which must be con- 
sidered for thin cadmium, since the data are 
taken at a higher energy where the effect of the 
thermal neutrons can be largely neglected. Cer- 
tainly, the effect of delayed emission does not 
require any appreciable correction above 0.15 ev. 
These data fit the sample resolution curve given 
in the previous section very well, hence a good 
idea of the amount of the resolution correction 
can be obtained (compare curves 0 and 2, Fig. 
19). The shape of the curve is not seriously 
altered by this effect, and probably at no point is 
the transmission value altered by more than 0.08. 
In general, the resolution correction is of the 
order of magnitude of the probable error in the 
measurement of the transmission value. 

What has been referred to as an absorption 


edge in cadmium is seen to be poorly defined 
since there is no energy where one can say that 
the transmission decreases rapidly. It is probably 
best merely to say that for the thickness used 
0.45 ¢/cm?), the energy for a transmission of one- 
half is 0.47 ev, and that below 0.2 ev the trans- 
mission is negligible (3 percent or less). 

The data for indium and rhodium are also free 
from the uncertainties introduced by the thermal 
distribution since the interesting parts of the 
curves are in the region of higher energy. On the 
other hand, the shapes of the curves are definitely 
affected by the resolution of the apparatus. 
Although the smallest ‘‘on’’ times were used for 
these measurements, the energy resolution is 
much poorer in this region than in the region of 
the cadmium resonance. When one compares the 
data with the curves shown in Fig. 17 it is 
evident that the absorption line is more poorly 
resolved than in the hypothetical case. Although 
the ‘‘on”’ times, actually used, correspond to case 
1 the resolution realized corresponds more nearly 
to case 2. The thickness of the absorbers used was 
such that the transmission at resonance was 
almost certainly negligible, while the data show a 
transmission at resonance probably over one- 
half. This indicates that the width of the line is 
somewhat less than was chosen for sample plots 
to illustrate the resolution where the smallest 
value of T considered was 0.2 ev. Since the 
corrections due to the low resolution were rather 
large and the data of an exploratory nature, it 
did not seem reasonable to try to make any 
resolution corrections. It appears that the width 
is probably less than 0.2 ev, and that the energy 
of the minimum in the transmission is about 1.0 
ev. The observed shape of the resonance, how- 
ever, is not entirely due to the resolving power of 
the apparatus since the resolution used was such 
that a very sharp resonance could at most affect 
four adjacent points. The breadth of the absorp- 
tion lines observed for rhodium and indium is, 
therefore, not negligible in comparison with the 
resolution width, even though the correction for 
resolution is large. No information about the 
cross section can be inferred from the data since 
absolute calibration of the transmission was not 
made. 

The energies of the indium and rhodium reso- 
nances have been determined by the boron 
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method. Goldsmith and Rasetti® found values of 
1.16 ev and 1.38 ev for rhodium and indium, 
respectively. Later measurements by Hornbostel, 
Goldsmith, and Manley" on the self- and mutual 
absorption of the resonance neutrons in rhodium 
and indium give values for [ of 0.13 and 0.07 ev 
for rhodium and indium, respectively, and the 
difference in the energies of the resonances as 
0.15 ev. While it is difficult to see how the time of 
flight measurements could be in error by more 
than 0.2 ev, in view of the uncertainties in both 
methods the results are in good agreement. The 
widths from absorption measurements are well 
within the upper limit set by these results, and 
the difference in energy of 0.15 ev is within the 
experimental error of this work. 

The time of flight measurements give a direct 
determination of the resonance energy and the 
approximate agreement with the results of the 
boron method constitute a check:on the 1/z 
absorption in boron. A direct measurement of the 
boron absorption was made by Fertel et al. with 
results inconsistent with the 1/v law. The most 
serious disagreement is at low energy where they 
found much too low a cross section. This disagree- 
ment is in the same energy region and in the 
same direction as their disagreement with the 
results of this work on cadmium, so that one is 
tempted to suspect some systematic error in their 
results at low energy. There were some experi- 
mental difficulties, and the authors state their 
regret that they were unable to pursue their 
investigations further. In view of the many good 
reasons for expecting the 1/v law to hold, in 
addition to the check in this work, their results 
cannot be considered as a basis for serious doubt 
of the 1/v absorption in boron, although it would 
be desirable to repeat the boron absorption 
experiment. 

The study of the delayed emission was under- 
taken primarily to ascertain the effect upon the 
time of flight measurements and the experi- 
mental conditions were the same except for the 
detector position. It also provides some direct 
evidence of the existence of thermal neutrons for 
some considerable period after the burst although 


( Ms. Goldsmith and F. Rasetti, Phys. Rev. 50, 328 
1 
%” Hornbostel, Goldsmith and Manley, Phys. Rev. 58, 18 


(1940) 


the geometry of the experiment was not well 
suited to an accurate determination of the mean 
life. It was found that the effective mean life of 
thermal neutrons coming out of the block is 
170 usec. This, however, is somewhat shorter 
than the true mean life because of diffusion of 
neutrons from the inside of the block to the 
surface where they escape. 

The effect of the diffusion can be estimated,* 
if one assumes as initial conditions the stationary 
neutron distribution due to a uniform production 
of thermal neutrons in an infinite slab of paraffin 
of thickness 2a. Under these conditions the 
neutron current at the surface of the slab is 


2q exp 
a odd 


q is the number of thermal neutrons produced per 
cm per sec., 7 is the true mean life of neutrons in 
the paraffin and y = 7A /2a where is the diffusion 
length, approximately 2 cm. For the source used 
a is 10 cm. The current 7 was plotted logarith- 
mically against ¢/7 and the best straight line 
approximation to this function was plotted in the 
time interval in which observations were taken. 
The slope of this line indicated that the observed 
decay or effective mean life, if the above con- 
ditions had been fulfilled, would have been 20 
percent less than the true mean life in the 
paraffin. This gives a mean life of 210 wsec., but 
the geometrical correction in the actual con- 
ditions used (Fig. 3) is probably larger than this, 
thus indicating a still longer true mean life. The 
effect of the lack of homogeneous sources would 
be in the opposite direction. These results there- 
fore cannot be used accurately to determine the 
true mean life. A determination of the mean life 
of neutrons in water was made by Frisch, v. 
Halban and Koch" who arrived at the result 
270 usec. from a measurement of the capture 
cross section. Reduced to paraffin this would be 
230 usec. in agreement with the corrected value 
above. 

The writers are grateful to H. A. Bethe and G. 
Placzek for numerous stimulating discussions 
during the course of these experiments. 


* The writers are indebted to G. Placzek for this calcu- 


lation. 
"Frisch, v. Halban and Koch, Kgl. Danske Videns. 


Selskab. Math.-fys. Medd. 15, 10 (1937). 
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Long-Lived Radioactive Carbon: C'' 


S. RuBEN AND M. D. KAMEN 
Radiation Laboratory and the Department of Chemistry, University of California, Berkeley, California 


(Received December ?, 1940) 


A long-lived activity has been obtained by deuteron bombardment of graphite. This activity 
has been shown by chemical methods to be isotopic with carbon. The same activity has been 
obtained by slow neutron irradiation of nitrogenous compounds (NH4NOs, NH4NQOs, etc.). 
The radiations consist of low energy negative electrons with a range of 19+2 mg/cm* aluminum 
(145+15 kev upper energy limit). No soft or hard gamma-rays could be detected. The nuclear 
reactions must be C™ (d, p) C™ and N™ (n, p) C4. A sample has been followed in a screen wall 
Geiger counter for 9 months and has shown no decay. From cross section and yield estimates 
the half-life of C" is probably 10°-10° years. The production, isolation, and detection of C' for 


tracer work is described. 


HERE are 5 known isotopes of carbon, 
stable C® (~99 percent abundance), stable 
C8 (~1 percent abundance), C" (8.8-second half- 
life),t C™ (21-minute half-life),? and C™ (half-life 
10*-10° vears).* The production and properties of 
C¥ will be described in this paper. In view of the 
possible application of C™ as a tracer,’ a dis- 
cussion of relative yields and specific activities 
obtainable is included together with a_ brief 
description of detection methods. 

It is of interest to review briefly the history of 
C¥, In 1934 Kurie, using a Wilson cloud chamber, 
studied the fast neutron disintegration of nitro- 
gen and observed in addition to the known® 
alpha-particle emission, another type of transmu- 
tation,® less frequent, in which singly charged 
heavy particles were emitted. Kurie suggested 
these were protons due to 


no (1) 


Later Bonner and Brubaker,’ irradiating nitro- 
gen with slow neutrons, observed a few heavy 
tracks (~1.06-cm range) which they attributed 


1 Barkas, Creutz, Delsasso, Fox and White, Phys. Rev. 
57, 562 (1940). ; 

?Cockcroft, Gilbert and Walton, Proc. Roy. Soc. 
(London) A148, 225 (1935). 

°S. Ruben and M. D. Kamen, Phys. Rev. 57, 549 (1940); 
58, 194 (1940). 

*C", despite its short half-life has already been used 
rather extensively as a tracer, cf. G. T. Seaborg, Chem. 
Rev. 27, 199 (1940). 

° N. Feather, Proc. Roy. Soc. (London) A136, 709 (1932). 
(1938) N. D. Kurie, Phys. Rev. 45, 904 (1934); 46, 330 
(1935) W. Bonner and W. M. Brubaker, Phys. Rev. 48, 469 


to alpha-particles. Burcham and Goldhaber* 
demonstrated these tracks were protons and 
suggested they originated from reaction (1). 
Bonner and Brubaker® agreed with Burcham and 
Goldhaber and moreover estimated Q;=0.58 
+0.03 Mev. From the known mass difference 
between neutron and proton, C™ was estimated 
to be unstable with respect to N“ by ~0.17 Mev 
so that it could be expected that C™ would 
be a beta-emitter of relatively long half-life. 
McMillan” found a long-lived soft radiation 
from metal scraped from inside the cyclotron 
vacuum chamber and suggested it might be due 
to C'™ formed by the reaction 


Qo. (2) 


Unfortunately, the sample was accidentally lost 
before any chemistry was performed. However, 
the upper energy limit of the beta-radiation was 
measured by Libby and Lee" using a screen wall 
Geiger counter and found to be ~10 kev. (It 
seems quite likely this activity was due to H* 
which was discovered by Alvarez and Cornog.") 
Shortly after, the writers had occasion to prepare 
N® (by deuteron bombardment of graphite) for 
use as a tracer in nitrogen fixation experiments. 
The same graphite target was used in all bom- 
bardments and after ~ 1000 u amp. hr. of 8-Mev 


8 W. E. Burcham and M. Goldhaber, Proc. Camb. Phil. 
Soc. 32, 632 (1936). 

* T. W. Bonner and W. M. Brubaker, Phys. Rev. 49, 223 
(1936); 49, 778 (1936). 

10E. M. McMillan, Phys. Rev. 49, 875A (1936). 

"W. F. Libby and D. D. Lee, Phys. Rev. 55, 245 (1939). 
(1939) W. Alvarez and R. Cornog, Phys. Rev. 56, 613 
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deuterons it was burned to CO: which was 
precipitated as CaCO;. A fraction of this ma- 
terial, when mounted on the inner surface of a 
cylinder and counted in a screen wall counter of 
the type developed by Libby," gave an effect of 
16+2 counts/min. It was possible to show this 
effect was not due to natural contamination of 
the chemical reagents used. It seemed desirable, 
however, to obtain stronger samples before 
drawing any positive conclusions. It was ap- 
parent that a prolonged bombardment was 
necessary and, due to the heavy schedule under 
which the cyclotron was operating, this was not 
possible until recently. In the meantime B.O3; 
was bombarded in a gas-tight target chamber 
with 32-Mev alpha-particles accelerated in the 
60” cyclotron in a search for C™ produced by 


Qs. (3) 


It had been observed that when B.O; is bom- 
barded with deuterons the major fraction of the 
newly formed C" is expelled from the target as 
CO.""5 After ~10 4 amp. hr. of alpha-particle 
bombardment, the gas from the target chamber 
was introduced into an ionization chamber con- 
nected to an FP-54 electrometer tube and 
amplifier circuit.!® No evidence for any long-lived 
activity was found. 

At this time Pollard’? and Bower and Bur- 
cham!§ presented additional evidence for the 
production and instability of C™. Pollard in- 
vestigated the protons from reactions (2) and (3). 
Using 3-Mev deuterons, the longest proton range 
was 81.1+2 cm yielding Q2=8.21 Mev. Bom- 
bardment of boron with ThC alphas produced 
protons with a 97.6-cm range ; Po alpha-particles 
gave a proton group of 42.1-cm range. These 
experiments yielded an average Q;=0.66 Mev 
and using these values he proposed the mass ot 
C"™ to be 14.00780+0.00020 and thus unstable 


183 W. F. Libby, Phys. Rev. 46, 196 (1934). 

ja a Shinohara and Ridenour, J. Chem. Phys. 3, 133 
1 

4 Ruben, Kamen and Hassid, J. Am. Chem. Soc. 61, 661 
(1939); 62, 3443 (1940). 

16 We are indebted to Dr. E. Segré for the loan of the 
electrometer and for aid in many of the measurements 
made with it. 

17E. Pollard, Phys. Rev. 56, 1168 (1939). 

18 J. C. Bower and W. E. Burcham, Proc. Roy. Soc. 


(London) A173, 391 (1939). 


with respect to N“ by ~300 kev. His attempts to 
detect an activity ascribable to C™ after 
~Ouzamp. hr. bombardment of graphite with 
3-Mev deuterons were unsuccessful. Bower and 
Burcham, using 0.8-\ev deuterons, found a 48- 
cm group of protons which they attributed to 
reaction (2) with Q2=6.1 Mev. Recently Bennett, 
Bonner, Hudspeth and Watt'® have confirmed 
this observation, using targets enriched in C®, 


ISOLATION OF C™ AFTER PROLONGED DEUTERON 
BOMBARDMENT 


By resorting to the method of internal targets*° 
it was finally found possible to produce appreci- 
able amounts of an activity which was proven 
chemically to be isotopic with carbon. A probe 
target, covered with a small quantity of Aquadag, 
was inserted into the vacuum chamber and 
subjected to a 40-uamp. beam of 7-8 Mev 
deuterons for ~120 hr. The probe was frequently 
removed and examined and the Aquadag renewed 
whenever necessary. The target prepared in this 
way proved unsatisfactory for prolonged bom- 
bardment and was finally removed and examined 
for activity. For a second bombardment, the 
graphite was electroplated with copper and hard 
soldered to the probe surface. The graphite was 
of a type prepared for use as an anode in a high 
power rectifier tube. This target proved more 
satisfactory and was able to withstand 500-600 u 
amperes of ion current without appreciable 
deterioration. This probe was placed far into the 
vacuum chamber to intercept high currents 
(~500 » amperes) of 3-4-Mev deuterons until 
13,000 u amp. hr. had been accumulated. 

Since high specific activities are necessary for 
many tracer experiments it is apparent that the 
bombardment of carbon enriched in C" is 
desirable. Professor H. C. Urey and Dr. D. W. 
Stewart very kindly sent us 0.63 gram of ~23 
percent C as amorphous carbon prepared by 
the reduction of enriched NaCN with magnesium 
at ~ 400°C. A portion of this carbon was placed 


19 Bennett, Bonner, Hudspeth and Watt, Phys. Rev. 58, 
478 (1940). Holloway and Moore (Phys. Rev. 57, 1086A 
(1940)) were unable to detect the protons from C(d, p)C™. 
They placed a lower limit of 1.2 10-*5 cm? on the value 
for the cross section of this reaction at 1 Mev. 

20R. R. Wilson and M. D. Kamen, Phys. Rev. 54, 1031 


(1938). 
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in a steel mold and subjected to high pressure*! 
(40 tons/in.*). The material failed to cohere and 
examination under a microscope revealed a gran- 
ular structure quite different from graphite found 
to press well. Since no equipment for very high 
temperature work was available, the carbon was 
sent to the Edgewater Works of the National 
Carbon and Carbide Corporation, at the sug- 
gestion of Dr. H. C. McPherson. In these labo- 
ratories, after considerable effort, Dr. \[cPherson 
succeeded in preparing several satisfactory tar- 
gets. The procedure employed by Dr. McPherson 
involved dilution of the amorphous carbon with a 
suitable amount of pitch binder. The mixture was 
pressed into a mold to make a slab 2 mm thick 
and 1 cm square. The slab was baked at 1200°C 
and then “graphitized”’ (heated in vacuum) at 
2800°C to drive off all volatile matter. The 
resultant material contained 9.3 percent C™, a 
reduction by a factor of more than two from the 
starting material. (Such loss is unavoidable with 
small quantities of carbon). The slab was sprayed 
with molten copper (Schoop process) to provide a 
tenacious metallic backing for soldering the ma- 
terial to the water-cooled probe surface. It was 
necessary in some instances to repeat the spraying 
since the slabs showed a tendency to embrittle 
and disintegrate when reheated. The slabs which 
survived this treatment possessed about } the 
heat conductivity of the best graphite obtainable 
and proved to be adequate as target material. 
One slab mounted with phos-copper solder on a 
probe withstood 10 hours’ bombardment with 
300 » amperes of 3-4-Mev deuterons and showed 
no visible deterioration. 


CHEMICAL IDENTIFICATION 


The graphite was burned in a stream of oxygen 
and the resultant gases passed over hot cupric 
oxide and finally absorbed in Ca(OH)». The 
precipitated CaCO; was filtered and carefully 
washed. The filter paper containing the precipi- 
tate, when placed over a thin-walled aluminum 
Geiger counter (0.1 mm wall thickness), failed to 
show any activity. A portion (3) of this sample, 
spread out in a thin film inside a cylinder which 
was introduced into a screen wall counter, was 


_ ™ We are indebted to the Department of Civil Engineer- 
ing for cooperation in these experiments. 


active (~400 counts/min.). In the counting 
arrangement used, the cylinder containing the 
sample could be slid to and fro over the counting 
region, thus enabling alternate readings of sample 
activity and background to be taken. The pro- 
cedure was repeated using ordinary CaCQs pre- 
pared in the same way to test for activity in the 
reagents. The CaCO; precipitate obtained in this 
way showed no activity. 

The active CaCO s was treated with acid and 
the COz collected. After drying, this was intro- 
duced into a Geiger counter and gave an effect of 
~2000 counts/min. The 5-fold increase in ac- 
tivity was due to increased solid angle and elimi- 
nation of self-absorption. The active gas was 
absorbed in Ca(OH)» and reconverted into CO» 
several times with no loss in specific activity. 
However a preliminary measurement indicated 
the range of the radiations in aluminum to be 
approximately 12 mg/cm’; which is close to the 
known value for S**. The procedure used to this 
point did not exclude the possibility that part or 
all of the activity was due to S*®O».. Carrier 
(inactive) SOs was mixed with the active gas and 
the mixture treated with KMnQ, in sulfuric acid 
solution to oxidize the SO. to SO,-~. There was 
no loss in gaseous activity. In another experiment 
the mixture of the gaseous activity and inactive 
SOz was absorbed in alkali and the SO;~-~ oxi- 
dized to SO,;-~ with I;~. The COs» was liberated 
by acidifying the solution and absorbed in 
Ca(OH)». The CaCO ; was active, there being no 
loss in specific activity. 


PRODUCTION OF C™ By NEUTRON ‘TRANSMU- 
TATION OF NITROGEN 


It was realized from the very beginning that 
deuteron bombardment of carbon was not the 
best method for preparing C™ in high specific 
activity suitable for tracer work. The transmu- 
tation of N“ by slow neutrons (reaction 1) offered 
many attractive features: (1) Possibility of ob- 
taining very high specific activity since the target 
material is not carbon but nitrogen. (2) Produc- 
tion of Cis achieved as a by-product of cyclotron 
operation and thus does not interfere with other 
research. (3) Any strong neutron source will 
suffice. It is not necessary to go to high bom- 
barding voltages (3-4 Mev as is the case for 
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deuteron bombardment). (4) The nitrogenous 
material, also containing some hydrogen, can 
profitably replace water or paraffin as a shielding 
material.” 

While the slow neutron experiments are by no 
means complete, the results obtained thus far 
conclusively indicate this to be the most practical 
method for C™ production. 

Two five-gallon carboys filled with saturated 
solutions of ammonium nitrate were placed in the 
region opposite the deflector of the Berkeley 
medical cyclotron. The solutions were irradiated 
with neutrons during a period of ~6 months. 
Approximately 40,000 .amp. hr. of deuteron 
bombardment at 16 Mev of various targets 
(principally beryllium and phosphorus) took 
place in this time. A chemical analysis was per- 
formed in which carbon monoxide, carbon dioxide, 
methane, carbon, methanol, cyanide and formal- 
dehyde were added as carriers and then separated 
and analyzed for radioactivity. After vigorous 
and prolonged shaking the gaseous compounds 
were pumped off. The gases were oxidized to 
carbon dioxide and absorbed in Ca(OH)». The 
CaCO3 was very active, containing 10° counts/ 
min. By the chemical methods described above 
the activity was shown to be isotopic with carbon. 

Charcoal and other sediment was separated 
from the NH4NO; solution by filtration. This 
proved to be completely inactive. Cyanide was 
precipitated as AgCN which also was inactive. 
Methanol and formaldehyde were distilled from 
an aliquot (~10 percent) of the NH4yNO; and the 
distillate divided into 2 portions. From one 
portion the insoluble hydrazone of formaldehyde 
prepared by the addition of excess 2,4 dinitro- 
phenyl hydrazine was found to be inactive. The 
methanol from the other portion was volatilized 
and passed over hot cupric oxide. The resulting 
CO, was not active. For future experiments a 
more complete and elaborate chemical analysis 
has been devised. 

In addition to NH,NOs3 solutions, several 
samples of solid NHsNO; (~1 Ib.) as well as 
concentrated NH,4NOz solutions have been ir- 
radiated with neutrons. C™ was found only in the 


22 C, Lapointe and F. Rasetti, Phys. Rev. 58, 554 (1940), 
have found the total capture cross section of nitrogen for 
thermal neutrons to be 1.2 10~*4 cm*. Hydrogen has a 
cross section of 0.2 10-*4 cm?. 
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COUNTS PER MINUTE 


ALUMINUM /cu) 
Fic. 1. Absorption of radiation from C“ in aluminum. 


gaseous fraction. The solid NHyNQOs was irradi- 
ated in sealed soft glass containers. No volatile 
activity could be pumped off the solid even after 
addition of carrier gases. However, when the solid 
was dissolved in water the gaseous fraction con- 
tained 1000 counts/min. Evidently the gases are 
firmly trapped and are released only when the 
crystals are completely dissolved. 

A large number of other nitrogen compounds 
are being irradiated not only to determine the 
best conditions for large scale** C™ production 
but also to investigate the chemistry of the 
recoiling carbon as a function of the chemical 
environment. Experiments are now in progress 
with urea, melamine, dicyandiamide, etc.** 


PROPERTIES OF THE LONG-LIVED RADIATION 


The absorption of the radiation in Al was 
measured by wrapping the screen wall of the 
Geiger counter with thin Al foils. With the first 
CaC*O; sample (prepared by deuterons on 
carbon) which had an activity of only 400 
counts/min. (background was 120/min.), the 
activity coming through 12 mg/cm? Al was 0+6 


23 1t is a matter of common knowledge that NH4NO; is 
thermodynamically unstable. However, it seems certain 
that its solution at room temperature offers no explosion 
hazard. 

* We are indebted to the American Cyanamid Corpora- 
tion for the gift of large quantities of these materials. 
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counts/min. However with a very thin layer 
(~0.1 mg/em*) of CaC*O; (produced from 
NH,NO3) which had an activity of 5000 counts/ 
min., 0.5 percent+0.1 percent of this effect was 
still observable at 13 mg/cm*. The first estimate 
of the range was too low and it seems a more 
accurate value is 19+2 mg ‘cm?*. The absorption 
curve is shown in Fig. 1. Since the activity coming 
through 19 mg/cm* is <0.1 percent of the initial 
effect it is safe to conclude that y-rays of ~100 
kev or more are absent. If the radiation consisted 
wholly of positrons the annihilation y-rays should 
have given an effect of 100 counts/min. coming 
through the Al.*5 A comparison of the absorption 
in Cellophane and Al indicates the absence of soft 
y-rays. 3.26 mg/cm? of Cellophane absorbed 
69+4 percent while the same amount of Al 
absorbed 71+4 percent of the radiation.*® There- 
fore it seems safe to conclude negative electrons 
are being emitted and the unstable isotope is C™. 
If the upper energy limit of the beta-ray spec- 
trum is taken to be 145,000+15,000 ev and the 
mass of N" to be 14.007526,?7 the mass of C" is 
calculated to be 14.00768. It is also possible, 
with the neutron-proton mass difference,?’ to 
determine Q; =0.64 Mev. The estimate of Bonner 
and Brubaker’? of 0.58+0.03 Mev is in agreement 
with this value. 

The activity of a sample has been followed for 
9 months and no decay observed. From this fact 
alone the half-life must be longer than 30 years. 
It is possible to make rough estimates of the half- 
life from the yields and cross sections. The 
C8(d, p)C™ experiments indicate 10*-10' years 
while the p)C™ data suggest 10*-10° vears. 
A lower limit of 1000 years is not unreasonable. 
Using the same technique employed in the H’ 
Goldhaber con- 
cludes® the half-life to be longer than 300 years. 

The upper energy limit of ~ 150 kev coupled 
with the long half-life (~10* years) is quite 
unexpected. It would have been supposed that 


>This estimate is based on measurements with a 
positron emitter (Na®2) under the same conditions. 

26 A more detailed investigation of the spectrum by Dr. J. 
Halpern and Mr. R. C. Raymond, using a Wilson cloud 
chamber, is now in progress. 

*7 Hahn, Fliigge and Mattauch, Ber. 73, 1 (1940). 
usan, D. O'Neal and M. Goldhaber, Phys. Rev. 57, 1086 

). 

*°Conference on Applied Nuclear Physics, Massa- 

chusetts Institute of Technology, 1940. 


the decay of C™ to N™ would be an allowed 
transition analogous to the He®-Li® case. Dr. P. 
Morrison has calculated that on this basis, a 
maximum half-life value of the order of a day is 
to be expected. Since the actual half-life is quite 
certainly a factor of 105-10° greater it follows 
that C'-N"™ decay involves a forbidden transi- 
tion. In consequence,” it must be supposed that 
configurations with high angular momentum 
make appreciable contributions to the ground 
state of C™. A similar* situation has been found 
in the case of Be", a nucleus which might be ex- 
pected to behave like He’, i.e., decay with a short 
half-life. Actually, despite the high upper energy 
limit of the emitted electrons (~580 kev), Be" 
seems to have a very long half-life (about 10*°-10° 
years). 

It is of interest to note that because of ‘their 
long half-lives the production of macroscopic 
amounts of C™ and Be” is at hand and in the 
near future it may be possible to determine their 
nuclear spins, etc. If the half-life of C™ is as- 
sumed to be 10‘ years, then 0.1 microgram is 
present in the sample prepared from the trans- 
mutation of nitrogen. 

In conclusion, a few remarks may be made 
regarding the measurement of C™. For small 
samples (up to ~50 mg), counting as gaseous 
CO: seems to be the simplest and most sensitive 
method. Carbon dioxide, with a trace of foreign 
gas (HO, C2H;OH, Ox, etc.) has been found® in 
this laboratory to be a quite satisfactory counter 
gas. However, since CO: requires relatively high 
counting voltages, the source of very steady high 
potential may be a limiting factor. For the 
measurement of large samples of low specific 
activity the use of thin carbonate (LisCOs, etc.) 
layers inside a screen wall counter is the most 
efficient technique. 

In many respects it is unfortunate that the 
half-life of C™ is so long, since, for the production 
of strong samples, long time exposures are re- 
quired. On the other hand, it is possible to 
concentrate C™ by various isotope separation 
methods (chemical exchange, thermal diffusion, 
ultracentrifugation, etc.) and thus material of low 


% We are indebted to Professor J. R. Oppenheimer and 
Dr. L. I. Schiff for valuable discussions on this point. 
3tE. M. McMillan and S. Ruben, unpublished. 


T. H. Norris, private communication. 
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specific activity (including the end products of 
tracer experiments), can be reconcentrated. A 
given sample of C™ (and C"’) may, therefore, be 
used repeatedly. When the use of C™ and C% 
becomes widespread enough, it should be feasible 
for a commercial company to perform this task 
at a nominal cost.” 


% The Eastman Corporation, at the suggestion, and with 
the cooperation, of Professor H. C. Urey, has begun the 
concentration of N' on a small scale. It is reasonable to 
believe that the interest and demand for carbon tracers 
will far exceed the supply of the radioactive isotopes and, 
therefore, the large scale concentration of C™ is of con- 
siderable importance. 

Note added in proof.—The estimate for the C" half-life 
was made by three different methods (1) the yield of C¥ 
coupled with an extrapolation of Amaldi’s curve (cf. 
Bethe’s review) for the C'(d, p)C™ reaction and Pollard’s 
value for the yield ratio C¥(d, p)C%/C®(d, p)C®, (2) the 
C¥8(d, p)C™ cross section was estimated by comparison 
with the B'(d,”)C" yield at bombarding energies of 
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7-8 Mev, (3) the C™ yield from the slow neutron irradia- 
tion of NH,NO;3. Method (1) gave 2-104 years while (2) 
yielded 2-108 years. Schultz and Watson, Phys. Rev. 58, 
1047 (1940), have found Pollard’s C¥(d, — 
yield ratio was high by a factor of 5-10. This seems to 
bring the first two C" half-life estimates into agreement at 
10° years. The N'4(n, p)C™ experiments suggest a still 
longer half-life for C. Hence it seems reasonable to regard 
10° years as a lower limit. 
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Ionization and Dissociation of Diatomic Molecules by Electron Impact 


Homer D. HaGstRuM* AND JOHN T. TATE 
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An extension of the mass spectrometric method of study of ionization and dissociation 
processes occurring in diatomic molecules under electron impact and the results of further 
investigation of the molecules Hz, CO, NO, Ne, and O:2 are discussed. The shapes of the peaks 
which appear in the graph of resolved ion current as a function of ion accelerating potential 
are analyzed and accounted for in terms of the wave-mechanical theory of the diatomic molecule 
and the characteristics of the mass spectrometer. It is indicated how the study of the peak shape 
coupled with appearance potential measurements has made possible conclusions concerning the 
nature of the dissociation process and the structure of the molecule. Differences between pre- 
vious mass spectrometric data and those from other sources have been clarified either by ob- 
taining better data with the mass spectrometer or by formulation of a satisfactory interpretation 
of the difference based on conclusions of an analysis of the ion peak shape. Values of dissociation 
energies, ionization potentials and interpretations of the processes that occur are given. 


INTRODUCTION 


OLLISIONS of electrons of low velocity with 
diatomic molecules may produce charged 
fragments of the molecule. The data on the 
minimum electron energy necessary to produce 
these fragments, the identification of them and 
of the processes involved may permit one to 
determine the heat of dissociation of the mole- 
cule. The relative abundance of the ion products 


* Now at the Bell Telephone Laboratories, New York. 


and the probability of a process as a function of 
the electron energy are also of interest. 

In the usual mass spectrometer the ions of a 
given m ‘e value are segregated and the depend- 
ence of the ion current on the energy of the 
bombarding electrons and on the gas pressure in 
the collision chamber is determined.'~* In the 


1For a review of all work previous to 1931 see H. D. 
Smyth, Rev. Mod. Phys. 3, 347 (1931). Reviews have also 
been given by H. Kallmann and B. Rosen, Physik. Zeits. 
32, 521 (1931); W. de Groot and F. M. Penning, Handbuch 
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dissociation of diatomic molecules some frag- 
ments are given kinetic energy. This initial 
kinetic energy cannot be determined easily by 
the mass spectrometer. By means of a modified 
apparatus Lozier® studied the minimum electron 
energy at which an ion with a definite kinetic 
energy would appear. He was able to verify 
theoretical predictions on the dissociation of 
hydrogen. The retarding potential method of 
Lozier has been improved*'* and a number of 
diatomic molecules studied. The method is 
limited in that the me value of the ion cannot 
be determined. 

It is the purpose of this paper to show how an 
analysis of the peaks in the curves of ion current 
observed with a mass spectrometer yields infor- 
mation on the ionization and dissociation of di- 
atomic molecules which previously could be 
obtained in part from the retarding potential 
methods and in part from the mass spectrometer. 
Information on the potential energy curve in- 
volved in the dissociation process, a method of 
distinguishing between dissociation and ioniza- 
tion processes which yield ions with the same 
m/e values and new interpretations of measure- 
ments on the appearance potentials and relative 
abundance of ions have been obtained. 

Data on He, Os, CO, Ne, and NO are given. 
Earlier discrepancies are cleared up and some 
known dissociation processes are studied quan- 
titatively for the first time. 


EXPERIMENTAL ARRANGEMENT 


The mass spectrometer is similar to that 
described by Tate, Smith and Vaughan* and by 
Nier."' The dimensions are shown in Fig. 1. The 
ions are formed in the electron beam of ribbon 
form and are accelerated to the slit S;. Additional 


der Physik (J. Springer, Berlin, 1933), Vol. 23/1; H. 
Sponer, Molekiilspektren (J. Springer, Berlin, 1936). 
. Vaughan, Phys. Rev. 38, 1687 (1931). 
al 7, Tate, P. T. Smith and A. L. Vaughan, Phys. Rev. 
48, 525 (1935). 
usd O. Nier and E. E. Hanson, Phys. Rev. 50, 722 
: W. W. Lozier, Phys. Rev. 36, 1285 (1930). 
. T. Tate and W. W. Lozier, Phys. Rev. 39, 254 a. 
. W. Lozier, Phys. Rev. 44, 575 (1933); 45, 840 (L) 
(1934). 
®W. W. Lozier, Phys. Rev. 46, 268 (1934). 
*E. E. Hanson, Phys. Rev. 51, 86 (1937). 
_™R. Buchdahl, Phys. Rev. 57, 1071A (1940). 
"A. O. Nier, Phys. Rev. 52, 933 (1937). 
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a \ 
A+ \ 
8+ 0 30cm S, *0.037(0 025) cm 
j 
ai 0+ 0 50cm Sp* 0022(0 006) cm 
| 
A Se 200 cm 0080 cm — 


Fic. 1. Geometrical details of the mass spectrometer. 
Slits S; and S; were widened during the course of the work 
and both values of slit width are given. The width of slit 
S; is adjustable between the limits given but was usually 
kept about half open. The lengths of the slits are S; 0.8 cm, 
S, 0.8 cm, S; 1.0 cm. A drawing-out field between the back 
plate and S, of about 30 volts per cm was used. The 
dimensions of the defining slit of the electron beam are 
0.01 X0.5 cm. The plane of the beam is half-way between 
the back plate and the first ion slit. 


acceleration is given to them between the plates 
S, and S». The ion current passing through 5S; is 
plotted as a function of the accelerating field 
between the electron beam and Sb». 

For a particular curve or peak it is essential 
that the field of the electromagnet be maintained 
constant and that there be no drift in the ac- 
celerating electric fields. The electromagnet is 
operated by storage batteries and the electric 
potentials are obtained from batteries or a 1000- 
volt stabilized rectifier." No difficulty has been 
encountered in obtaining very steady operating 
conditions. The ion current is measured by a 
direct-current amplifier with a Western Electric 
electrometer tube D96475 in a modified Barth 
circuit described by Penick." 

A porcelain leak" admits gas to the system in 
which a pressure of 10™ mm Hg is ordinarily 
maintained. A tungsten filament was used with 
Ho, Ne and COs, and an oxide-coated filament 
with NO and O¢. Hydrogen, oxygen and nitrogen 
were obtained from commercial tanks and were 
admitted to the apparatus through a liquid-air 
trap. The carbon monoxide was prepared by 
treating oxalic acid with concentrated sulphuric 
acid and washing of the resultant gases with 
concentrated potassium hydroxide and alkaline 


( aE. V. Hunt and R. W. Hickman, Rev. Sci. Inst. 10, 6 
1 

13 —D. B. Penick, Rev. Sci. Inst. 6, 115 (1935). 

4 W.R. Smythe, Rev. Sci. Inst. 7, 435 (1936). 

1 E. H. Archibald, The Preparation of Pure Inorganic 
Substances (Wiley and Sons, New York, 1932). 
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Fic. 2. Peak shape for 
an ion which is formed 
with no initial velocity. 
(type 0). This particular 
peak is for N2* from No. 
The abscissas give ion 
accelerating potential 
measured from the center 
; of the peak which corre- 


sponds to a total acceler- 
ation through approxi- 
mately 500 volts. 


pyrogallol solutions. The gas was dried by passing 
it through phosphorus pentoxide and a liquid-air 
trap. The NO was prepared'® by the action of 
sulphuric acid on a solution of potassium nitrite 
and potassium iodide. The gases formed were 
passed through a dry-ice trap to remove water 
and the higher oxides of nitrogen and the NO 
collected in a liquid-air trap and fractionally 
distilled a number of times. 

For each gas a study was made of the abun- 
dances and appearance potentials of ions from 
impurities which have the same m/e as ions from 
the molecule under investigation. It was always 
possible to regulate the amounts of impurities 
present by baking. out or further purification of 
the gas sample so that their effects were entirely 
negligible. 

The minimum energy of the bombarding elec- 
tron for the observation of an ion is called its 
appearance potential. In comparing appearance 
potentials, the amplifier is adjusted to give the 
same rate of increase of ion current above the 
break for each ion. The ionization potential of 
argon was used to obtain absolute values of 
electron energies. In studying the shape of ion 
peaks, the bombarding electrons were accelerated 
through a potential drop of 100 volts. 


ANALYSIS OF PEAK SHAPES 


Certain features of the shapes of the curves of 
ion current observed with a mass spectrometer 


ane 3. Fs Fate 

have been noted by previous investigators.* 4 
The peaks for certain ions “‘trail out” on the high 
mass side. This characteristic was attributed to 
ions which possess some initial kinetic energy. 
Hustrulid, Kusch and Tate!® reported observa- 
tion of a number of satellite peaks for ions from 
benzene vapor. The separate peaks were assumed 
to be formed by ions with relatively high kinetic 
energy in a narrow range. 

The shape of an ion peak depends mainly upon 
two factors: (1) the characteristics of the instru- 
ment, (2) the distribution in initial kinetic energy 
of the ion. The ions with different initial velocities 
are not focused by our instrument at the same 
point. This characteristic permitted a study of 
these initial velocities from which some informa- 
tion on the structure of the diatomic molecule 
may be obtained. 

Instrumental factors affecting the width of an 
ion peak are the finite width of slits and dia- 
phragms, space charge effects, shape of electron 
beam and the fact that ions are formed at differ- 
ent points in the accelerating field below S, 
(Fig. 1). If the ion products have no initial energy 
other than that of thermal agitation and if ions 
with the small range of kinetic energy for this 
case are collected equally efficiently, the curve 
of ion current has the symmetrical form shown 
in Fig. 2. A curve of this form we shall designate 
as type 0. It is obtained when ions with essen- 
tially zero kinetic energy are formed as for ex- 
ample in a process involving only ionization. 

An ion which is the result of dissociation may 
have initial energy of several electron volts. 
Since ions of the same m/e value must have the 
same velocity at the entrance slit to be focused on 
the exit slit, ions with different initial velocities 
will require different accelerating potentials to 
be brought to a focus. Thus a curve of ion cur- 
rent vs. accelerating potential will have a form 
different from that of Fig. 2 and will vield in- 
formation on the initial velocities of the ions. 
From this information data on the structure of 
the molecule may be obtained. 

The effect of the initial ionic velocity apart 
from instrumental factors on the collection of 
ions will now be considered. A complete solution 
which takes all factors into account becomes 


16 A. Hustrulid, P. Kusch and J. T. Tate, Phys. Rev. 54, 
1037 (1938). 
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exceedingly involved. But one may rather easily 
carry through a simplified calculation to deter- 
mine the percent of ions formed with different 
kinetic energies which pass through the first slit 
from a given point in the electron beam. The 
results are shown in Fig. 3. The method of 
obtaining the results is described in the figure 
caption. The final over-all curve showing the 
dependence of collection on the initial kinetic 
energy is curve 3 in Fig. 3. It is evident that ions 
of very low or zero kinetic energy are much more 
efficiently collected than those of higher energy. 


Distribution in kinetic energy of ion products 


Condon" has shown that the wave-mechanical 
formulation of the Franck-Condon principle pre- 
dicts the dependence of the probability of tran- 
sition between two molecular vibrational states in 
which a photon is emitted or absorbed upon the 
value of the so-called overlap integral : 


CS vovidr}. (1) 


In this expression fy is the wave function of the 
upper state, Y, is the wave function of the lower 
state and ¢ is the nuclear separation. The case of 
absorption for which yy is the eigenfunction of a 
state above the asymptote in the continuum 
corresponds to a transition induced by electron 
impact resulting in dissociation of the molecule. 
The success of the application of the Franck- 
Condon principle to the prediction of possible 
electron impact processes in hydrogen'*-*° sug- 
gests the use of expression (1) in a more detailed 
discussion of electron impact processes in di- 
atomic molecules. 

The probability of transition produced by 
electron impact is a function of the energy of the 
transition (£,), the energy of the bombarding 
electron (E.) and of the overlap integral. For 
values of E, well above that of E, it appears that 
the probability depends mainly on the overlap 
integral. The value of E, was 100 electron volts 
which is well above the usual value of 20 to 35 
electron volts for E;. The form of the function 


17E. U. Condon, Phys. Rev. 32, 858 (1928). 

18 E. U. Condon, Phys. Rev. 35, 658A (1930). 

1” FE. U. Condon and H. D. Smyth, Proc. Nat. Acad. Sci. 
14, 871 (1928). 
(1932). Bleakney, Phys. Rev. 35, 1180 (1930); 40, 496 
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does not alter radically with E;. This is evidenced 
by measuring as a function of E, the efficiency of 
formation of ions collected at different points on 
an ion peak which involve different values of E;. 

The value of the probability of transition is 
found by methods suggested first by studies of 
the continuous absorption spectra.*!~* The prob- 
ability of finding the molecule with a given sepa- 
ration of the atoms is proportional to y,*. The 
probability of electronic transition with a given 
nuclear separation maintained is then taken as 
proportional to y,2 for that separation. If the 
upper potential energy curve representing the 
state to which the molecule is taken is known 
and the F-C principle is assumed to be valid, 
then the curve representing the distribution in 


“o 8 4 5 6 7 8 8 
KINETIC ENERGY IN VOLTS 


Fic. 3. Efficiency of collection of ions by the mass 
spectrometer as a function of initial kinetic energy. 
Investigation of the motion of ions in the collision chamber 
has been simplified by neglecting space charge and the 
magnetic field and by assuming that the ions are formed 
with velocity vectors pointing uniformly in all directions 
with zero components along the electron beam. (a) and (b) 
show the limiting paths of ions which pass through the first 
ion slit from two _ in the electron beam. Ions initially 
directed inside the angles 6; and 6; pass through the slit. 
Curve 1 gives qualitatively the relative efficiency for case 
(a), curve 2 for case (b). The second ion slit and analyzer 
diaphragms (see Fig. 1) perform a further selection by 
limiting the spread in angle at the first slit of the trajectories 
of ions which are collected. This discriminates against ions 
formed near the edges of the electron beam which, on the 
average, possess the greatest kinetic energy of the ions 
getting through the first slit. Curve 3 is the proposed 
“over-all” efficiency curve for the spectrometer used in the 
present work, 


mt G. Winans and E. C. G. Stueckelberg, Proc. Nat. 
Acad. Sci. 14, 867 (1928). 

2 A. S. Coolidge, H. M. James and R. D. Present, J. 
Chem. Phys. 4, 193 (1936). 

*3 G. E. Gibson, O. K. Rice and N.S. Bayliss, Phys. Rev. 
44, 193 (1933). 
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Fic. 4. Construction of the distribution in initial kinetic 
energy of C* ions from CO by the reflection method. 
Curve 1 is a part of the potential curve of the ground 
electronic state of the molecule; curve 2 is the square of the 
wave function of the lowest vibrational level (harmonic 
oscillator); curve 3 is a part of the upper potential curve 
to which transitions will yield curve 4, the distribution in 
kinetic energy. The ordinates of curves 2 and 4 are arbi- 
trary. The position of curve 3 has been chosen so that a 
resultant ion peak of type I is obtained (see Fig. 5). 


initial energy of the fragments of the molecule 
can be obtained. 

The method is illustrated in Fig. 4. Curve 1 is 
the potential energy curve of the ground elec- 
tronic state of the molecule, curve 2 is the square 
of the wave function (Wf ,?), curve 3 is the upper 
potential energy curve. Curve 4 is then shown 
with ordinates representing the kinetic energy of 
an ion at dissociation. The ordinate of this scale 
will be proportional to that representing the 
total kinetic energy at dissociation of the mole- 
cule and will have its zero determined by the 
position of the asymptote (infinite nuclear 
separation) of the upper potential energy curve. 
The abscissas denote the probability of produc- 
ing an ion with a given kinetic energy, and their 
values are taken as proportional to the value of 
yi? for a particular value of r. Thus the prob- 
ability is taken as proportional to the number of 
molecules having a nuclear separation such that 
when the electronic transition takes place and 
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the molecule is raised to a new energy level, that 
is to a point on the upper curve, it will break up 
and yield an ion having a given kinetic energy. 
For example, in Fig. 4 the value of the prob- 
ability that the C* ions have kinetic energy of 2 
electron volts is proportional to the height of the 
yi? curve at the nuclear separation value, which 
on curve 3 would yield an ion of that energy. In 
Fig. 4 the position of the upper potential energy 
curve is such that the most probable nuclear 
separation in the lower state (maximum yw”) re- 
sults in zero kinetic energy of the dissociated 
parts. Curve 4 then has its maximum at zero 
kinetic energy. The values of negative kinetic 
energy indicated by the same construction for 
points of greater nuclear separation correspond 
to transitions in which excitation but no dissocia- 
tion may take place. 

The position and shape of the upper potential 
curve are in general not known. The position of 
the curve is more or less arbitrary but there are, 
of course, general restrictions on its shape. It is 
so chosen that the resultant distribution obtained 
by the method just described is in accord with 
that which may be inferred from experimental 
measurement of the ion peak shape. This method 
of ‘‘working backwards”’ has also been used to 
obtain the course of the potential curve from the 
observed intensity distribution in continuous 
spectra.” 


Construction and classification of ion peak 
shapes 


Curve 4 of Fig. 4 indicates the distribution in 
energy of the ions at the moment of formation. 
Those having small velocities are collected at the 
first slit (S,) the most efficiently (Fig. 3). The 
distribution in energy of the ion current passing 
through S; may be obtained by multiplying each 
point of curve 4 (Fig. 4) by the value of the 
corresponding efficiency (Fig. 3, curve 3). Curve 
3 of Fig. 5 is the resultant curve. Because of 
finite widths of the collecting slits and of those 
in the diaphragm, and the slight differences in 
the acceleration of ions starting at different 
points in the electron beam, the actual observed 
ion peak is broader and has no discontinuity. It 
is shown in Fig. 5 (curve 4). The abscissas for 


*4N.S. Bayliss, Proc. Roy. Soc. 158, 551 (1937). 
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this curve are determined by changes in the 
values of the electric field between S; and S» 
(of the order of 1000 volts/em) which, with a 
given magnetic field, sweeps the ion beam over 
the collecting slit S3. In curve 4 of Fig. 5 the zero 
has been placed at the maximum and it is as- 
sumed that then the abscissas as shown measure 
the initial kinetic energy of the ions. 

Ion current peaks of the form of curve 4 (Fig. 
5) are called type I. They result when the ions 
formed with zero or very near zero kinetic energy 
are most probable. Examples are H* ions from 
H2, N* ions from Ne and N* ions from NO. 

Figures 6 and 7 show the construction of a 
peak for ions which result from transition to a 
potential curve lying at larger nuclear separation 
than that of Fig. 4. The resultant peak shape 
possesses two maxima, one due to the maximum 
in the “original” distribution corresponding to 
the maximum transition probability, the other 
to the rise in the efficiency of collection by the 


10 6 4 
vouTs 

Fic. 5. Construction of the peak shape from the distri- 
bution in initial kinetic energy of C+ from CO. Curve 1 is 
the distribution from Fig. 4 (curve 4); curve 2 is the 
instrumental efficiency curve; curve 3 is the predicted peak 
shape obtained by multiplying curves 1 and 2. Curve 4 is 
the observed peak shape. The abscissas of curve 4 are 

determined as they were in Fig. 2. 
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Fic. 6. Construction of the distribution in initial kinetic 
energy of O* ions from NO by the reflection method. See 
caption of Fig. 4 for the designation of the curves. The 
position of curve 3 has been chosen so that a peak of type II 
is obtained (see Fig. 7). 


spectrometer near zero kinetic energy. Such a 
peak shape, to be designated as type II, has been 
observed for O* ions from NO and is to be ex- 
pected whenever ions of kinetic energy near zero 
are present but with reduced abundance when 
comparison is made with type I. 

If the upper potential curve lies at still larger 
values of r the probability of formation of ions 
of very low or zero initial kinetic energy may be 
zero. The distribution in kinetic energy possesses 
a minimum kinetic energy greater than zero 
because the upper potential curve crosses its 
asymptote outside the Franck-Condon region. 
By virtue of the shape of the function y,?=0 for 
the ground state (curve 2 of Fig. 4) the prob- 
ability of formation of ions possessing kinetic 
energies near the minimum is relatively low. The 
instrumental efficiency is relatively constant over 
the range of values of initial kinetic energy 
possessed by the ions; hence the shape of the 
peak, designated as type III, approximates that 
of the distribution at formation and the max- 
imum appears at about the kinetic energy cor- 
responding to the maximum transition prob- 
ability. The O* ions which result from the dis- 
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Fic. 7. Construction of the peak shape from the distri- 
bution in initial kinetic energy of O* from NO by method 
A. See caption of Fig. 5 for the designation of the curves. 


sociation of CO exhibit a peak shape of this type 
(Fig. 8). This type occurs when ions of kinetic 
energy less than one volt (approximately) are 
absent. 

A distribution may possess a minimum kinetic 
energy greater than zero for a second reason. If 
the potential curve to which transition takes 
place has a maximum at a value of r greater than 
that existing at transition, the total kinetic 
energy of the products must not be less than the 
height of this maximum above the asymptote. In 
this case there need be no restriction on the 
probability of formation of ions possessing kinetic 
energies near the minimum. If the point on the 
curve, transition to which results in the products 
of minimum energy, lies near the center of the 
Franck-Condon region the probability of transi- 
tion to points near it is relatively large. The 
distribution in kinetic energy of the products 
then cuts off abruptly at the minimum energy. 

A situation in which there is a maximum in the 
potential energy curve but in which the minimum 
kinetic energy is only a few tenths of an electron 
volt may appear as type I or II because of the 
rise of the instrumental efficiency mear zero 
energies. 


Variation of peak height with drawing-out field 


If the ions have very little initial kinetic 
energy, the ion current collected will saturate at 
a small drawing-out field. If, however, the ions 
initially have velocities in random directions, 
more and more ions will be directed to the slit 
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and the current will continue to increase with 
increasing field. Also, a small retarding field will 
reduce the current to zero if the initial kinetic 
energy is zero but ions with initial kinetic energy 
can be collected against a retarding field. Thus 
curves of ion current against the drawing-out 
field will yield some information of the process 
of formation of ions. The height of peaks of 
types O, I and II should not vary with the field. 
The height of peak type III should increase with 
increasing field. Figure 9 shows experimental 
results. 


Appearance potentials and relative probability 
of dissociation processes 


From the construction of curves for types I 
and II, it is evident that the peaks are formed by 
ions which have the least kinetic energy of those 
formed in the dissociation process, even though 
such ions are not the most abundant of those 
formed in the electron beam. If the upper po- 
tential energy curve is of the usual type without 
a maximum, the appearance potential of those 
ions which appear at the peak of the experimental 
curve locates a dissociation limit of the molecule 
since the initial kinetic energy of the ions is zero. 
Dissociation energies obtained from the appear- 
ance potentials of H* from Hz, N* from Ne, N* 


C* From CO 
4 2 ° 2 
VOLTS 
O* From CO 
VoLTs 


Fic. 8. Peak shapes for the C* and O* ions from carbon 
monoxide at 100 volts electron energy. The abscissas give 
ion kinetic energy in volts. Ion accelerating potentia 
increases to the right, the maximum of the C* peak 
corresponding to about 300 volts. It is explained in the text 
how the scale of ion kinetic energy for each peak was 
calibrated. 
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from NO agree with those obtained from retard- 
ing potential methods and from spectroscopic 
data. 

If comparison with other data indicates the 
appearance potential is higher than a dissocia- 
tion limit, it must be concluded that ions of zero 
kinetic energy were not formed. Since peaks of 
type I or II are considered, the probability of 
formation of ions near the minimum energy must 
be fairly large. This situation may be explained 
by assuming a potential energy curve with a 
maximum between the F-C region and the 
dissociation region. Then it is possible to have a 
large number of ions formed with the minimum 
energy which may be some tenths of an electron 
volt above zero. Comparison of data obtained 
in our experiments with those from retarding 
potential methods and spectroscopy aids in de- 
termining the position of the zero on the abscissa 
scale in experimental curves (curves 4, Figs. 
5, 7, 8). The location of the zero of these curves 
is sometimes a difficult matter but a very im- 
portant one in the interpretation of the results. 

The conclusion that the maximum of a peak 
of type I corresponds to ions of zero initial 
kinetic energy if such ions are formed in the 
process has been confirmed in another manner. 
A test was made to see if the maximum of the 
ion peak of N* from Ne appears at the “‘position”’ 
on the scale of accelerating potential for ions of 
zero initial velocity and m/e=14. In order to 
avoid complications due to contact potentials 
and uncertainties of calibration of the voltage 
divider which applies the ion accelerating poten- 
tial to the slit system of the spectrometer the 
following procedure was adopted. Argon and 
nitrogen were admitted into the spectrometer at 
the same time. The voltage interval (about 350 
volts) between the maxima of the N+ and the N+ 
ion peaks was measured. The magnetic field was 
then adjusted so that the A+ peak appeared at 
the position on the voltage scale previously 
occupied by N.* and the interval between A* 
and At* determined. As the ratios of the m/e 
values of A+ and At++ and of N+ and N* are 
equal these intervals should be equal if the 
maximum of the N* peak corresponds to the ion 
of zero initial velocity. It was found that the two 
voltage intervals agreed within 0.2 volt which is 
within the estimated experimental error. 
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Fic. 9. Variation of peak height with the drawing-out 
field in the collision chamber of the mass spectrometer. The 
abscissas give the potential applied across the collision 
chamber between the back plate and the first ion slit (0.3 
cm apart). Curve 1 is for Cot from CO (type 0) and curve 
2 for OF from CO (type III). The curve for a peak of type 
I or II would approximate curve 2 for negative fields and 
curve 1 for positive fields. When the drawing-out field is 
small, stray fields apparently ‘‘reach into” the collision 
chamber causing peculiar focusing effects. For this reason, 
the shapes of the curves for potentials near zero (between 
the vertical dashed lines) have no significance in the 
present discussion. 


Ions of initial velocity near zero are not formed 
in a process which gives rise to an ion peak shape 
of type III. Since ions of appreciable initial 
kinetic energy focus at the maximum of such a 
peak the appearance potential measured at the 
maximum may be considered only as an upper 
limit for the energy level of the asymptote of the 
potential curve. A more accurate determination 
of the position of the dissociation limit can be 
made if it is possible to determine the kinetic 
energy corresponding to some point on the dis- 
tribution. How this was accomplished for the 
single process of this type observed in this work 
is discussed in the section of the paper which 
deals with CO (Process F). 

The study of the shapes of ion peaks indicates 
that neither the height of the peak nor the area 
under the curve is a direct measure of the abun- 
dance of an ion. Instrumental and other factors 
determine the height and area. A peak of type I 
may be ten times as high as one of type III even 
though the total number of ions formed is the 
same in each case. Nevertheless it is possible to 
estimate the relative abundance of the ions 
formed. 


EXPERIMENTAL RESULTS 


In the presentation of the experimental results 
the following symbols will be used : J(X.)—ioniza- 
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ry H* From 


'2 0 8 6 4 2 ° 2 
VOLTS 

Fic. 10. H* ion peak from H, for electrons of 100 volts 
energy. The abscissas give the kinetic energy of the ions in 
electron volts (accelerating potential in volts) measured 
from the position of the maximum. Total accelerating 
potential increases to the right, that corresponding to the 
maximum of the peak being about 300 volts. The peak is a 
superposition of two peaks of types I and III. 


tion potential of X which may be an ion, an 
atom or a molecule. A(X*)—appearance poten- 
tial of X*. This is understood to refer to the 
product of zero initial kinetic energy unless 
otherwise specified. D(X)—heat of dissociation 
of X. EA(X)—electron affinity of X. E(X)—ex- 
citation energy of X. 


Hydrogen 

The restudy of the hydrogen molecule was 
undertaken primarily as a check on the operation 
of the mass spectrometer and the analysis of 
peak shapes. Bleakney?® has demonstrated ex- 
perimentally that two modes of formation of H* 
as a primary ion product from H: under electron 
impact can be distinguished. These entail transi- 
tion to the *2,* and *,* states of the molecule 
ion. By retarding the ions in the collision chamber 
of the spectrometer Bleakney showed that they 
possess initial kinetic energies in the ranges 
predicted by theory. This work was corroborated 
by Lozier® who used the retarding potential 
method. 

From the course of the potential curves one 
should expect H+ peaks of types I and III for 
the ions resulting from transitions to the *2,* 
and 2Z,,+ states of H2*, respectively. The resultant 
peak observed with the mass _ spectrometer 
(Fig. 10) is a superposition of two such peaks. It 
is to be distinguished from a peak of type II 
which may also possess a double maximum 
although only a single potential curve is involved. 

An attempt has been made to construct from 
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the theoretical potential curves the H* peak by 
the method described. Qualitatively the con- 
structed and observed peak shapes agree very 
well. It is impossible to predict the relative 
heights of the two maxima of the peak, however, 
as the relative probability of transition to the 
two states at the same nuclear separation is not 
known. Considering this limitation, it is believed 
that the results of this study indicate that the 
essential features of the peak shape have been 
accounted for. 

Appearance potential measurements at the 
maxima of the H+ and H* peaks yielded the data: 


A(H.*) =15.4+0.1 volts, 
A(H*) =18.0+0.2 volts. 


These data and the values of D(H2) and D(H_"*) 
predicted by them agree within experimental 
error with the results of previous electron impact 
work,” band spectroscopic*®> and thermochem- 
ical*® investigations and wave-mechanical calcu- 
lation.** It is believed that data for other mole- 
cules taken in the same manner are of com- 
parable accuracy when interpreted correctly even 
though knowledge of the structure of the mole- 
cule does not permit so rigorous a check as in the 
case of hydrogen. 


Carbon monoxide”’ 


A summary of data on CO is given in Table I. 
A number of processes in which the ions may be 


TABLE I. Summary of electron impact data for CO. A ppear- 
ance potential in volts. 


A(COt) A(COt*) A(Ct) A(OF) A(O>) 


Hogness, Harkness* 13.9 22.8 24 
(nflass spectrometer) 
Vaughan? 13.9 43.0 22.5+0.2 9.5+1.0 
(mass spectrometer) -+0.2 +1.0 
Tate, Smith, Vaughan’ 27 
(mass spectrometer) 1 
Tate, Smith 14.1 
(total ion current) +0.1 
Lozier$ 20.9+0.1 9.5 40.1 
(retarding potential) 22.8+0.1 20.9+0.1 
present work 14.1 42.0 20.940.2 23.1— 9.5+0.2 
(mass spectrometer) -+0.1 9+0.2 


+0.5 22.840.2 23.5 20. 


*T. R. Hogness and R. W. Harkness, Phys. Rev. 32, 936 (1928). 


26 See H. Beutler and H. O. Jiinger, Zeits. f. Physik 100, 
80 (1935) for a summary of values obtained for the 
dissociation and ionization energies of H2 from spectroscopic 
observations and theoretical calculation. 

26 R. F. Bichowsky and L. C. Copeland, J. Am. Chem. 
Soc. 50, 1315 (1928). 

274 part of this work on CO was presented at the 
Washington Meeting of the American Physical Society, 
April 29, 1939. H. D. Hagstrum and J. T. Tate, Phys. Rev. 
55, 1136A (1939). 
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formed have been identified. These will be 
considered now. 

Process A: CO+e--CO*+2e.—-CO* ions were 
found to appear first at 14.1+0.1 volts. This 
datum agrees very well with the results of previ- 
ous work (see Table I). Further consideration of 
Vaughan’s data by Tate and Smith** gave 14.0 
instead of 13.9 volts. 

Process B: CO+e—-CO*+++3e.—A value of 
A(CO**) =42.0+0.5 volts has been obtained. It 
is in substantial agreement with Vaughan’s value 
of 43.0+1.0 volts. 

Process C: CO+e—C+0O-.—The onset poten- 
tial of the resonance process yielding O~ ions was 
found to be 9.5+0.2 volts (Fig. 11), in good 
agreement with the results of previous investiga- 
tions. The fact that this potential measured with 
the ion accelerating potential of the mass spec- 
trometer adjusted to the maximum of the peak 
agrees with the extrapolated value of Lozier for 
ions of zero initial kinetic energy indicates that 
such ions are produced in the process in ap- 
preciable abundance. This result confirms Lo- 
zier’s observation that O~ ions of kinetic energy 
from zero to approximately one volt are formed. 
The potential curve to which transition takes 
place must lie very near the asymptote in the 
Franck-Condon region and perhaps crosses it 
there. This would indicate that the curve 
possesses a minimum and Lozier suggests that a 
careful search with a mass spectrometer might 
reveal CO~ ions. A search was made for such ions 
but none could be detected. If they are formed 
their abundance cannot be greater than 0.01 
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ELECTRON ENERGY IN VOLTS 
Fic. 11. Efficiency of formation of O~ ions from carbon 
monoxide as a function of the energy of the impacting 
electrons. The ion accelerating potential was adjusted to 
the maximum of the peak. 


*8 J. T. Tate and P. T. Smith, Phys. Rev. 39, 270 (1932). 
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A[A‘] 


Fic. 12. Typical 
appearance poten- 
tial curves for ions » 
from carbon mon- 
oxide. The C* curve 
was obtained with 
the ion accelerating 
potential held con- 
stant at the value 
corresponding to 
the maximum of the | 


peak for higher elec- 
tron energies. | 


20 2 22 2 24 
ELECTRON ENERGY IN VOLTS 


percent of the CO* ion. Massey*® suggests that in 
low pressure experiments of this type it will not 
be possible to observe the formation of negative 
molecular ions because they will dissociate by the 
reverse process before making sufficient collisions 
to get rid of their excess vibrational energy. 

Process D: CO+e—-Ct+0O-+e.—The first ap- 
pearance potential of C* ions at 20.9 volts 
(Fig. 12) coincides within experimental error with 
the second appearance potential of O~ ions. (Fig. 
11). This potential agrees very well with Lozier’s 
value. 

Lozier demonstrated quite conclusively that 
the numbers of positive and negative ions formed 
for electron energies between 20.9 and 22.8 volts 
electron energy are equal. This has been verified 
within experimental error in the present work 
for the C+ and O- ions. Any doubt as to the 
identification of the ions involved has thus been 
removed. 

Process E: CO+e—-C*++0-+2e.—That this 
process is much more probable than process D 
is seen from the rates at which the C+ ion current 
rises with electron energy below and above 22.8 
volts. Thus the C* peak observed at 100 volts 
electron energy is made up in the main of ions 
from this process. The shape of the peak is shown 
ir Fig. 8. It is seen to be of type I. The good 
agreement of the position of the second break at 
22.8 volts in the curve for C* against clectron 


Py H.S. W. Massey, Negative Jons (Cambridge, 1938), p. 
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ELECTRON ENERGY IN VOLTS 
Fic. 13. One set of appearance potential curves for O* 
ions of 1.5, 2.0, 2.5 and 3.0 volts initial kinetic energy. It 
was found that no better data could be obtained by at- 
tempting to take readings at smaller voltage intervals 
because of the extremely low ion current and the fact that 
the position on the peak could be checked only before and 
after the data for each curve were obtained. 


energy with Lozier’s extrapolated value indicates 
that ions of zero initial velocity are formed in the 
process and that for this reason the potential 
curve involved cannot possess a maximum. 
Vaughan’s conclusion that the C* ions possess 
kinetic energy is untenable if it is to be inter- 
preted as meaning that ions of zero initial kinetic 
energy are not formed. 

Process F: ions from 
CO have been observed by previous investigators 
but have not been studied in any detail. The only 
quantitative data on them were given by Tate, 
Smith and Vaughan’ who report the peak to be 
10 percent as high as that of Ct at 100 volts 
electron energy and the appearance potential as 
27+1 volts. It thus appears to be of interest and 
importance to ascertain the reason for the low 
peak height and to attempt to fit energy deter- 
minations into the interpretation of other 
processes. 

The dissociation process from which the O* 
ions result could be either CO+e—C~-+Ot+e or 
CO+e—-C+0++2e. Failure to detect C~ ions 
eliminates the first of these possibilities. 

The O+ peak is plotted below that of C+ from 
CO to the observed relative abundance in Fig. 8. 
The shapes of the two peaks are entirely differ- 
ent, that of O*+ being of type III. It is evident 
that the low peak height for O* ions in com- 
parison with that for C* ions is due to the nature 


AND jJ. T. TATE 


ao 


KINETIC ENERGY IN VOLTS 


Fic. 14. Relation of the appearance potential and kinetic 
energy of ions from the dissociation of CO. Curve 1 is for 
C* and curve 2 for O*. The lines are drawn with the values 
of slope predicted by the laws of conservation of energy and 
momentum applied to the dissociation process (see reference 
cited in footnote 7). Two sets of data are plotted for O*. 
In placing curve 2 the data for lower kinetic energies were 
favored since they were obtained at the higher intensity 
portions of the peak (see Fig. 8). 


of the distribution in kinetic energy of the 
products rather than to any significant difference 
in the probability of occurrence of the processes 
from which they result. The parts of the peaks 
corresponding to initial kinetic energies of more 
than 2 volts are about equally high. Thus the 
assumption that the O* ion is not present in any 
appreciable abundance because the total peak 
height is small is erroneous. 

The location of the dissociation limit for this 
process requires a different technique from that 
used for process E since ions which focus at the 
maximum of the peak possess an unknown 
amount of initial kinetic energy. The scale of 
kinetic energy was calibrated by using as the 
zero point the position of the maximum of the 
small O* ‘impurity peak”’ which remained when 
the electron energy was reduced to a value below 
the appearance potential of any ions from CO. 
The shape of this peak and the fact that the ions 
result in the main from the dissociation of water 
vapor impurity are taken as good evidence that 
the peak is composed of ions possessing almost 
no initial velocity. The appearance potentials 
(Fig. 13) of O* ions of 1.5, 2.0, 2.5 and 3.0 volts 
initial kinetic energy were measured by ad- 
justing the ion accelerating potential to the 
proper positions on the peak of Fig. 8. 

Extrapolation of the relation of appearance 
potential and kinetic energy to zero kinetic 
energy (Fig. 14) yields 23.5 volts for the level of 
the dissociation limit above the ground state of 
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the molecule. An interesting comparison may be 
made here with the datum of Tate, Smith and 
Vaughan. Their appearance potential, presum- 
ably measured on the maximum of the peak, 
would, according to the lower curve of Fig. 8, be 
that of ions of approximately 1.5 volts initial 
kinetic energy. From Fig. 13 it is seen that this 
ion should be formed first at about 27 volts 
which is just the value given by Tate, Smith and 
Vaughan. 

The low intensity of the O* ion current and 
the limitations of the method of calibration of 
the scale of kinetic energy make it necessary to 
consider the value of 23.5 volts as an upper 
limit. Any possible correction can be no more 
than a few tenths of a volt, however, and it is 
estimated that the true dissociation limit lies no 
lower than 0.4 volt below 23.5 volts. 


Heat of Dissociation of CO-(D(CO)) 


The energy relations set up for the processes 
in which O*, O~ and C+ are formed limit the 
values which can be assigned to D(CO). The 
largest value of D(CO) compatible with the 
dissociation limit of process F lies in the range 
9.5 to 9.6 electron volts. Interpretation® of 
various processes which led to D(CO) = 11.6 elec- 
tron volts seems untenable. The value 9.6 volts 
remains as the only value of D(CO) which yields 
a consistent interpretation of all the electron 
impact processes observed in CO. 

We have, for example, in process D 


A(C*) = D(CO)+1(C) +E(C*) 
+EA(O-)+E(O-), 


For process E we have 


TABLE IT. Summary of electron impact data for NO. A ppear- 
ance in volts. 


A(NO*) A(NO**) A(N*) A(O*) A(O>) 
Hogness, Lunn* 9 22 21 


(mass spectrometer) 

Tate, Smith™ 9.5 
(total i ion current) +0.1 

Tate, Smith, Vaughan* 9.5 41 22+0.5 21 
(mass spectrometer) +0.1 +1 

Hanson® 21.7 +0.2 3.2 X0.5 
(retarding potential) 

Present work 9.5 41.3 19.9+0.2 20.5 7.0+0.3 
(mass spectrometer) +0.1 +0.3 21.840.2 240.2 20.0+0.3 


*T. R. es and E. HL. Lunn, Phys. Rev. 30, 26 (1927). 


A[At] A[(NO*] 


ELECTRON ENERGY IN VOLTS 


A(O*] 


l 


ELECTRON ENERGY IN VOLTS 


Fic. 15. Typical appearance potential curves for ions 
from nitric oxide. The slight slope of the O* curve below 
the break is due to amplifier drift and a small amount of 
ion current from dissociation of water vapor impurity. 


A(Ct) =D(CO)+1(C)+ 
22.8=9.6+-11.3+0+ 1.96, 
For process C we have 


+ E(O-), 
9.5=9.6+0—2.2+2.2. 


In process F the relations are 


A(O*) = D(CO) + E(C) + E(0*), 
23.1=9.6+13.5+0+0. 


The assignments of E(C)=E(Ct)=0 and 
EA(O-) = E(O-) =2.2+0.2 seem reasonable from 
spectroscopic data and earlier data on impact 
phenomena. It should be noted that the value 
D(CO) =9.6 volts depends on the appearance 
potentials in processes C, D and E which are not 
subject to the same error involved in A(O*). 
They were determined from peaks of a different 
type. 

The value D(CO)=9.6 volts is not in agree- 
ment with either of two values proposed by band 
spectroscopists, however. Schmid and Geré*° 
interpret the three predissociation limits at 9.61, 


*® R. Schmid and L. Gerd, Zeits. {. Physik 99, 281 (1936). 
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11.11 and 11.60 volts as corresponding to three 
dissociation limits of the molecule leading di- 
rectly to D(CO)=6.921 volts.*! Schmid*® has 
attempted to reconcile the electron impact data 
of Lozier with this value of D(CO) by assuming 
tetravalent carbon products. He calculates the 
level of C+(*#P)+O(P) as 23.5 volts above the 
ground state. This is outside the experimental 
error for A(C+) at 22.8 volts. Furthermore if 
A(C*) is in error it would in all probability be 
higher than the actual dissociation limit. 
Herzberg,* on the other hand, considers only 
the predissication at 11.11 volts as corresponding 
to a dissociation limit of the molecule. The pre- 
dissociation at 11.60 volts is considered to be 
only an upper limit for the level of some asymp- 
tote and that at 9.61 volts as accidental pre- 
dissociation. Herzberg’s final conclusion is that 
D(CO) is most probably 9.144 volts* with the 
limit at 11.11 volts corresponding to C(*P) 
+O('D). It is possible to interpret only process 
C of the dissociation processes observed in CO 
under electron impact on this basis, however. 
This can be done by assuming unexcited O- and 
C in the !S state. If D(CO)=9.144 volts it is 
necessary to assume an experimental error in 
the appearance potentials for processes D and E 
of about 0.5 volt. This conclusion seems a little 
severe on the basis of the agreement of inde- 


4 
O” From NO 
2 
° 
5 10 25 30 


15 20 
ELECTRON ENERGY IN VOLTS 


Fic. 16. Efficiency of formation of O~ ions from nitric 
oxide as a function of the energy of the impacting electrons. 
The ion accelerating potential was adjusted to the maxi- 
mum of the ion peak at each reading. 


| The new conversion factor of 8066.0 cm~ per electron 
volt has been used in this paper. For this reason a number 
of the values cited differ slightly from those given in the 
original papers. See Table 2 in G. Herzberg, Molecular 
Spectra and Molecular Structure, I. Diatomic Molecules 
(Prentice-Hall, New York, 1939). 

® R. Schmid, Zeits. f. Physik 99, 274 (1936). 

% G. Herzberg, Chem. Rev. 20, 145 (1937). 
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TATE 


N* From NO 


POSITIVE CURRENT 


o 8 6 4 
voLTs 

Fic. 17. Peak shapes for N* and O* ions from nitric 
oxide. Position 1 of the O* peak relative to the scale of 
kinetic energy corresponds to the construction of the peak 
by method A, position 2 to method B (see text and Fig. 7). 
The maximum of the N* peak corresponds to an acceler- 
ating potential of about 300 volts. 


pendent electron impact experiments in He, Ne 
and NO by different methods with each other 
and with band spectroscopic data which is much 
more conclusive than that for CO. 

All thermochemical data on the heat of sub- 
limation of carbon, L(C), published prior to 
Jan. 1, 1934 have been reviewed by Bichowsky 
and Rossini** who conclude that the value most 
probably lies in the range 110 to 170 kcal. 
D(CO) =9.6 volts (221.3 kcal.) yields L(C) =5.9 
volts (136.3 keal.). 

Goldfinger and Jeunehomme* have calculated 
vapor pressure curves for C from three values for 
L(C) (170.2, 124.8 and 108.3 kcal.)*! which are 
in agreement with the predissociation limit in 
CO at 11.11 volts. It was found that the direct 
measurements of the vapor pressure of C de- 
scribed in the literature lie mostly between the 
curves calculated for L(C) = 170.2 kcal. and 124.8 
kcal. Although Z(C) = 170.2 kcal. gives the better 
agreement, it was eliminated for spectroscopic 
reasons and L(C) =124.8 kcal. corresponding to 
D(CO) =9.144 volts taken as the most probable 
value. It should be pointed out that the value 


4 R. F. Bichowsky and F. D. Rossini, The Thermochemis- 
try of Chemical Substances (Reinhold, New York, 1936). 

36 P. Goldfinger and W. Jeunehomme, Trans. Faraday 
Soc. 32, 1591 (1936). 
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60 
ELECTRON ENERGY IN VOLTS 


Fic. 18. Efficiency of formation of ions from Ne. Curve 1 
is for Ne* ions and curve 2 for ions of m/e=14. A; indicates 
the onset potential of process C and A, that of process B. 
The ordinates of each curve have been arbitrarily chosen 
and have no relative significance. 


L(C) = 136.3 kcal. based on D(CO) =9.6 volts is 
perhaps in better agreement with vapor pressure 
measurements as it is somewhat higher than the 
“spectroscopic”’ value chosen by Goldfinger and 
Jeunehomme. 


Nitric oxide*® 


Of particular interest in studying NO is the 
process NO-N*++QO- and the dissociation 
process in which O* is formed. The results of this 
paper and of earlier work are given in Table II. 
A number of processes have been identified. 

Process A: NO+e—NO*++2e.—The value 
I(NO) (Fig. 15) agrees very well with previous 
measurements. 

Process B: NO+e—-NO** + 3e.— A (NO**) was 
found to be 41.3+0.3 volts. 

Process C: NO+e—-N+0O-.—The efficiency of 
formation of O~ ions as a function of electron 
energy is shown in Fig. 16. The onset potential of 
the resonance process measured at the maximum 
of the ion peak is 7.0+0.3 volts. The fact that 
this potential is higher than that obtained by 
extrapolation by Hanson® indicates that ions of 
zero initial kinetic energy are not formed in the 
process. 

Process D: NO+e—N*++0O-+e.—This process 
has been identified for the first time in this work. 
The value of A(N*) is found at 19.9+0.2 volts 


** A part of this work on NO was presented at the New 
York Meeting of the American Physical Society February 
24, 1940. H. D. Hagstrum and J. T. Tate, Phys. Rev. 57, 
561A (1940). 


(Fig. 15). 20.0+0.3 volts was obtained for 
A(O-). The appearance potential A(N*) was 
measured with the same ion accelerating field as 
was used to obtain the N* peak in process / and 
it is believed that the ions (N*) had zero kinetic 
energy here as they did in process E. 

The onset potential of this process is very 
nearly equal to the electron energy calculated 
from D(NO)* and J(N) at which the unexcited 
products N* and O should appear. (19.84 
volts).** *! If Lozier’s value (2.2 volts) is assumed 
for EA(O)* it is necessary that E(N*+)+E(O-) 
equal approximately 2.2 volts. The nearest state 
in N* is the 'D at 1.90 volts. Identification of this 
value with the 2.2-volt excitation is a possibility 
but better numerical agreement is obtained if the 
analogy with the process CO+e--Ct+O--+e is 
carried out and the excitation assigned to the O~ 
ion. This excited state of O~, lying very near the 
positive energy continuum, is the only one which 
the theory of the negative ion renders at all 
probable.** 

Process E: NO+e—-N*t+0+2e.—The 
peak shown in Fig. 17 is to be attributed mainly 
to ions from this process. The shape of the peak 
is of type I. The appearance potential of 21.8+0.2 
volts (second break in the C* curve of Fig. 15) 
measured at the maximum of the peak is in 
better agreement with Hanson's extrapolated 


A(Net*) A(N*) A(Net*) A(N**) 

Vaughan? 15.8+0.1 24.54-0.1 40 +1(?) 

(mass spectrometer) 47+1 
Tate, Smith® 15.7 

(total ion current) 
Lozier’ 24.27+0.1 

(retarding potential) 
Tate, Smith, Vaughan’ 15.65+0.1 

(mass spectrometer) 
Present work 15.7+0.1 24.340.2 49.540.5 6442 


(mass spectrometer) 


37 D(NO) =5.29 volts calculated from D(N2)=7.384 
volts, D(O2) =5.082 volts and the heat of formation of NO 
from thermochemistry. See G. Herzberg, Molecular Spectra, 
etc., Table 36 and R. F. Bichowsky and F. DL. Rossini, 
reference 34. 

38 Atomic excitation and ionization energies have been 
taken from R. F. Bacher and S. Goudsmit, Atomic Energy 
States (McGraw-Hill, New York, 1932). 

39 For a discussion of the possibility of excited states in 
negative ions see O. Oldenberg, Phys. Rev. 43, 534 (1933) 
(especially footnote 9 in this reference); W. W. Lozier, 
Phys. Rev. 46, 268 (1934) and H. S. W. Massey, reference 
29, p. 19. 
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value of 21.7+0.2 volts than previous mass 
spectrometric data. This agreement between the 
data of the two electron impact methods is again 
taken as evidence that ions of zero initial velocity 
are formed in the process. The N* product is 
unexcited whereas the O atomic product is in its 
'D state. The energy difference between the two 
breaks in the N* curve of Fig. 15 (1.9 volts) is 
thus to be interpreted as the excitation of the 'D 
state rather than the electron affinity of the 
oxygen atom. 

Process F: NO+e—N+0O*+2e.—The O* peak 
is shown in Fig. 17 drawn to the same scale used 
in the curve for N*. The peak is assumed to be of 
type II. In curve 1 the ions which focus at the 
maximum are assumed to have no kinetic energy. 
But the appearance potential of 20.5+0.2 volts 
(Fig. 15) is not compatible with the known value 
of D(NO) and excitation energies. One must 
assume then that the dissociation products are 
unexcited but have kinetic energy. The total 
kinetic energy turns out to be 20.5—18.9 or 1.6 
electron volts of which the O* ion possesses 0.75 
electron volt. The correct position of the O* peak 
is then that shown in curve 2 in which the 
maximum is at 0.75 volt. 

The experimental curve 2 can be most satis- 
factorily constructed from the curve showing the 
distribution in kinetic energy of the ions at the 
instant of formation by assuming that the upper 
potential energy curve has a maximum of 1.6 
volts. 

Herzberg and Mundie*® have assumed a po- 
tential maximum of 1.3 volts to explain a pre- 
dissociation limit. It appears that the possible 


TABLE IV. Summary of electron impact data for Oz. A ppear- 
ance potential in volts. 


A(O2*) A(O2t*) A(O*) A(O7) 
Hogness, Lunn* 13 20 
(mass spectrometer) 
Kallmann, Rosen** 13 19.5 
(mass spectrometer) 
Tate, Smith 12.5 6 
(total ion current) 
Lozier® 2.9+0.2 
(retarding potential) 12.0+0.2 
Present work 12.3+0.1 50.040.5 19.2+0.2 
9+0. 


(mass spectrometer) 


* T. R. Hogness and E. G. Lunn, Phys. Rev. 27, 732 (1926). 
** H. Kallmann and B. Rosen, Zeits. f. Physik 61, 61 (1930). 


 G. Herzberg and L. Mundie, J. Chem. Phys. 8, 263 
(1940). 
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AND J. T. TATE 
existence of potential curves with maxima as high 
as is necessary to postulate here is not in doubt. 


Nitrogen 


Data from previous electron impact studies in 
Ne are given in Table III. Restudy of this mole- 
cule with a mass spectrometer was undertaken to 
extend the experimental study of peak shapes and 
to obtain, if possible, a definite interpretation of 
the increases in the efficiency of production of 
ions of m/e=14 observed by Vaughan.? 

Process A: Nz+e—-N.++2e—-The value of 
15.7+0.1 volts obtained for J(N2) is in agree- 
ment with previous electron impact data listed 
in Table III and with the value of 15.581 volts 
determined from the limit of a Rydberg series of 
bands in the molecular spectrum.‘ *! 

Process B: No+e—N2**-+3e.-The efficiency 
of formation of ions of m/e=14 as a function of 
electron energy is shown in Fig. 18. The second 
break at 49.5+0.5 volts does not agree very well 
with the position of either break reported by 
Vaughan? but is undoubtedly to be identified 
with the more pronounced break in his curve at 
47 volts. It is believed that the first break in 
Vaughan’s curve at 40 volts may be the onset of 
the formation of CO** ions. CO was an acknowl- 
edged impurity in the apparatus. Furthermore 
the energy difference between the breaks (7 
volts) is just the difference found in this work 
between A(CO**) and the second break in curve 
2 of Fig. 18. No definite identification of the 
ionization or dissociation process responsible for 
this increase in the ion current has as yet been 
made. It will be indicated in the discussion of 
process C that the break of 49.5 volts may now 
be ascribed to the onset of double ionization of 
the molecule. Such a process is to be expected 
from the results for other molecules. 

Process C: No+te—-N*++N+2e.—The peak at 
m/e=14 for 100-volt electrons is of type I and 
is very similar in appearance to the peaks for C* 
from CO and N+ from NO. (see Figs. 8 and 17). 
The appearance potential measured at the max- 
imum of the peak is 24.3+0.2 volts. Agreement 
with Lozier’s extrapolated value of 24.27 volts‘ 


41 R.S. Mulliken, Phys. Rev. 46, 144 (1934); A. van der 
Ziel, Physica 4, 373 (1937); R. E. Worley and F. A. 
Jenkins, Phys. Rev. 54, 305 (1938). 
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indicates that ions of zero initial velocity are 
formed. 

The ion peak at 100 volts electron energy is 
also made up of ions from the process which 
begins at 49.5 volts. These ions may be either 
N.** or N* from a second mode of dissociation of 
the molecule. The efficiency of formation of ions 
of 2 volts initial kinetic energy as a function of 
electron energy was determined by setting the 
ion accelerating potential at a value 2 volts less 
than that corresponding to the maximum of the 
peak. The curve shows no second break at any 
point and resembles curve 1 of Fig. 18. This result 
indicates that the process starting at 49.5 volts 
is double ionization of the molecule. If the ions 
were N** a second break should have been ob- 
served at about 53.5 volts electron energy in the 
efficiency curve for 2 volt ions. 

Process D: ex- 
tremely small ion peak was observed at the ion 
accelerating potential corresponding to m/e=7. 
The abundance of the ion relative to N»* is only 
0.006 percent. These ions could either be N+*+ 
resulting from dissociation of the molecule or 
N.t*+*, Since no N2*** ions were found at the 
m/e=9.35 position it is likely that the ions of 
m/e=7 are N**. Furthermore the shape of the 
ion peak appears to be of type III indicating that 
the ions possess considerable initial kinetic 
energy. This could be possible only if the ions are 
products of dissociation rather than multiple 
ionization. The appearance potential of the ions 
at the maximum of the distribution was found to 
be 64+2 volts. Lack of knowledge of the kinetic 
energy of the ions precludes any further interpre- 
tation of the process. 


Oxygen"” 


The work on oxygen revealed the formation of 
O,** and O** and gave additional information 
on processes in which other ions were formed. 
The results of this and other recent experiments 
are shown in Table IV. 

Process A: value of 
12.3+0.1 volts was obtained for J(Q,). This is 
lower than the value of 12.5 volts obtained by 


_ © A part of this work on O2 was presented at the Wash- 
ington Meeting of the American Physical Society, April 26, 
1940. J. T. Tate and H. D. Hagstrum, Phys. Rev. 57, 
1071A (1940). 
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Fic. 19. Efficiency of formation of O~ ions from oxygen. 
The spectrometer was adjusted to collect ions of zero initial 
kinetic energy. 


Tate and Smith*® but is in better agreement 
with the value of 12.2 volts obtained by means 
of the cycle = D(Oz)+1(O)—D(O.*) and 
spectroscopic determinations of D(Q.), J(O) 
and D(O-*).* 

Process B: relative 
efficiency of formation of ions of m/e=16 from 
Oz by electrons of energies up to 100 electron 
volts has been measured. The curve is similar to 
that for ions of m/e=14 from N» shown in Fig. 18 
(curve 2). In this case, however, the second break 
occurs at 50.0+0.5 volts. It has been demon- 
strated by the method used for process B in N» 
that this break indicates the onset of the pro- 
duction of O.**. 

Process C: O2+e—-O0+0-.—The efficiency of 
formation of O- ions as a function of electron 
energy is shown in Fig. 19. The first appearance 
potential at 3.0+0.4 volts is undoubtedly to be 
identified with the value of 2.9 volts obtained by 
extrapolation by Lozier*® for ions of zero initial 
kinetic energy. At first sight this result appears to 
contradict Lozier’s observation that ions of 
kinetic energy less than 1.3 volts are not formed. 
It is possible, however, that the ions observed in 
this work and those observed by Lozier result 
from transitions to different potential curves 
having the same asymptote. Massey“ has 
interpreted the resonance process observed by 
Lozier as transition to the unstable 2, state of 
O.~. It is suggested that the process observed in 


43 R. S. Mulliken and D. S. Stevens, Phys. Rev. 44, 720 


(1933). 
“H.S. W. Massey, reference 29, pp. 28 and 52. 
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this work with a mass spectrometer involves 
transition to the II, state. If the potential curve 
for this state has the position proposed by 
Massey ions of zero initial kinetic energy could 
be formed. Furthermore if the ions from this 
process possess initial kinetic energies no greater 
than a fraction of a volt they would not be ob- 
served by the retarding potential method. 

Process D: O2+e—-0++0+2e.—The shape of 
the O* ion peak at 100 volts electron energy is of 
type I. Ions of kinetic energy near zero are rela- 
tively numerous. The appearance potential 
measured at the maximum of the peak (19.2+0.2 
volts), however, is 0.5 volt greater than the po- 
tential at which ions of zero initial velocity should 
appear. [D(O2)+J(O) =5.08+ 13.62=18.7 volts ]. 
On the basis of the accuracy of similar measure- 
ments in other molecules this difference is con- 
sidered to be well outside the experimental error. 
It cannot be attributed to excitation of the 
products, as neither O nor O* possesses an ap- 
propriate state. The appearance potential must 
then correspond to ions of 0.25 volt initial kinetic 
energy. The peak shape analysis is compatible 
with a potential curve which has a maximum 0.5 
volt above the asymptote (see the discussion of 
process F in NO). 

Process E: O2+e—-0++0-+e.—The second 
appearance potential of O- ions from Oz: is 
18.9+0.4 volts. (Fig. 19). This value is 0.2 volt 
higher than the potential energy of the products 
O++0O above the ground state of O2. There can, 
however, be little doubt that two values should 
be identified. The accuracy of the datum for this 
process is not to be compared with that for 
process D which is about 20 times as probable. 

Since it is not possible to identify the appear- 
ance potential with any other dissociation limit 
of the molecule, the products Ot+O- are excited 
by very nearly the amount EA(O)=2.2 volts. 
Since the lowest state of O* is the *D° at 3.3 volts 
the O- ion must assume all of this energy. This 
is an additional example of the formation of the 
O- ion in the ‘‘loosely bound” state.*® 


AND J. T. TATE 


’ The curve of O* ion current as a function of 
electron energy exhibits only the single break at 
19.2 volts which has been identified with process 
D. Since O- ions from process E are observed at 
18.9 volts one questions why O* ions which are 
formed in equal numbers in the process are not 
observed at the same potential. This is to be 
explained as follows. The measurement of A (O*) 
was made at the maximum of the O* peak which 
has been shown to correspond to ions of 0.25 volt 
kinetic energy. When the electron energy was 
varied the ion accelerating potential was not 
changed so that ions of 0.25 volt kinetic energy 
were being collected when the onset potential of 
process E was measured. Since the dissociation 
limits for O*(4S°)+O(@P) and O*(¢S°)+O- (ex- 
cited by 2.2 volts) almost coincide, only one 
appearance potential should be observed for O+ 
ions of any initial velocity. 

Process F:  ex- 
tremely small ion peak was observed at the 
m/e=8 position. Its abundance was of the order 
of 0.01 percent of O2*. Failure to find any ions at 
the position corresponding to O,*** indicates 
that the ions of m/e=8 are more likely O+* than 
O.*+*++, No appearance potential or peak shape 
measurements were attempted. Because of the 
effect of O2 on the filament of the spectrometer 
the gas pressure and the electron beam intensity 
could not be run high enough to make such 
measurements possible. 
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Experimental Study of the Diffuse Scattering of X-Rays 
by Potassium Chloride Crystals 


STANLEY SIEGEL 
Ryerson Physical Laboratory, University of Chicago, Chicago, Illinois 
(Received December 23, 1940) 


Diffusely scattered radiation has been studied with a KCl crystal and copper radiation. 
Observations have been made in the immediate region of a strong reflection. The shapes and 
half-widths of the maxima have been studied as functions of the direction of incidence. Observa- 


tions have been found to agree with the theory. 


GENERAL 


T the December, 1939, meeting of the 

American Physical Society in Chicago, W. 
H. Zachariasen proposed a new theory of the 
diffuse scattering of x-rays by crystals. In one 
aspect of the theory, it was predicted that the 
diffusely scattered x-radiation would appear as 
sharp diffraction maxima. In a subsequent paper 
at the same meeting, S. Siegel and W. H. 
Zachariasen reported an experimental study of 
these maxima. The results of the theory' and 
experiment have been published. 

The difference between this theory and the 
earlier work of Debye, as well as the experimental 
studies by Jauncey and his collaborators, has 
already been treated.! 

The existence of streaks in Laue photographs 
has been known for many years, and in 1923, 
Faxén* suggested that these streaks might be 
caused by a thermal effect in the crystal. 
Wadlund‘ observed very prominent streaks in 
some photographs and Zachariasen® interpreted 
these as being due to a two-dimensional lattice 
effect. Recently Preston® has shown the effect of 
temperature on these streaks, thus confirming 
Faxén’s suggestion. Preston observed streaks 
when general radiation was used ; however, when 
monochromatic radiation was used, maxima, 
similar to Laue spots, appeared. Preston inter- 
preted these on the basis that the thermal vibra- 
tions of the lattice break the crystal into groups 
composed of an atom and its nearest neighbors, 


!W. H. Zachariasen, Phys. Rev. 57, 597 (1940). 
an Siegel and W. H. Zachariasen, Phys. Rev. 57, 795 


3H. Faxén, Zeits. f. Physik 17, 266 (1923). 

4A. P. R. Wadlund, Phys. Rev. 53, 843 (1938). 

5 W. H. Zachariasen, Phys. Rev. 53, 844 (1938). 
°G. D. Preston, Proc. Roy. Soc. 172, 116 (1939). 


and these thermal motions change the interatomic 
distance from one group to another, thus pro- 
ducing new maxima. Raman, apparently unaware 
of Preston’s work, or of the work carried out at 
Ryerson Laboratory, published photographs 
showing these spurious maxima.’ He interpreted 
the maxima as a specular reflection of x-rays 
resulting from stratifications of density in the 
crystal which are not static but dynamic in 
character. Other articles have since appeared. 
Knaggs and collaborators® have published photo- 
graphs of these maxima. Lonsdale and collabo- 
rators® have made a detailed qualitative study of 
this effect. The conclusions reached by them are 
those predicted by theory, many of the results 
however, having been already studied quanti- 
tatively. More recently, another paper by Raman 
and Nilakantan has appeared.'° 

Although the agreement between theory and 
observations in our original experiment was quite 
good, there was some discrepancy between theo- 
retical and observed widths of the diffuse 
maxima. Thus it was decided to carry out the 
observations in more detail, to extend them as 
close as possible to the Bragg angle, and to in- 
crease the accuracy of observations with better 
equipment. The original experiment was _ per- 
formed with a crystal of sodium chloride. Since 
the theory was developed for crystals containing 
atoms of only one kind, a good approximation to 
such a case would be potassium chloride, which 
was used throughout the present experiment. 

7C. V. Raman and P. Nilakantan, Nature 145, 667 
(1940); Proc. Ind. Acad. Sci. 11, 379 (1940). 

8]. E. Knaggs, K. Lonsdale, A. Miiller, A. R. Ubbelohde, 
Nature 145, 820 (1940). 
(1940) Lonsdale, I. E. Knaggs, H. Smith, Nature 146, 332 


10 C, V. Raman and P. Nilakantan, Proc. Ind. Acad. Sci. 
12, 141 (1940). 


371 


at 
‘SS 
re 
ot 
ch 
It 
as 
ot 
ry 
of 
on 
X- 
)+ 
X- 
1e 
er 
ait 
eS 
in 
or 
‘h 
d 
| 
It | 
e 
| 
|_| 


372 STANLEY SIEGEL 


4 


Fic. 1. A section 
of the reciprocal 
lattice. 


METHODS FOR OBSERVING THE MAXIMA 


As the theory indicates, the streaks are pro- 
duced by the continuous radiation in the x-ray 
beam while the spots are individual diffraction 
maxima due to the characteristic radiation. In 
the present paper only the diffuse spots will be 
dealt with. These maxima are observed whenever 
Bragg’s law is nearly satisfied. 

Figure 1 shows the conditions under which they 
may be observed. This figure and the following 
description have been given in Zachariasen’s 
paper.! Vectors ky and k are the directions of 
incidence and scattering of the x-ray beam. 
|ko| = |k| =1/A. By is a vector in the reciprocal 
lattice, and its shortest separation from the 
vector k—Kp is to. Let k vary. Its terminus will 
describe a sphere of radius 1/\ about the terminus 
of —k», thus producing the sphere of reflection. 
The points on the sphere of reflection which are 
closer to B,, than to any other lattice point define 
a sector of the solid angle 47, so that the same 
vector By is associated with all directions k 
within the sector. As k varies within a sector, 7 
varies with it since k—k)+*9=By. Considering 
ky and B, to be fixed, to becomes a minimum 
value, tmin When k is parallel to By+ko. Thus 
Tmin — |Buat+ko| |k| or AT min —A sin 26z, 
where 6, is the Bragg angle and 0;=6,+A, A 
being a small deviation from the Bragg angle. 
The intensity of the scattering reaches a maxi- 
mum whenever 7o is a minimum, and the 
direction of the intensity maximum, k,, is 
[1+Atmin =By+Ko. The scattering angle cor- 
responding to the maximum is shown to be 
20m ~ 20+2A sin? 6g. For scattering directions k 
near a maximum k,,, (Aro)? = (Armin)? ++ x? where x 


is the angle between k and k,,. The intensity 
distribution is symmetrical about the maximum, 


and the half-width at half-maximum is x, 


=X|Tmin| = |A| sin 26z. 

Thus, in a transmission photograph there will 
always be some reflection planes for which 7,,i, 
will give appreciable intensity, enough to be 
conveniently recorded. 

Transmission photographs were obtained by 
Preston and later by Raman. They obtained 
Laue photographs using apparently large pin- 
holes to define the primary beam, and employing 
sensitizing screens to hasten recording. But for 
quantitative work it is necessary to use a narrow 
pinhole or slit to keep the divergence of the beam 
at a minimum and to avoid the use of sensitizing 
screens, for such screens widen the spots, causing 
them to become too diffuse for any width or 
intensity measurements. The author tried to 
obtain transmission photographs using a small 
pinhole, monochromatic copper radiation, and a 
thin KCI crystal. One would expect weak reflec- 
tions with Cu Ka and KCI, except at a scattering 
angle of 90°, in which case the absorption would 
become quite large. After a 36-hour exposure only 
very faint spots were discernible. The transmis- 
sion method does not appear very amenable to a 
quantitative study. 

The best method for observation is obtained by 
reflecting the primary beam from the crystal 
surface and making observations on a strong 
reflection. In this case the crystal can be set to 
produce any desired 7,,in values, and, depending 
on the nature of the crystal surface, observations 
can be made on diffuse maxima having relatively 
large intensity values. It must be remembered, 
however, that in any method of observation, the 
divergence of the primary beam must be smaller 
than the widths of any of the observed diffuse 
maxima. This immediately rules out the use of 
wide slits or sensitizing screens, and gives rise to 
small recording intensities. This point will be 
elaborated. 


APPARATUS AND EXPERIMENT 


A Bragg spectrometer was used in making 
observations. All such observations could be 
duplicated to within 1’ of arc, so that any error 
was less than that value. The beam was defined 
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by two slits separated a distance of 8.0 cm. The 
width of each slit was 0.005 cm, the height being 
(0.3 cm. The distance between the second slit and 
crystal was 5.0 cm, and between crystal and 
photographic plate, 6.04 cm. A copper x-ray tube 
was used throughout the experiment and the 
radiation was filtered with a nickel foil. Special 
care was taken in orienting the crystal so that the 
same facial area would always receive the same 
illumination. This was accomplished by the usual 
optical methods. The photographic plate was 
mounted normal to the incident beam. 

Before observations of any kind are made, it is 
of primary importance to have as “perfect” a 
crystal surface as possible. The theory calls for a 
well-defined beam. The divergence of the beam 
reflected from the crystal face must, of course, be 
narrower than the widths of any of the recorded 
maxima. The maxima are sharpest at the Bragg 
angle and become broad and diffuse as the crystal 
is rotated from its Bragg position. Hence the 
narrowest maximum is determined by the point 
of observation closest to the Bragg angle for the 
characteristic radiation producing the maximum. 
Observations cannot be made at the Bragg angle, 
for the diffuse temperature maximum and the 
Bragg reflection coincide. Thus the crystal is 
slightly rotated away from its Bragg angle so 
that the Bragg reflection is completely suppressed 
and will not interfere with the emitted tempera- 
ture maximum. It is found however, that most 
crystal surfaces are rough and warped, and this 
allows the Bragg reflection to persist over quite 
an angular range about the Bragg angle. 

Experimental difficulties have been encoun- 
tered at this point. In the original experiment 
with NaCl it was found that the Cu Ka line 
reflected from the (200) face could persist even 
though the crystal was rotated 60’ from its 
Bragg position. The diffuse maxima cannot be 
studied under such conditions for they are very 
nearly superimposed upon the persistent Cu Ka 
line, and it was only by properly etching the 
surface that the line radiation disappeared when 
the crystal was 30’ from the Bragg position. If 
the intensity of the line radiation is plotted as a 
function of the deviation of the crystal away 
from the Bragg position for the line, then one 
obtains a “surface distribution curve.”’ Etching 
the surface produces ‘‘surface distribution curves” 


of various shapes and widths. Large primary 
slits allow illumination of a larger facial area, and 
this increases the persistence of the Bragg 
reflection. 

Therefore, if observations are to be made on 
the sharp intense temperature maxima near the 
Bragg angle, a good crystal surface showing a 
rapid cut-off of the Bragg reflection must be used. 
Such a surface will exhibit a narrow distribution 
curve. After determining the point of observation 
closest to the Bragg angle, slit widths can be 
approximately adjusted to account for the 
narrowest diffuse maximum which will be emitted 
at that point. 

In the present experiment an extraordinary 
KCI crystal section was cleaved from its parent. 
The parent was an artificial crystal. The cleavage 
(200) surface was so perfect that the Bragg angles 
for each component of the Cu Ka doublet were 
quite sharply determined. Figure 2 shows the 
appearance of the doublet on rotation of the 
crystal through the Bragg angles for the doublet. 
In Fig. 2(a) Ka; is strongly reflected while Kag is 
barely visible. In Fig. 2(b) Kae has increased in 
intensity so that its intensity is greater than that 
of Ka. Figure 2(c) shows that Ka, has still 
decreased in intensity. The angles indicated in 
the figure are arbitrary. A surface distribution 
curve for this crystal has been made. Two curves 
are possible, one for each component of the 
doublet. The deviation of the crystal away from 
the average Bragg angle for the doublet is indi- 
cated by A. For A=6’ the intensity of the 


Fic. 2. Showing the change 
in intensity of each com- 
ponent of the Cu Ka doublet 
over a three-minute interval 
in the crystal angle. These 
lines are magnified. The lines 
tend to break up when the 
crystal is moved too far from 
the Bragg positions. 
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Devvuss Spot 
Fic. 3. Typical photograph for A=29’. The exposure 
time is 32 hours. The diffuseness of the temperature 
maximum can be seen extending into the Laue spot. A diffuse 
temperature halo can be seen surrounding the direct beam. 


emitted doublet became weak enough to make 
observations on the temperature maximum. A 
5-hour exposure on the photographic plate 
yielded a temperature maximum and only a very 
faint doublet which did not interfere. 

The necessity for such a narrow beam made the 
use of the ionization-chamber method impossible 
for recording the intensities of the diffuse maxima. 
For studying the widths of the lines, the analyzing 
slit in front of the ionization chamber must 
measure the differential intensity and must thus 
be narrower than the primary slits. Even then 
the intensity of a diffuse line for A=40’ say, may 
be about 10-4 that of the Cu Ka line emitted at 
the Bragg angle. 

In making observations, the crystal was 
oriented so that the cleavage face would reflect at 
the average Bragg angle for the Cu Ke doublet. 
All diffuse maxima were then due to the line 
radiation. The intensity of the continuous radi- 
ation is small in comparison with that of the 
Cu Ka lines, so that the diffuse scattering of the 
continuous radiation can be neglected. Since the 
crystal was quite perfect, observations began for 
A=6’ and extended up to A=41’. For A=6’ the 
exposure time to produce the maximum was 5 
hours, and for A=41’ the exposure time was 71 
hours. The exposure time varies approximately as 
1/A*. Thus any observations beyond A=41’ 
would require extremely long exposures. 


EFFECT OF TEMPERATURE 


All observations in this experiment were made 
at room temperature. However, the diffuse 
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maxima have been studied over a wide range of 
temperatures. Lonsdale and workers have ob- 
served them at liquid-air temperatures. They 
find the diffuse spot pattern disappears almost 
completely at that temperature. Raman and 
Nilakantan have observed them from room 
temperature up to 275°C and find an increase in 
intensity at higher temperatures. Preston has 
made observations at 500°C and finds increased 
intensity. The observations include a variety of 
crystals. 


EXPERIMENTAL RESULTS 


Figure 3 shows a typical photograph. For this 
case 4=29’ and the exposure time is 32 hours. 
The diffuseness of the temperature maximum 
may not be apparent from the reproduction 
given, but a photometer tracing reveals that it 
extends well into the Laue spot next to it. The 
geometrical patterns are the Laue reflections. 
The Laue reflection next to the diffuse line is very 
uniform and regular, indicating the quality of the 
crystal face. 

According to theory, the direction of scattering 
corresponding to an intensity maximum of the 
diffuse radiation lies in the plane of incidence, 
and the scattering angle is 26,,=26g+2A sin? 
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Fic. 4. Shape of the diffraction maxima for different 
values of |A]. 
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If the incident beam is a narrow pencil, the 
intensity near a maximum is 


K 
(A sin 208)*+x°+¢" 


where x and ¢ are the angular deviations from 
the maximum in the plane of incidence and 
normal to it. If A, x, and g remain small, K can be 
treated as a constant. However, in this experi- 
ment a slit was used. The intensity measurements 
were made in the plane of incidence, so that the 
above equation must be integrated with respect 
to g. The result is 


J2« tan (¢n/y)/y, y=[(A sin 268)?+ x7}, 


where ¢n=//2R and / is the height of the slit. 
R is the distance from the crystal to the point of 
observation. In the present experiment ¢, 
was 85.2’. 

Figure 4 shows the shape of the diffuse temper- 
ature maxima for various values of A sin 26%, 6 
being the average Bragg angle for Cu Ka;,». The 
experimental curves were calculated from photo- 
meter tracings of the plates. The agreement 
between observed and theoretical shapes of the 
maxima is quite good. 

In Fig. 5 x; represents the half-widths at half- 
maximum. x; is plotted as a function of A sin 26. 
The circles represent single measurements. In 
some cases it was necessary to average the values 
obtained from more than one tracing. These 
latter points are indicated by crosses. According 
to theory the maxima do not appear at an angle 
26, but at an angle 26g+2A sin? 63. This means 
the lines are shifted by a matter of some seconds. 
This has been observed on some of the plates, 
but it is too small to measure accurately. It will 
change the experimental half-widths by a very 
small amount only. 


20'} 
Fic. 5. Variation } / 

of half-width with '5- - 

THE ORY 
EXPT. 
10 18 20 25 30 

4 SIN2@,——> 
SUMMARY 


All measurements have been made on maxima 
due to the line radiation present in the incident 
beam, and only a simple crystal, KCI, has been 
studied. The crystal has been found to exhibit a 
perfect surface for this study. The observations 
extend from A= +6’ to A= +41’. The shapes and 
widths of the maxima have been found to agree 
quite well with the theory. 

A study of the variation of the peak intensities 
of the maxima with A sin 26, has previously been 
made.? This variation of intensity is the back- 
ground present about each reflection in an 
oscillation photograph. For every Laue reflection 
there is associated with it a diffuse temperature 
maximum. However, maxima may arise in direc- 
tions other than those of the Laue directions in 
the manner indicated. The intensity of the 
maxima has been shown to increase with in- 
creased rise in temperature. The theory indicates 
that there should be a slight increase in intensity 
of the diffuse radiation in a region about the zero 
angle of scattering. Figure 3 shows that it is 
manifest as a halo about the transmitted direct 
beam. 

The author wishes to thank Dr. Zachariasen 
for suggesting the problem and for his valuable 
aid and criticism during the investigation. 
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The Measurement of Self-Diffusion in Zinc* 


FLoyp R. Banxks** 
Randal Morgan Laboratory of Physics, University of Pennsylvania, Philadelphia, Pennsylvania 
(Received October 3, 1940) 


Measurements of self-diffusion along the C axis of zinc single crystals have been carried out. 
Long-lived radioactive Zn® was used as an indicator. Measurements were made at six different 
temperatures, enabling a determination of the activation energy, which is 19,600 calories per 
mol. The diffusion coefficient at 410.4°C is 9.89 10-4 cm*/day. 


INTRODUCTION 


IF FUSION in metals is important both to 
the metallurgist and to the physicist.' It is 
important to the former because phenomena 
such as carburization of steel, production of 
corrosion- and heat-resisting substances, rates of 
hardening and transformation, and heat treat- 
ment of castings, all depend on diffusion rates. 
Self-diffusion is especially important because of 
the information it may reveal concerning creep 
in solids. It is important to the physicist because 
it will elucidate the mechanism of rate processes 
in solids, since it involves atomic movements. 
Inquiries into the subject of self-diffusion (in 
which the atoms of a lattice are replaced by 
atoms of their own kind) led von Hevesy? and 
his collaborators to the use of radioactive isotopes 
as indicators. The property of radioactivity of an 
atom A permits its labeling or “tagging’’ and 
hence affords a means of studying its progress 
through the lattice of nonradioactive solvent A. 
By such means Hevesy, Seith, and Keil? studied 
self-diffusion in lead and Seith and Keil*4 
studied self-diffusion in lead and bismuth, with 
thorium B (which is an isotope of lead) and 
thorium C (which is an isotope of bismuth). 
Prior to the advent of artificial radioactivity 
the use of radioactive isotopes as indicators in 


* This investigation has been carried out in cooperation 
with Dr. H. M. Day, Research Laboratory of the General 
Electric Company, Schenectady, New York. The prepa- 
ration of accurately oriented, strain-free single crystals, 
and the precision heat treatments will be described in 
publications dealing with other phases of the diffusion 
problem. 

** Now at Morgan State College, Baltimore, Maryland. 

'™R. F. Mehl, Trans. A. I. M. E. Inst. of Metals Div. 122, 
11 (1936). 

* Hevesy, Seith and Keil, Zeits. f. Physik 74, 197 (1932). 

* W. Seith, Zeits. f. Elektrochemie 39, 538 (1933). 
(1933) Seith and A. Keil, Zeits. f. Metallkunde 25, 104 


studies of self-diffusion was limited to the 
heavier elements since these were the only ones 
which had active isotopes in nature. With 
artificially produced radioactive isotopes, how- 
ever, means of studying the lighter elements, 
such as zinc and copper (elements of special 
interest to the metallurgist) are available. Self- 
diffusion in gold® and copper’ has been studied 
by using artificial radioactive isotopes. The 
present paper deals with self-diffusion in zinc. 


DESCRIPTION AND DISCUSSION OF 
PRESENT INVESTIGATION 


Radioactive Zn® is used as an indicator in 
this study. It is electrolytically deposited on a 
flat, polished, and etched surface of a single 
crystal of high purity zinc; the sample is then 
held for a fixed time at an elevated temperature 
to allow diffusion to proceed in the crystal, and 
finally cut into sections on a lathe at room tem- 
perature. The activity of each section is measured 
with a Geiger counter or an electroscope. From 
the weight of each section, the thickness can be 
computed, and from the thickness and activity 
the diffusion coefficient is obtained. 

Zn® (half-life eight months) can be obtained 
free from inactive Zn by proton bombardment 
of copper: Cu®(p, 2)Zn®. The Zn® produced by 
this proton-neutron reaction emits positrons and 
conversion electrons of a strong gamma-ray. 
Zn® can also be produced from Zn* by deuteron 
bombardment: Zn"™(d, p)Zn®. 

It is not possible to determine the diffusion 
coefficient by a method such as that of Hevesy, 


5 A. Sagrubskij, Physik. Zeits. Sowjetunion 12, 118 (1937). 
®H. A. C. McKay, Trans. Faraday Soc. 34, 345 (1938). 

7B. V. Rollin, Phys. Rev. 55, 231 (1939). 

8 Steigman, Shockley and Nix, Phys. Rev. 56, 13 (1939). 
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Fic. 1. Counter measurements on crystals of group A. 
Logarithm of concentration for each section plotted against 
square of mean depth in crystal. Temperature of heat 
treatment—412.5°C; time of heat treatment—12.5 hours 
for crystals I and II, 24 hours for crystal III, 6 hours for 
crystal IV; number of sections—ten, ten, ten, and six, 
respectively. 


Seith and Keil,?~* ® viz., evaporation of Zn® onto 
the surface of a crystal of ordinary zinc and 
measurement of the radioactive decay before and 
after diffusion, because there is no fixed range 
for electrons. It is possible, however, to solve 
the diffusion equation and take into account the 
approximately exponential character of the ab- 
sorption of a continuous energy spectrum of 


electrons (the emitted radiation) in the sub- 


stance employed, as Steigman, Shockley and 
Nix® did. This would lead to a very uncertain 
interpretation of the data in the case of Zn®, 
for the ionizing radiation from this radioactive 
element is complicated by the emission of strong 
gamma-radiation. Methods may be devised, how- 
ever, which do not depend on calculating the 
amount of diffusion which occurs from the de- 
crease in intensity of the radiation measured at 
the surface. The procedure outlined at the 
beginning of this section is such a method. 

Zn® is separated from the copper from which 
it is produced by dissolving the surface of the 
copper containing the Zn® in 1—1 nitric acid. 
Chlorides are made of the Zn® and copper, and 
the copper precipitated by hydrogen sulphide 
from a 0.3 normal HCI solution. Actually, there 
will be more copper ions present in the filtrate 
than Zn® ions, because of the degree of efficiency 
of the precipitation process and because there is 


°G. v. Hevesy and W. Seith, Zeits. f. Physik 56, 790 
(1929). 


such a small quantity of Zn® initially. Several 
precipitations are made until further bubbling of 
hydrogen sulphide through the solution brings 
down no more copper. The precipitate is tested 
in each case for its activity to ascertain whether 
any of the Zn® is brought down with the copper. 
After this procedure practically all of the Zn® 
remains in the solution, which is then concen- 
trated down to dryness, and the residue dissolved 
in 5 ce of water. This solution serves as the eclec- 
trolyte in the electrolysis process. 

Single crystals were grown from high purity 
Zn and were used in the form of cylinders about 
0.5 inch in diameter and 0.3 inch long. All the 
crystals had the same crystallographic orienta- 
tion (diffusion interface perpendicular to the C 
axis), and their flat, etched, and polished sur- 
faces were perfect reflectors. Each crystal served 
as the cathode during the electrolysis, a portion 
(0.5 cc in 3 ce of water) of the Zn®™ solution men- 
tioned above was the electrolyte, and a spiral 
platinum wire was the anode. The time of the 
electrolysis varied from forty-five minutes to six 
hours, the average voltage used was 4.5 volts, 
the average current 10 milliamperes. None of 
these factors seemed to be critical. During the 
electrolysis the crystals generally became tar- 
nished. 

A rotating cathode was used to provide 
stirring of the electrolyte to have a uniform 
distribution of the Zn® ions throughout the 
solution, to prevent bubbles from gathering on 
the surface of the cathode, and to permit 
selective deposition of Zn® in the presence of 
copper. A rotating cathode!" is indispensable 
for selective deposition of Zn®, for it is well 
known that copper will plate out before zinc. 
A high speed rotating cathode obviates this. 
Also, the cathode was dipped in and out of the 
solution during the electrolysis to prevent the 
formation of bubbles, for they would destroy 
the intimacy of the contact between the surface 
of the zinc crystal and the Zn®. 

After electrolysis the impregnated crystals 
were held for a fixed time at elevated tempera- 
tures to allow diffusion to proceed. Variations 
in temperature during the course of the heat 

Ww” C, G. Fink and F. A. Rohrman, Trans. Am. Electro- 


chem. Soc. 58, 403 (1930). 
" J. Steigman, Phys. Rev. 53, 771 (1938). 
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Fic, 2. Electroscope measurements on crystal 1V, group 
A. Logarithm of concentration for each section plotted 
against square of mean depth in crystal. Temperature of 
heat treatment—412.5°C ; time of heat treatment—6 hours; 
number of sections —six. 


treatment were always less than 0.1°C. At room 
temperature the crystals were cut into sections 
approximately 0.002 inch thick on a lathe, the 
cut being parallel to the original face. The 
sections were in the form of ribbons, as cut by 
the lathe tool. It was not necessary to have the 
section of predetermined or accurately equal 
thickness; for each section was weighed, and the 
thickness computed from the weight, the known 
density of zinc (taken here to be 7.10 g/cc), and 
the area of the face from which the sections were 
turned. General freedom from edge effects was 
ensured by turning off the cylindrical surface of 
the crystal to a depth of approximately 0.005 
inch before the sections were taken off. 

The activity of each section was measured 
either with a Lauritsen” electroscope or with a 
Geiger-Miiller counter. For the electroscope 
measurements each section was placed in the 
clectroscope can. Most of the measurements were 
made with a thin glass-walled Geiger counter, 
directly under which each section was placed 
(in a small pill box) in approximately the same 
position. The background of this counter was 
approximately 17 counts per minute. After diffu- 
sion the activities of the first sections were at 
least 100 counts per minute above the back- 
ground, the activities of the last sections at 
least 5 counts per minute above the background. 
About 2000 counts per section were taken, so 
that the precision ranged from one percent for 
the first sections to nine percent for the last. 

It was not necessary to compare each measure- 
ment with a standard to allow for radioactive 
decay, as McKay* and Steigman, Shockley, and 


" C.C. Lauritsen and T. Lauritsen, Rev. Sci. Inst. 8, 438 
(1937). 
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Nix® did. Because of the long half-life of Zn® 
(eight months) the decay during any series of 
measurements is negligible. Also, it was not 
necessary to correct for the absorption of elec- 
trons in matter since each section was only 
approximately 0.002 inch thick (except for 
crystal IV, group A, where each section was 
0.004 inch thick), and geometrically similar to 
every other section. 

The electroscope and Geiger-counter measure- 
ments provided independent methods of deter- 
mining the respective activities (and hence 
diffusion coefficients D) and the close agreement 
between these two methods (e.g., see results of 
crystal IV, group A) was construed to indicate 
reliability of the present technique in the study 
of the phenomenon in question. 

The diffusion equation is 


dc/dt = Dd*c/dx* (1) 
if there is no evaporation to or from the surface. 
Here c is the concentration of Zn®, ¢ the time for 
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Fic. 3. Counter measurements on crystals of groups B 
and C. Logarithm of concentration for each section plotted 
against square of mean depth in crystal. Temperature of 
heat treatment—400°C and 374.4°C, respectively; time of 
heat treatment—6 hours and 7.5 hours, respectively; 
number of sections—ten. 
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the diffusion, « the axial coordinate in the sample, 
and D the diffusion coefficient. D is assumed to 
be independent of c. The solution of this equation 
is given in this form: 


c=(rDt) exp [—x*/4Dt), (2) 


in which the initial thickness of the active layer 
is assumed infinitesimal. If In ¢ is plotted against 
#, which represents the square of the average 
depth of the section considered below the surface 
of the crystal, a straight line results of intercept 
—In (rDt)', and of slope —}Dt. Thus, if ¢ and x 
are known, D can be obtained directly from the 
graph, for a given ¢. The concentration in each 
section is obtained in relative units from the 
activity by dividing by the thickness of the 
section (c,=A,/t,). This introduces an error of 
less than 1 percent (for the small distances 
dealt with in the experiment), which is well 
under the experimental error. 
From Eq. (2) 
(3) 
4 In 


From the half-value of the concentration, i.e., 
the point where the concentration has fallen to 
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Fic. 4. Counter 
measurements on 
crystals of groups 
Dand F. Logarithm 
of concentration for 
each section plotted 
against square of 
mean depth in crys- 
tal. Temperature of 
heat treatment— 
364.7°C and 355°C, 
respectively; time 
of heat treatment— 
8 hours and 9 hours, 
respectively; num- 
ber of sections— 
ten. 
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Fic. 5. Counter measurements on the only crystal of 
group /. Logarithm of concentration for each section 
plotted against square of mean depth in crystal. Tempera- 
ture of heat treatment —410.4°C; time of heat treatment— 
3.5 hours; number of sections—ten. 


one-half of its original value we have 
Dt=0.3618;, (4) 


where Z,? is the square of the abscissa where the 
half-value ordinate occurs. 

is determined from the graph (of In ¢ vs. 
and ¢t is known. Therefore D is obtained directly. 


RESULTS 


The rate of diffusion was determined at six 
different temperatures for six groups of single 
crystals: Group A consisted of four crystals, 
group £, one, and all other groups of two 
crystals each. 

Figure 1 displays the graphs of the counter 
measurements for the crystals I-IV of group A. 
Figure 2 is the graph for crystal IV, group A, 
for electroscope measurements. The agreement 
between the counter and electroscope measure- 
ments on the sections of this crystal is quite 
satisfactory. Figures 3, 4, and 5 are the graphs 
for crystals of groups B, C, D, E, and F. 

The temperature dependence of D is given by 


D=Ae~-URT, (5) 
where A is the activation constant, Q the activa- 


tion energy, R the gas constant, 7 the absolute 
temperature, and e is the logarithmic base. 


In D=In A—(Q/RT). (6) 
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DIFFUSION COEFFICIENT (CmM7/DaY) 


Fic. 6. Diffusion coefficient as a 
function of reciprocal temperature. 
The crystals of group A were heat 
treated before the precision furnance 
control had been developed. 
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If In D is plotted against 1/7 a straight line of 
intercept In A, slope —Q/R results. Q, the 
activation energy is obtained directly from the 
graph. Figure 6 shows the temperature de- 
pendence of D (a plot of D against 1/T). The 
value obtained for Q is 19,600 cal./mol. 

Q can be calculated also from the Dushman- 
Langmuir" equation 


where 6 is the interatomic distance, N is Avo- 
gadro’s number, / is Planck’s constant, D the 
diffusion coefficient, Q the activation energy, R 
the gas constant, 7 the absolute temperature, 
and e the logarithmic base. The average calcu- 
lated value of Q is 22,200 cal./mol. 


DISCUSSION OF RESULTS 


Different boundary conditions prevailing at 
the interface between the diffusing solute and 
the ordinary zinc crystal will alter the results 
because of the possibility of the trapping of 
active zinc on the surface. Some of the plated 
zinc may be in good contact and start to diffuse 
immediately, some may be in poor contact and 
start to diffuse later, and some may not start 
at all. As a matter of fact, it is seen from the 
graphs that the first point of the curve for each 
crystal seems to be somewhat high. It is assumed 


%S. Dushman and I. Langmuir, Phys. Rev. 20, 113 
(1922). 


157 158 160 


that the greater part of the diffusion takes place 
normally. Justification for this assumption is in 
the agreement between the experimental and 
theoretical curves. 

Possible sources of error in the determination 
of D are: (1) variations in the temperature of 
heat treatment (practically negligible since the 
furnace control was good to 0.1°C; (2) errors in 
weighing (less than 1 percent); (3) slightly 
different geometric orientations of the sections 
with respect to the counter in measuring activi- 
ties; (4) assumption of independence of D and 
concentration c; (5) assumption that the thick- 
ness of the initial layer is not comparable with 
the depth of penetration (and hence the evapora- 
tion is not important) ; (6) statistical fluctuations 
in the counts made on each section ; (7) location 
of the average line through the points; and (8) 
variation in time at the temperature of the heat 
treatment (a maximum of 3 percent). 


CONCLUSION 


The method described presents a very satis- 
factory way of determining the coefficient of 
self-diffusion D, and the activation energy Q, in 
zinc. It is especially useful in that it provides a 
simple and direct procedure for getting an 
average value for D. 
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Louis N. Ridenour, who initiated the experi- 
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The radiation emitted when conduction electrons fill vacancies in the L shell is spread over 
a band whose shape represents the energy distribution of the conduction electrons. This band 
has been studied spectrophotometrically for pure sodium, magnesium and aluminum. The 
observed shapes are in general agreement with theoretical predictions and confirm the results 
of previous experimental studies. In each case a “tail” has been observed at the low frequency 
end of the band, for which alternative explanations are offered. With the x-ray continuum as a 
standard of intensity, the total intensity of the Z band increases nineteen-fold from sodium to 


aluminum. 


HE electronic structure of a metal can be 
investigated by a study of the x-radiation 
emitted by the metal when one of its conduction 
electrons falls into a vacancy in an inner shell. 
Since the conduction electrons have energies 
throughout a range of several electron volts, this 
x-radiation is a band of corresponding width. 
The distribution in energy of the intensity in this 
band is related to the distribution in energy of 
the conduction electrons; thus the definite upper 
limit of the energy of the conduction electrons 
is evidenced by the fact that the x-ray emis- 
sion band is sharply bounded at its high fre- 
quency end. 

X-ray bands of this origin have been observed 
repeatedly. The recent work of Bearden and 
collaborators! provides examples in the region of 
1.5A. In principle, it is preferable to observe the 
structure of the bands emitted in the ultra-soft 
x-ray region (A>10A), for here the frequency of 
the radiation is not large in comparison with the 
frequency spread of the emission band and high 
resolution technique is not necessary. Of several 


: ai Friedman and W. W. Beeman, Phys. Rev. 58, 400 


reports of bands emitted in this region, those of 
O'Bryan and Skinner® and of Skinner® are out- 
standing in that they present and discuss reliable 
photometric descriptions of the radiation from 
very clean surfaces of certain metals. In the 
emission of ultra-soft x-rays, the surface of the 
target must be extremely clean, since all photons 
except those from the surface are absorbed in the 
target. Foreign elements on the target not only 
introduce their own characteristic spectra, but 
also affect the characteristic radiation of the 
target element itself. The numerous reports of 
ultra-soft x-rays by Siegbahn and Magnusson‘ 
and others at Uppsala must be regarded as 
qualitative, since their targets were admittedly 
contaminated with carbon, nitrogen and oxygen. 

The present work was undertaken as an at- 
tempt to confirm and extend previous work in 
the photometry of the band of radiation emitted 
as conduction electrons fall to Zi; or Li levels 
made vacant by electron bombardment. 

2H. M. O'Bryan and H. W. B. Skinner, Phys. Rev. 45, 
370 (1934). 

3H. W. B. Skinner, Prog. Rep. Phys. Soc. 5, 257 (1939). 


* M. Siegbahn and T. Magnusson, Zeits. f. Physik 62, 435 
(1930); 87, 291 (1934); 88, 559 (1934); 95, 133 (1935). 
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EXPERIMENTAL 


The vacuum spectrograph comprised a glass 
grating of radius of curvature about 1.5 meters 
and effective width 2.5 cm, ruled with 30,000 
lines per inch and adjusted for about 5.7° grazing 
angle of incidence; a vertical slit, made of two 
fragments of a razor blade; and a plateholder, 
embracing the wave-length range from zero to 
about 740A, which could be raised or lowered 
without breaking the vacuum. (This motion of 
the plateholder caused no wave-length displace- 
ment of spectrum lines, within the wave-length 
accuracy demanded in the present work.) A 
horizontal slot before the plate made it possible 
successively to bring adjacent zones of the plate 
into the same part of the x-ray beam, and an 
externally controlled shutter was used to begin 
and end the exposures. The dispersion varied 
from about 1.0A/mm at 100A to 2.2A/mm at 
740A. 

The x-ray tube was similar to that of O’Bryan 
and Skinner.’ It consisted of a 2-liter Pyrex bulb 
having arms with tapers into which were inserted 
the filament and target assemblies. The x-ray 
filament was a tungsten helix, 2 cm long, wound 
on a 60-mil mandrel out of 10-mil wire 40 cm in 
length. The filament was partly surrounded by a 
tantalum shield to prevent direct illumination of 
the slit by the filament. An evaporator for 
aluminum, inserted in a third sidearm, consisted 
of a 20-mil tungsten filament in the form of a 
helix, the coils of which were charged with small 
loops of aluminum wire; magnesium was evapo- 
ated from tungsten baskets, while in the case of 
sodium a special distilling furnace was employed. 
The metals used in coating the target were of high 
purity. In no case was the impurity content 
greater than 0.1 percent. 

The target was a thin-walled copper tube of 
square section. It was free to rotate so that any 
one of the four faces could be placed in position 
for the exposure, the target being so turned that 
the radiation entering the slit made an angle of 
45° with the normal to the face in use. The target 
was placed at 6 mm from the filament and 14 cm 
from the spectrograph slit, which was protected 
by a “‘fore-slit.’’ None of the grease (Apiezon L) 
used to lubricate the tapers at the three arms 
was nearer than 15 cm from the target; therefore 


the grease was not appreciably warmed by the 
heat of the filament. During all aluminum and 
magnesium runs the target was water cooled. 

In the case of sodium it was necessary to keep 
the target temperature below 0°C in order to 
prevent the sodium from melting. This was ac- 
complished by circulating a liquid (alcohol), 
cooled in passing through a copper coil sur- 
rounded by solid carbon dioxide and alcohol. The 
circulation was maintained by a gear pump. The 
temperature of the alcohol leaving the target 
assembly was about —45°C. Provision was made 
for keeping the greased taper of the target near 
room temperature in order to prevent “‘freezing”’ 
of the metal to the glass. 


PROCEDURE 


The x-ray exposures were so planned that the 
relative intensity could be found as a function of 
the distance along the length of the spectrogram. 
The target was bombarded for a couple of hours 
to outgas it, the shutter being closed; two or 
three times during this conditioning the target 
was coated by heating the evaporator. When the 
ionization gauge read about 4X10~-* mm Hg, the 
target was recoated and the shutter opened to 
begin the exposure. During the exposure, the 
tube was connected directly to a transformer 
giving 1450 volts, r.m.s.; by adjustment of the 
filament current (about 5.4 amp. a.c.) the d.c. 
milliammeter in series with the target was held 
at 200 ma. 

In the case of aluminum, the exposure time 
was 50 minutes with the plateholder up, and 100 
minutes with the plateholder down; for an ad- 
ditional 50 minutes the plateholder was up and 
the target potential was off, to equalize the 
slight light fog from the x-ray filament upon the 
two zones of the plate. There was no visible 
fouling of the target, but to insure ideal surface 
conditions the target was recoated at 5-minute 
intervals. (It was found that aging the aluminum 
target for 10 minutes gave a general reduction of 
the intensity by a few percent, with no obvious 
change in the wave-length distribution of rela- 
tive intensity.) Since during the exposures the 
pressure fell typically from 4X10~ to 1X10 
mm Hg, it was considered necessary to break up 
the total of 200 minutes of exposure into ten 
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Fic. 1. Densitometer curves of a one-hour and a two-hour exposure of the aluminum Z band on the same plate, 
illustrating the sensitometric method. The wave-length range is 162 to 218A. 


cycles to prevent favoring either the “up”’ or 
the “down” portions. During each cycle the 
plateholder was up (with x-rays turned on) for 
5 minutes, up (with x-rays off) for 5 minutes, 
and down (with x-rays on) for 10 minutes. This 
routine, here described for aluminum, provides 
two juxtaposed exposures which have had the 
same treatment as regards light fog, and which 
differ only in respect to their x-ray illumination. 
Two good aluminum plates were made in 
this way. 

The seven good magnesium plates were ex- 
posed by a routine similar to the above, but 
there were six cycles of 4 times 10 minutes each, 
giving x-ray exposure times of 60 and 120 
minutes; the target was recoated at 10-minute 
intervals. For sodium, all four faces were simul- 
taneously coated at hourly intervals, so that 
each face was used for 15 minutes before recoat- 
ing; there were eight one-hour cycles, and the 
four target faces were in each cycle used in 
sequence, permuted variously to avoid beginning 
each cycle with the same face. Thus on each 
sodium plate the x-ray exposure times were 120 
and 240 minutes; two plates were so made, the 
tube current being limited to 100 ma to prevent 
the coating of sodium from melting. 

The respective slit widths for sodium, mag- 
nesium and aluminum were roughly 0.32 mm, 
0.12 mm, 0.08 mm. The wave-length at the peak 
of the sodium band was determined by compari- 


son with the line spectrum of air from a disrup- 
tive capillary discharge to which the plate was 
exposed after the x-ray tube had been removed. 
Ilford QI plates were used in all photometric 
runs. These seem nearly as sensitive as Schu- 
mann plates, and their sensitivity is fairly free 
of local irregularities. 

The photometric process by which the intensity 
distribution in the band was determined from the 
spectrogram is reminiscent of a method used by 
Skinner and Johnston® for comparing spectrum 
line intensities. The process depends on the as- 
sumption that the plate sensitivity is independ- 
ent of wave-length throughout the band and also 
that the reciprocity law is valid. No assumption 
need be made as to the characteristic curve of the 
plate; in fact, the characteristic curve is deduced 
in the course of this process. 

Densitometer curves from the two exposures 
on the same plate are recorded together as in 
Fig. 1. Since the exposures stand in the ratio 2:1, 
it may be seen that the intensity at the point x2 
on the faint trace is half the intensity at x,, that 
the intensity at x; is half the intensity at x2, etc. 
On the basis of these observations a preliminary 
characteristic curve may be sketched, showing 
the relation between the intensity and the trans- 
parency of the fainter exposure. (For convenience 
the deflection of the densitometer galvanometer 


5H. W. B. Skinner and J. E. Johnston, Proc. Roy. Soc. 
A161, 420 (1937). 
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may be used in place of transparency.) A second 
preliminary characteristic curve may similarly 
be obtained from the points x6, x7, etc., x» being 
arbitrarily chosen. The final characteristic curve 
may then be obtained by superposing two or 
more such preliminary curves. The galvanometer 
deflection is now tabulated as a function of 
distance along the plate; with the aid of the 
characteristic curve these data now yield values 
of (x), the relative intensity per unit distance 
along the plate.* The values of I(x) when multi- 
plied by dx/dV give I(V); here V is the energy 
in electron volts of the photons incident at any 
position x on the plate and J(V) is the relative 
intensity of the radiation in the band per 
electron-volt interval. The quantity dx/dV is 
related to the dispersion of the spectrograph. It 
must be pointed out that this scheme of pho- 
tometry gives results which represent the x- 
radiation only and are not influenced by uniform 


light fogging. 


TABLE I. Intensity distribution in the L bands. 


SopIum MAGNESIUM ALUMINUM 
V v3I(V) V 

25.00 1.6 39.5 1.1 58.0 4 
25.50 2.4 40.0 1.3 59.0 6 
25.75 5.6 40.5 2.0 60.0 9 
26.00 7.5 41.0 4.2 61.0 13 
26.25 9.4 41.5 6.6 62.0 21 
26.50 12 42.0 10 63.0 33 
26.75 14 42.5 17 64.0 48 
27.00 18 43.0 27 65.0 60 
27.25 21 43.5 39 66.0 67 
27.50 28 44.0 52 67.0 71 
27.75 41 44.5 65 68.0 71 
28.00 53 45.0 75 69.0 69 
28.25 64 45.5 81 70.0 67 
28.50 74 46.0 86 70.5 68 
28.75 81 46.5 89 71.0 71 
29.00 87 47.0 91 71.5 82 
29.25 94 47.5 91 72.0 100 
29.50 97 48.0 89 72.2 97 
29.75 100 48.4 85 72.4 73 
30.10 97 48.8 87 72.6 29 
30.20 93 48.9 92 72.8 12 
30.25 79 49.0 97 73.0 2 
30.30 61 49.1 100 73.5 1 
30.35 40 49.2 83 74.0 0 
30.40 25 49.3 43 

30.45 17 49.4 15 

30.50 11 49.5 7.3 

30.55 6.1 49.6 3.9 

30.70 3.7 50.0 1.1 

31.25 1.9 50.4 0.4 

31.75 0.8 51.0 0.4 


6 This “‘staircase’’ method, originated by Professor L. G. 
Parratt for a somewhat different purpose, is described by 
him in a paper now in preparation. 
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RESULTS 


For each element there has been observed a 
characteristic band superimposed on a continu- 
ous background. (This background is assumed to 
be the true x-ray continuum, although it may 
be due in part to x-rays scattered from the 
grating.) The results on these spectra, averaged 
for each element, are presented in Tables I and II 
and Fig. 2. Table I gives, for several values of V, 
the intensity distribution as »~*J(V) where is 
the frequency of the Z radiation; the continuous 
spectrum, evident on both sides of the Z band, 
has been subtracted off. From the work of Jones, 
Mott and Skinner’ and Skinner,* it is expected 
that the curve v“*J(V) vs. V for an Ly wy band 
will be parabolic near the low frequency end of 
the band. These curves are shown for sodium, 
magnesium and aluminum in Fig. 2, together 
with the parabolas that best fit them. The curves 
resemble those of Skinner.* 

Certain numerical characteristics are compiled 
in Table II. For sodium, (peak) is the wave- 
length for the greatest value of J(A), which is 
the relative intensity per unit wave-length 
interval; while the quantity designated as 
V(peak) is the photon energy when J(V) is 
maximum. For the remaining elements V(peak) 
is calculated directly from the published values 
of A(grenze) which are here listed as \(peak). 
“Maximum V, extrapolated” is the intercept of 
the steepest tangent to the high V end of the 
vJ(V) curve. “Band width” is defined as the 
separation between this intercept and the foot of 
the ‘‘best parabola.”” For comparison, there are 
also listed the data of O’Bryan and Skinner? as 


TABLE II. Summary of results. 


MAGNE- ALUMI- 


SoDIUM SIUM NUM 

\ (peak) (A) 407.6+0.4 250.20* 170.49** 
Vv (ev) 29.9 49.1 72.0 
Maximum V, extrapolated (ev) 30.4 49.4 72.8 
V at foot of best parabola (ev) 27.35 43.0 62.2 
d at tip of tail (A) 487 307 212+ 
V at tip of tail (ev) 25.2 40 58— 
Band width (observed) (ev) 3.05+0.10 6.4 10.6 
Band width (Skinner) “) (ev) 3.0+0.2 7.640.3 13.2+0.5 
Sommerfeld band width®) (ev) 3.2 7.3 11.9 
Area ratio 

(tail to total band) 0.104 0.065 0.073 
Intensity ratio 

(peak to continuum) 1.7 7 11 
Intensity ratio 

(total band to continuum) 4.2 34 78 


* M. Siegbahn and T. Magnusson, Zeits. f. Physik 87, 291 (1934). 
** M. Siegbahn, Ergeb. d. exakt. Naturwiss. 16, 123 (1937). 


7 Jones, Mott and Skinner, Phys. Rev. 45, 379 (1934). 
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recalculated by Skinner,® and band width pre- 
dictions on the basis of the Sommerfeld free 
electron model. 

The tail at the low energy end of the bands, 
reported also by O’Bryan and Skinner, is interest- 
ing. Table II gives a comparison of the area, 


_(measured on the »v“*J(V) graph, between the 


curve and the best parabola) of the tail with the 
total area under the »~*J(V) curve. This ratio of 
areas is roughly independent of atomic number. 
The origin of the tails is uncertain. Skinner* sug- 
gests that they are caused by a softening of the 
radiation in the band, as by the Raman effect. 
Seitz® attributes them to ‘“‘excitation levels,” 
i.e., to extra levels below the valence band, that 
occur wherever a completed shell of an atom has 
been ionized; since this is, of course, the case at 
each point in the lattice where a photon is pro- 
duced, the excitation states might affect the 
emission band. 

It may be pointed out that thermal fluctua- 
tions in the lattice constant might, according to 
the Sommerfeld theory, cause a diffuseness at 
the foot of the parabola, observed as a tail. 
Alternatively one may propose a double jump 
hypothesis to account for the tails: that as one 
conduction electron falls to the ionized L shell a 
second conduction electron jumps into an un- 
occupied conduction level; in this ‘‘semi-Auger”’ 
process a photon is emitted, of which the energy 
is less than the energy loss of the first electron.* 

The bands have very unequal intensities for 
the different elements. A convenient measure for 
the relative intensities of the bands is afforded 
by the continuous spectrum, since the con- 
tinuum is not expected to vary greatly among 
the present three elements. The maximum /(V) 
is seen in Table II to increase rapidly in relation 
to the continuum, with increasing atomic num- 
ber. It is, perhaps, of theoretical interest to 
compare, for each element, the total band in- 
tensity {J(V)dV with the intensity per volt of 
the continuous spectrum. This ratio is also seen 

y A Seitz, Modern Theory of Solids (McGraw Hill, 1940), 

Note added in proof —A semi-Auger process, inde- 
pendently considered by Skinner and the authors, (H. W.B. 
Skinner, Phil. Trans. Roy. Soc. A239, 95 (1940)) (above) 
as resulting in the tails that appear on the low energy side 
of the emission bands, is rejected by Skinner in favor of 
an alternative hypothesis on the ground that this process 
would produce tails at both sides of the bands. This rejec- 


tion is however invalid: the semi-Auger process can de- 
crease but never increase the energy of an emitted photon. 


4 

Na 


-3-2-1 O!1234 5 electron volts 


Al 


-4-3-2-10 1234567891011 


Fic. 2. Intensity distribution in the L bands corrected 
for continuum, divided by v*. The dotted curves are the 
parabolas which best fit the data. The zero of the energy 
scale is at the foot of the parabolas. 


to increase greatly with atomic number. This 
increase is to be expected; in fact a semi- 
quantitative calculation involving the radiative 
transition probability, Auger effect and ioniza- 
tion probability affords a prediction that the 
foregoing ratio shall, between sodium and 
aluminum, show a 12-fold increase. 

It is a pleasure to acknowledge helpful con- 
versations with Professor Debye and Professor 
Parratt, and the assistance generously given by 
Professor Bethe in the interpretation of the 
bands. The design and operation of the x-ray 
tube were influenced by the good advice of Pro- 
fessor H. M. O'Bryan of Georgetown University. 
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Term Values in V II and Mn II 


D. S. BowMAN 
Allegheny College, Meadville, Pennsylvania 


(Received December 23, 1940) 


The term values for d’s and d*s have been calculated in Russell-Saunders coupling. For some 
of the terms of the same kind in d's, only the sum of the energies is obtained. The first approxi- 
mation to the separations of the levels is also found for d*s. The theoretical formulas including 
Ostrofsky’s results for d‘ are then applied to the observed values of VII and Mn II. 


HE term values for d*s and d*s have been 
calculated in Russell-Saunders coupling. 
Only the center of gravity of the levels with 
different values of J belonging to a single term 
have been found in ds. In d*s the separations of 
the levels have also been calculated. Separate 
values have been found for each of the terms in 
d*s. Terms of the same kind in d* lead to energy 
sums in d's rather than to the separate energies of 
each term. There are three cases, 2D, ?F, and 2G, 
in d° for which more than one term of the same 
kind arises. When an s electron is added, singlet 
and triplet terms in d*s are obtained from these 
terms of d*. Separate energies have not been 
found for these singlets and triplets. Only the 
sum of the energies of those terms of the same 
kind is obtained. 
The formulas for the energies in terms of the 


that the calculations were made. Now that the 
analysis of V II has been completed by Meggers 
and Moore,' it seemed desirable to present the 
theoretical calculations also. The septet and 
quintet terms of Mn II have all been found by 
Curtis,? and the theory gives in addition the 
positions of the singlets and triplets. To aid in 
the analysis further, Ostrofsky’s* formulas were 
used to calculate the missing terms in V II 3d 
with parameters found from the observed values 
of V II 3d‘. 

The terms in Russell-Saunders coupling can be 
calculated by the standard method given by 
Slater.4 Van Vleck® has given a different method 


1W. F. Meggers and C. E. Moore, J. Research Nat. Bur. 


Stand. 25, 83 (1940). 


2C. W. Curtis, Phys. Rev. 53, 474 (1938). I am grateful 


to Dr. Curtis for permission to include the unpublished °F 


Slater-Condon Fs and Gs were then used to find _ term. 


the term values in V IT 3d*4s and Mn II 3d*4s. 
It was to assist in the analysis of these spectra 


TABLE I. Term values of d*s for V II. 


3M. Ostrofsky, Phys. Rev. 46, 604 (1934). 
4J. C. Slater, Phys. Rev. 34, 1293 (1929). 
5 J. H. Van Vleck, Phys. Rev. 45, 405 (1934). 


“AP = Fy—6F:—12Fy 

3P = 

= 

5P = Fy—147F,—3G2 

1D = Fo t+ (193 Fo? — 1650 F 
8) = Fy +5 F.+3Fy4 (193 1650 +8325 — 
1D = Fot+ Fy— (193 1650 +8325 F 2)! 
3D = Fot+5F.+3 Fy— (193 1650 F2 +8325 — 
1F= F\+9F,—87F, 

3F= 

3F = 

5F= 

1G = Fo—11F2.4+13F 

3G = Fo— 11 F2+13F,—2G, 

= Fy— 6F,.—12F, 

= Fo— 6F.—12F,— 2G, 


Fo=3Fo(nd, nd) +3Fo(n’s, nd) =31155.58, F(nd, 
F,= F,(nd, nd) =90.31642, G:=G.(n's, nd) = 1335.338 


WAVE NUMBERS 


22273.54 22898 

20237.81 20228 

19147.03 19214 

13646.67 13873 

— 50274 

2G, 44139.29 47603 
25191.08 24535 

2G, 20600.84 21865 
34228.79 34058 

30637 .30 31388 

8903.58 8054 

2926.16 2713 

19112.93 19178 

16446.05 16508 

23391.09 22898 


20302.57 20228 


nd) =1195.557 
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TABLE II. Separations of the levels of the terms of d*s for 
VII, given in wave numbers. The subscript on the left of the 
terms tndicates the number of terms of the same kind which 
occur. The sum of the separations of the two *P, *D and *F 
terms are given. {na= 154. 


Oss. CALc. Oss. CAaLc. 


5G 


“15/20fnq «111.49 116 
4 
12/20tnd 


3 


5 
«194.58 193 
“4/4tnad 159.46, 154 
«121.75 116 


2 | 

2/4tna $2.19 77 4/3tnd 228.10 

3/3tnd 170.33 154 
3/4ind 146.88 116 2D, 

na 83.02 77 | : 3/6ind —55.03 77 

2/6tnd 53.37 51 
6/10fna 83.03 92 

«37.870 77s 219.94 231 

3/4tnd —62.16 116 
| 


TABLE III. Term values in V Il 3d* from Ostrofsky's 
formulas. The plus sign written as a superscript on the right 
indicates the higher of two terms of the same kind, the minus 
sign the lower. 


WAVE NUMBERS WAVE NUMBERS 

Oss. CaLc. Oss. CALc. 
1st _ 61919 1G 17704.72 18560 
44451.73 44843 Wy 18985.24 18188 
IGt 36218.81 35589 3p 18110.74 17854 
apr 30096.11 29842 1G 14365.92 14399 
apr 31962.80 29765 3p- 11502.76 13847 
if 26633.56 25590 3F- 13352.50 13769 
ip- 20774.66 23124 3H 12427.79 11932 
1s- 19696.34 20208 5p 0.00 —580 


Fo =5318.86, F2 =869.862, and F4=75.2701. 


of arriving at the same results, and in addition 
shows how to find the terms of a‘s from those 
of a*. Thus d*s may be most simply found from 
d’, which has already been calculated® by the 
methods of both Slater and Van Vleck. These 
calculations included the separation of the two 
*D terms. The separations of the different levels 
of d*s are included as they were used in identify- 
ing the *F terms. These levels are calculated from 
the diagonal matrix elements of spin-orbit inter- 
action by Goudsmit’s method.’ The nondiagonal 
elements of spin-orbit interaction between the 
different terms in Russell-Saunders coupling are 
neglected. 

Similarly d's may be found from d°. In order 
to find the terms of d®° from those given for d°s in 
Table IV, the G2(n’s, nd) parameters which 


°C. W. Ufford and G. H. Shortley, Phys. Rev. 42, 167 
(1932). Robert Serber, Phys. Rev. 45, 461 (1934). 
7S. Goudsmit, Phys. Rev. 31, 946 (1928). 
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TABLE IV. Term values of ds applied to Mu Il 3d*4s. 
The subscript on the left of the terms indicates the number of 
terms of the same kind which occur. Where more than one 
term of the same kind is indicated, the formulas represent the 
sum of the energies of the terms of the same kind. 


WAVE NUMBERS 


Ons. Cale. 

3] =Fo —24F 2 —90F 4 —3G2 38711 
= Fo —24F2—90F 4 41884 
3H = Fo —22F —30F 4 46913 
1H =Fy—22F2—30F 4 50086 
9G =Fo—25F 2 —-190F 4 27571 26620 
3G = Fo —25F 2 -190F 4 3296S 
2G =2Fo —10F —300F 4 117970 
=2Fo —10F 2 —300F 4 —2G2 124315 
5F = Fo —13F2—180F'4 —4G2 43475 43806 
=Fo—13F 2-180) 4 50151 
oF =2Fo—34F 2 —180F 96406 
2'F =2Fo—34F 2 — 18094 102751 
=Fo—18F2—225)' 4 —4G2 32828 32909 
3D) = Fo —-18F2—225F 4 39254 
#D =3Fo —10F2 —300F 4 —9G2 197459 
=3Fo —10F2—300F 206977 
=Fy—28F 2 —105F' 4 29912 30329 
36674 
=Fo+20F 2 —240F 84650 
1P =Fo+20F 2 —240F 87822 
7S =Fo —35F 2-315 F 4 0.00 -—37 
§S = Fo —35F 2 —315F 44+Ge 9472.86 9555 
38 = Fo —3F 2 —-195F 4 —3G2 57580 


= Fo —3F2—195F 4 —G2 60752 
Fo =10F nd) 4+-SFo(n’s, nd) =84173.83 


Fe =F nd) = 1355.969 
Fa =Fa(nd, nd) =91.49416, G2 =Go(n's, nd) =1586.311 


represent the interaction of the s electron with 
the electrons must be omitted. If Gz is omitted, 
pairs of equivalent terms are obtained, each pair 
representing a term of d° with the same value of 
L and a multiplicity, 2S+1, midway between 
the multiplicities of the two terms of d°s from 
which the pair arose. For example, d°*J= Fy 
—24F,—90F,;, where now Fo=10Fo(nd, nd). 
Laporte and Mack® have carried out similar 
calculations for d°s in connection with their work 
on tungsten. 

The results for the term values of d’s in terms 
of the Slater-Condon® Fs and Gs together with 
the observed and calculated values for V IT are 
given in Table I. Table IT gives the separations 
for d's. The agreement is fair except for the 
upper separation of the *Ds and the lower 
separation of the *Ps which were not used to 
calculate {,a. The assignments of the *Fs on this 
basis seem satisfactory. The observed and calcu- 


( §Otto Laporte and J. E. Mack, Phys. Rev. 58, 188A 
1940). 
® For a definition of these radial integrals and a complete 
account of the theory cf. Condon and Shortley, The Theory 
of Atomic Spectra (Cambridge University Press, 1935) 
hapters VI and VII. 
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lated terms for VII 3d* from Ostrofsky’s 
formulas are given in Table III. The parameters 
were found from */7, *G, *D, WU, and 'F as these 
terms were considered the most reliable. Finally 
Table IV gives the formulas for the term values 
of d's and the observed and calculated values 
from Mn II. 

Although the agreement of the observed and 
calculated term values is rough in terms of the 


accuracy with which spectra may be observed, 
some idea of the relative positions of the terms is 
a help in analysis as Meggers and Moore! have 
pointed out. As a further aid in the analysis of 
Mn II, it would be advisable to separate the 
different terms of the same kind in d° by the 
methods® used for d*. I am grateful to Professor 
C. W. Ufford for suggesting this problem and for 
his assistance in solving it. 
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Magnetization Near Saturation in Polycrystalline Ferromagnets 


T. Hotstern, New York University, New York, New York 


AND 


H. Prmakorr, Queens College, Flushing, New York 
(Received December 16, 1940) 


An important mechanism involved in the variation of 
the magnetization of ferromagnets near saturation is the 
rotation of the magnetization vector under the combined 
influence of the magnetic field and crystalline anisotropy 
torques. The effect of this mechanism in polycrystalline 
specimens has been calculated by Akulov and Gans; their 
result is Ma,= Mo[1—(c/M.*)H~*], where H is the sum 
of the external and demagnetizing fields, Mo the saturation 
magnetization, and ¢ a constant proportional to the square 
of the crystalline anisotropy constant. The Akulov-Gans 


I. INTRODUCTION 


CCORDING to present day theories! of the 
magnetization curve of a polycrystalline 
ferromagnetic in high fields, an important con- 
tribution to the deviation of the observed mag- 
netization, |M|,, from its saturation value 
|M.o|, arises from the crystalline anisotropy 
torques. These torques, acting on the magnetiza- 
tion vector of each crystal grain, tend to rotate 
this vector away from the direction of the 
external magnetic field, and toward the direction 
of the nearest easy axis of magnetization within 
the grain in question. Thus, in the equilibrium 
position, the direction of the magnetization 


1N.S. Akulov, Zeits. f. Physik 69, 822 (1931); R. Gans, 
Ann. d. Physik 15, 28 (1932). Cf. also, E. C. Stoner, 
Magnetism and Matter (Methuen, London, 1934), p. 405; 
F. Bitter, Introduction to Ferromagnetism (McGraw-Hill, 
New York, 1937), p. 222; R. Becker and W. Doring, 
Ferromagnetismus (Springer, Berlin, 1939), p. 167. 


derivation, however, is subject to a serious error; namely, 
neglect of the internal magnetic field arising from the 
magnetization itself. In the present paper, this internal 
field is taken into account; then, with the dual assumption, 
of randomness of orientation of crystallographic axes, and 
irregularity of shapes of the individual crystal grains, one 
obtains the formula, Ma, = Mo[1—(c’/.M.?)H~*]. Here, c’ is 
a slowly varying function of H; for H>>4rMo, c’=c; 
H<4nMo, c’=}3c. Applications of the last formula to the 
analysis of the experimental data are discussed. 


vector in any crystal grain is determined by the 
balance of the crystalline and magnetic torques 
acting on it, viz.: 


T=M,)X(H+h), (1) 


where My is the magnetization of the crystal 
grain under consideration; H the external mag- 
netic field (the demagnetizing field due to the 
magnetization on the surface of the polycrystal- 
line specimen is incorporated into H); h the 
internal magnetic field, i.c., the magnetic field 
due to the magnetization vectors of the totality 
of crystal grains; T the torque due to the crystal- 
line anisotropy forces within the crystal grain.! 
T in any grain, depends on the direction of Mo 
with respect to the crystallographic axes of that 
grain, and is independent of position within any 
given single grain.' The magnitude of T is 
proportional to the “anisotropy constant” K.! 
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It is obvious that h and Mo cannot be assumed 
constant even within a single crystal grain and 
must hence be considered as functions of position 
within the polycrystalline specimen. However, 
only the direction of Mo varies with position, 
since the magnitude of Mo is the intrinsic domain 
magnetization and is therefore constant. 

The problem of calculating the deviation of 
|M| from Mo| as a function of H was first 
considered by Akulov and by Gans! on the basis 
of the general mechanism outlined above. These 
authors obtained the result: 


BK? 


with B=0.0872 for Fe; 0.0762 for Ni. 

However, Akulov and Gans, in their derivation 
of (2), neglected completely the internal mag- 
netic field h, and thus the torque term Mo Xh in 
Eq. (1). The neglect of this torque term cannot 
be considered satisfactory since it is of the same 
order of magnitude as the term My) XH, for the 
range of external fields H used in the experi- 
ments.? It is not possible to dismiss the term 
with the observation that |h|<|H|, 
since H and M, are almost parallel, while h has 
a direction which in general, is quite different 
from that of My and H. 

In this paper, the deviation of |M| from 
|Mo|, considered as a function of H, will be 
obtained taking the internal magnetic field h 
into account. 


II. DERIVATION OF EXPRESSION FOR THE 
DEVIATION OF FROM |Mo| 


It is convenient to resolve My into components 
Mand m, parallel and perpendicular, respectively, 


to H. Thus, 


and Eq. (1) becomes 
T=mXH+MxXh+mxh, (4) 


while the observed magnetization |M|j is given 
by the integral of |M!| over the volume of the 


2? These experiments are discussed in Section IV; the 
magnitude of the fields used is ~2000—4000 gauss. Cf. 
reference 5. 


polycrystalline specimen, i.e., 


1 |M)| 
V Jy V vy |Mo|*/7 (5) 


Since only the magnetization observed at high 
external fields’ is under consideration, the devia- 
tion of |M|,, from |Mo! is small, and the follow- 
ing approximations may be considered as holding 
at each point within the specimen : 


{m|/|Mo|<1, (6) 
|M|—|Mo| |Mo—m|—|Mo|_ |m|* 7) 
|Mo| |M,| |Mo|? 


From (6), (7), it follows that with neglect of 
terms =|m|*/|Mo/? or higher, M is constant in 
magnitude, as well as in direction. 


Further, 
(8) 
since 
—div My 
h=-—grad | — dr’ 
—divm 
f —dr’, 
|r—r’| 


and hence h=m. 
The approximations (6)—(8) applied to Eqs. 
(5), (4) vield: 


1 
2|Mo|?V 


T=mXH+Mxh. (10) 


Further, in view of (6)-(8), T in any crystal 
grain can be taken as dependent only on the 
direction of the vector M or H with respect to 
the crystallographic axes of that grain, since My 
may be taken parallel to M to the first order in 
|m|/{|Mo!. T is hence independent of m. 

Thus to find |M|,, it remains to calculate 
from Eq. (10) m as a function of position within 
the specimen, and then to evaluate the integral 
in (9). This procedure will now be carried out. 

First, taking the vector product of both sides 


* These fields are such that |M| is well above the 
“knee” of the magnetization curve; e.g., H>1000 gauss 
for Fe. 
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of Eq. (10) with H, and solving for m, one My = Fey+ (42) ahxy, (18b) 
obtains, m,:=0, (18¢) 
-1 — (h. 2 
m=F+(4r)~'a[ h—(h-H)(H/|H}*)] (11) k-h, = —4rk-m,, (19) 
with the abbreviations kXh, =0. (20) 


F=(HXT)/|H)’, (12) 
a=4r|M (13) 


One also has the general field equations re- 
lating h and m, 


div h= —42 div divm, (14) 


In Eqs. (18), (19), (20), Fx, hy are the vector 
Fourier components of F and h, similar to m, 
being the vector Fourier component of m in 
(16); are the x components of Fy, hy; 
the z axis has been chosen parallel to H. 

To proceed further one eliminates hy, from 


(18), (19), (20), finding m, as a function of k. 


curl h=0. (15) One has as the solution of (19), (20) for hy, 
Introducing the Fourier resolution’ of m, h, = —4enn-m,, (21) 
m= {m, exp (2rik-r)dk, (16) with n=k/|k’, 
whence, substituting (21) into (18), and solving 


one obtains from (9) for Mxz, My, One obtains 


1 Fy2(1+an,?)+ Fyyan.n, (22) 
while from (11), (14), (15) the x, y, s components = Fry(1+an,’) 
of m, satisfy the relations: 1+a(n,?+n,’) 


Finally, substituting (22) into (17), one obtains 


M|,=|M | 
| | | al( (1+a(n,2+n,*))? 


(1+a(n,?+n,7))* 


| Fy: | *dk 


(1+a(n,*+n,7))? 


Now it will be shown in Section III below, that as far as the integrations over k in Eq. (23) are 
concerned, one may effectively take 


| =| Fey|?=3| Fel, (24a) 
| =0, (24b) 

and take | F,|? to depend only on the magnitude of k, i-e., 
| (24c) 


( De. , Fourier analysis method has been used by W. F. Brown, Jr., in a related problem, cf. Phys. Rev. 58, 736 
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Thus, 


1 
| w= | Mo| | (|k!) | k *d|k| “ff in 


1 1 
= f (1+ ) sin ) 
8|Mo\2V 0 (1+ sin* 


||T sin (T, H) | (1+a)!+a! 
=|Mo| (:- log ( “))): (25) 
(1+a)!—a! 


Equation (25) is the desired expression for the functional dependence of |M 4,—{!Mo| on H. In 
the limiting cases of a<1 and >1, Eq. (25) reduces to the following simple forms: 


{|T sin (T, 
H-?}, when i.e., 441 (26) 
0 
sin (T, 
(: 1), when a>>1, i.e., 44. (27) 
| ~ 


Equation (26) may be derived by neglecting the term proportional to a in Eq. (11), or equivalently, 
by neglecting M Xh in Eq. (10), or MyXh in Eq. (1). Equation (26) is thus identical with Eq. (2) 
derived by Akulov and Gans. Thus it is seen that the Eq. (2) or (26), of Akulov and Gans is only 
applicable, if H>>42M/)= 21,800 gauss in Fe, 6400 gauss in Ni. It should further be emphasized that 
only in the limiting cases 44. 4% Mo<H, does |Mo| vary, as const. /7~ (with different 
constants for the two limits). For intermediate field strengths, the dependence of |M\4 on H is 
much more involved, being given by Eq. (25). 


III. JusTIFICATION OF EQuATIONS (24) constant value of F in V; as F;. Thus, 


It now remains to justify Eqs. (24). One must 
consider the evaluation of the integrals in Eq. Fe=L Fi] exp (2aik-r)dr 
(23). These integrals are of the type: 7 Vi 


=> F; exp (27ik-R,) | exp [27ik- (r—R,) |dr 
i 


f W.(n) Feedk, f W,(n) Faydk, 


=) F; exp (2aik-R;) V;W;(k), (30) 
f Fee (28) 
where R; is a vector from the center of the poly- 
crystalline specimen to the center of the jth 
crystal grain, and W;(k) is the spatial average 
of exp [2rik-(r—R,) ] over Vj. 
The perfect randomness of the distribution of 
F(r) exp (2mik-r)dr, (29) the orientations of the crystallographic axes of 
Ss the various crystal grains within the polycrystal- 
where the volume integrals in the sum over j _ line specimen leads to the following ‘‘randomness 
extend over V;, the volume of the jth crystal properties’ of the vectors F;: 
grain. In V;, T and hence F, are constant both (A) There is no correlation between the magni- 
in magnitude and direction; one writes the tudes and directions of the vectors F; and 


Now 


b) | 
9) 
0) 
Ny 
in 
ik 
m 
1) 
| 
_| 
_| 


ie., 
Fi] ~ idm) (| J}, (31) 


where V/(V)y is the number of terms in the 
sum over 

(B) There is no correlation between the magni- 
tude and direction of F; and the position of the 
jth crystal grain within the specimen, i.e., 


where f(F;) is any function of F;, g(Rj) is any 
function of R;, and 


| w= (V/(Vi)w)* g(Ri) 


(C) There is no correlation between the direc- 
tion of F; and any fixed direction within the 
specimen, i.e., 


Fiy=3| 
or 
(Fiz)w = (Fiy)w = (3 | (33) 


From Eq. (30) one sees that Fy, oscillates 
rapidly when k varies over a region V~', since 
R; in exp (27ik-R,;) is ~ V!. On the other hand, 
from Eqs. (30), (32), (31), 


| Fe| ~(\exp R;) V;W;(K) | wv 


VjWills) | (34) 


and thus the variation with k of the order of 
magnitude of F, is independent of the rapid 
oscillations due to exp (277k-R,) and is deter- 
mined only by the W,(k). Also, from their 
definition, the W,(k) vary appreciably with k 
only when k ranges over a region ~V;~'; in 
particular, W,(k)=21 if |k|=V,-!, and W,(k) 
decreases if |k| exceeds V,~!. 
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Thus, turning back to the integrals of Eq. (28), 
one sees that the effective region of integration 
in k space is defined by the condition k =V;~!. 
Now W depends only on n=k, kK); one calcu- 
lates the variation in n, An, corresponding to a 
variation in k, Ak, and finds 


An=[k (Ak Xk) ]/ |k|* 


over most of the effective region of integration. 
Therefore, Y varies appreciably only when k 
varies over a region V;~'. 

Thus, it is possible to divide the k space into 
cells of size 2, which, if V'<2<V;~' possess the 
following properties: (1) Many oscillations of 
F,, | Fx|2, etc., take place in @. (2) The order of 
magnitude of | F,|? and of the value of VY, remain 
effectively constant over ©. 

It therefore follows that the integrals (28) can 
be written as 


f W,(k,/|k,|) Fidk 
2, 


Jo 


and hence the Fre in (28) are effectively replaced 
by their ‘‘smoothed over”’ values: 


F..dk’. 
Q 


2 2 
Now the “smoothed over” values of Fy:, 
x, satisfy, as is shown immediately below 
the following relations: 


of = | "| Fi (%dk, (36a) 
2 2 


| Fur Faydk=0, (36b) 
and thus, by use of Eqs. (35), (36a), (36b), the 
validity of Eqs. (24a), (24b), is established. 
To prove Eq. (36a), for example, one has: 


| FV}? [Wi(k) Pak +27 f exp (2mik- Rj.) Wi(k)dk 
Q if 2 


Fin W;(kKo) Fy2F iz V; W;(Ko) W(Ko) exp R;,) 
7 if 


XK sin (wAk,X sin (wAk,Y sin (wAk.Z j:) (37) 


in: 


th 
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where Kp is a vector from the origin of k space 
to the center of 2 and Ak,Ak,Ak. =. The last 
(approximate) equality in Eq. (37) is valid, 
since in 2 each W,(k) is approximately constant. 

From Eq. (32), the order of magnitude of the 
single sum is: 


~(VALW Fie 
(38) 
while the order of magnitude of the double sum, 
ViVi Wj (Ko) Wilko) 
Xexp (2riko: Rj) | iz, 
Vj ViW;(Ko) Wilko) 
Xexp (2miko- | F;z)?, 


~(QV)-\| (Findon (39) 


the last (order of magnitude) equality being a 
consequence of Eq. (31). 

Thus, from (38) and (39) the ratio of the 
orders of magnitude of the double sum and the 
single sum is (VQ2)-'<1. Hence, 


2 7 
and similarly 
2 


“Randomness properties” (B) and (C) now 
insure the equality of the right-hand sides of 
Eqs. (40) and (41), and Eq. (36a) follows 
immediately. 

To prove Eq. (36b), one observes that 
foFx2Fxydk can be expressed as and 


7 ig! 


that from Eq. (32), 


(42) 


the last (order of magnitude) equality being a 
consequence of “randomness properties” (A) 
and (C). 


Thus the ratio of the orders of magnitude of 
the single sum in Q-'foFy:Fy,dk, and of the 
single sum in 0-1 fy Fy .dk, is (cf. Eqs. (38) and 


(42)), ~(V / VIKA. Further, by an argument 
similar to that given in connection with Eq. 
(39), it can be shown that the double sum in 
foFx:Fx,dk is negligible compared with the 
single sum. Therefore, in comparison to the in- 


2 
tegral Q-' f)Fy.dk the integral Q-' fo Fy. Fy,dk 
can be taken to be equal to zero (to within 


terms Vii, and (VQ)-'<1). 

It now remains to justify Eq. (24c). For this 
purpose it is necessary to make use of another 
‘randomness property”; viz.: The shapes of the 
individual crystal grains are not correlated in 
any way with a fixed set of directions in the 
crystal. Then, apart from the rapid oscillations 
of |F,|* with k which are smoothed out by the 
integrations over 2 in Eq. (35). or (36) or (23), 


| Fx|?= 


Ff exp (2rik-r)dr 
i V; 


can only depend on |k|. Hence the validity of 
Eq. (24c). 


IV. Discussion OF EXPERIMENTS 


The experimentally observed variation’ of the 
magnetization of polycrystalline ferromagnetics 
with external magnetic field at high fields is 
customarily described by an empirical relation 
of the form :® 


d|M | obs 
dil 


=All’*+BH+*+CH, (43) 


where A, B, C, are assumed to be independent 
of 77. 

It has further been observed that C is inde- 
pendent of the metallurgical history of the 
specimen ;> it therefore must arise from the 
variation of the intrinsic domain magnetization 
with the field. The magnitude, temperature de- 
pendence, and actual slow variation with field, 


5 A. R. Kaufmann, Phys. Rev. 55, 1142 (1939); 57, 1089 
(1939). H. Polley, Ann. d. Physik 36, 625 (1939). Only Ni 
and Fe but not Co, were investigated. 
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of C, have been treated in a recent paper by 
the present authors.*® 

A and B depend on the metallurgical history 
of the polycrystalline specimen ;' A is conditioned 
by the plastic deformation of the material,® 
(small for well-annealed specimens) and is tem- 
perature independent,® whereas, B is determined 
by the crystalline properties and the elastic 
state,® and is furthermore strongly temperature 
dependent.® The theoretical interpretation of the 
term AJI- is still in doubt; however, a recent 
paper by W. F. Brown, Jr.,‘ ascribes this term 
to the existence of anisotropy torques which are 
concentrated along lines of dislocation produced 
by plastic deformation. 

Finally, the term BH- arises from the mecha- 
nism discussed in this paper, i.e., from the 
rotation of the magnetization vector away from 
the direction of JZ and toward the various 
directions of easy magnetization (in the various 
individual crystal grains) by crystalline aniso- 
tropy torques constant within any one single 
grain. Comparing the identical Eqs. (2) and (26), 
one sees that 


K=[{|Tsin (44) 


whence from Eqs. (26), (27), d|M)},,/dH may 
be written in the form BH- in the two limiting 
cases, i.e., 


when (45) 

when (46) 


For intermediate values of H, d|M|,/d/J is 
obtained by a differentiation in Eq. (25), and 
can be written in the form BH~ only if B itself 
is considered as a (slowly varying) function of 77. 

Now in all the experiments so far performed, 
H was considerably less than 47, so that, in 
evaluating K from the experiments, one should 


6 T. Holstein and H. Primakoff, Phys. Rev. 58, 1098 
(1940). 


use the relation :’ 


K=[8"|M,|? 
X (observed coefficient of H)}'. (47) 


On the other hand, these experiments have 
previously been interpreted in terms of the 
Akulov-Gans equation (2) or (26) or (45); thus 
the values of K, previously deduced from the 
existing experimental data, are smaller than they 
should be by a factor equal to 2!. If these 
previously deduced values of K be all multiplied 
by 2, they agree somewhat better than before 
with the values of K obtained directly from the 
magnetization curves of single crystal grains; 
cf. especially the experimental data given in 
Polley.’ 

The analysis of the present paper can be 
tested by experiments on the magnetization 
curves of polycrystalline specimens, in which, 
both the limiting cases, 44.Mo>H, 44M 
and also the case of intermediate H (described 
by the complete Eq. (25)) are studied. Now, 
both for iron and nickel, K is so small, that 
unless /7 is considerably less than 42.WZo, the 
whole term BH is negligible. For cobalt, how- 
ever, the crystalline anisotropy energy and 
hence K, are so large, that the portion of 
d|M|.»;/dH arising from the mechanism under 
discussion, could easily be observed as a function 
of H for values of H7 less than, of the order of, 
and considerably greater than, 47JJo.° In this 
way, a complete experimental test of the Eqs. 
(25), (26), (27) can be obtained. 


7 This formula, as well as the rest of the analysis in the 
present paper and in the derivation of the Akulov-Gans 
equation (2), is based on the assumption that the crystalline 
anisotropy torques are just those present in a strain-free 
material. The existence of strains creates additional 
anisotropy torques, which, if constant over regions com- 

arable in size to a single crystal grain (if not, cf. W. F. 

rown, Jr., reference 4) can be easily incorporated into the 
present treatment. The only resulting difference in the 
final formulae is that K in Eq. (44) and hence in Eqs. 
(45), (46), (47) has to be amplified to include the “‘strain 
anisotropy” constant, and that the numerical value of 8 
must be changed somewhat. 

5 Reference 5, p. 643, Fig. 12. 

®Since cobalt has a hexagonal lattice, the numerical 
value of 8 for Co will differ from the values quoted in Eq. 
(2), which are only valid for the cubic lattices of Ni and Fe. 
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Phase Series 


E. J. HELLUND* 
(Received October 18, 1940) 


An expansion theorem for the calculation of the parameters occurring in the scattering cross 
sections has been developed. The series converges exponentially and is expressible in terms of 
known functions, provided the interactions are expressed as a power series in r and 1/r. 


INTRODUCTION 


HE solution of problems involving the 

scattering of one particle by another is 
attended by extreme difficulties. In general, 
evaluation of cross sections can be obtained only 
in infinite series, each term of such a series con- 
taining a parameter (phase), which must be 
determined from a linear differential equation of 
the second order. The evaluation of this param- 
eter must proceed, usually, by numerical integra- 
tion of this equation. Such calculations are 
excessively laborious and are, furthermore, sub- 
ject to numerous errors. It has seemed desirable, 
therefore, to derive an expansion for the phases, 
in series form, the convergence of which could be 
accurately determined. Such an expansion pos- 
sesses a particular advantage which lies in the 
latitude permissible in the choice of function to 
represent the interaction between particles. This 
is of special importance in the evaluation of 
integrals involving the interaction potential and 
other known functions. Such integrals invariably 
arise, as the distortion of the wave function of 
one particle by another involves the values of 
the interaction for all distances of separation. 
Apart from the points mentioned, however, a 
series solution holds little practical superiority 
over the numerical integration of the differential 
equation. 


1. First approximations 


The problem to be solved is, then, the determi- 
nation of J(@)dw, the number of particles scat- 
tered into a given solid angle dw per unit time, 
when the incident beam is such that one particle 
crosses unit area per unit time. The series ex- 


ona Now with the Battelle Memorial Institute, Columbus, 
io. 


pansion for /(@) then has the form :'* 


1(8) = | f(0)|*, (1a) 


where 


f(8) =(1/2ik) (2/+1) 
XLexp (276,) —1]P.(cos (@)), (1b) 


hk? =82°mE/h?, (1c) 
and 6, is a constant to be determined from the 
equation, 
d*(rp) 8r*mV(r) 
dr? h? r? 


One selects that solution finite at the origin 
(or that solution of lower order singularity at 
the origin,’ and determines 6, from its behavior 
at infinity. 6, is then defined by the asymptotic 


solution, 
y~(1/r) sin (kr (3) 


The determination of 6, is then the immediate 
point of attack. Various formulae have been 
given for 6, which are due, in order, to Born,‘ 
Jeffreys,> and Massey and Mohr.*® 


2 
V(r) J (kr)rdr, (4) 
h? Jo +4 


1H. Faxén and J. Holtsmark, Zeits. f. Physik 45, 307 
(1927): 

2N. F. Mott and H. S. W. Massey, Theory of Atomic 
Collisions (Clarendon Press, Oxford, 1933), first edition, 

24 


3 Reference 2, p. 30. q 
* Reference 2, p. 28. 
5 Jeffreys, Proc. London Math. Soc., Series 2, 23, Part 6; 
or reference 2, p. 92. 
6H.S. W. Massey and C. B. O. Mohr, Proc. Roy. Soc. 
A144, 202 (1934). | 
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where 
F(r) = (822m /h?)(E— V(r)) —l1+1)/r?. 


Formula (4) is valid for small phases, that of 
(5) for large phases and for (6) Massey and Mohr 
claim accuracy for both large and small phases. 
Actually, all three expressions must be regarded 
as first approximations, not valid generally for 
large phases. An exact solution of (2) in finite 
terms has never been given, except for special 
choice of V(r). Moreover, any expression of the 
solution of (2) by a simple integral of some func- 
tion of V(r), may be seen by actual substitution 
to solve this equation for V(r) replaced by a 
function of V(r) and dV(r)/dr. The accuracy of 
(5) and (6) for large phases is then dependent on 
the choice of V(r), and unless one has some inde- 
pendent criterion for their accuracy, their use in 
calculations is somewhat hazardous. It is possible, 
however, to give a series for 6;, the first term of 
which is the expression (4). 


2. Derivation of the first-order equations 


Let us suppose, considering one dimension 
only, that it be required to determine the con- 
ditions necessary that a function f of variables 
r, &, n satisfy Eq. (2). € and 7 are functions of y 
so chosen that 


f(r, & (7) 
af(r, &, n)/dr=9(rp)/dr (8) 


for r=r, and r=r,;+Ar;. » must then possess 
first derivatives defined by the following system :7 
of dg dy 
de dr On > 


dé d 
dr Ondr dr 


where /7,(r) has been set for the coefficient of (7) 
in (2). In particular, if one takes f as a linear 
function of — and n, 


f=tgitnes, (11) 


7This procedure has some resemblance to that of 
A. Schuchowitzky and M. Olewsky, Physik. Zeits. 
Sowjetunion, 11.5 pp. 498-512 (1937). They however, use 
finite intervals, and after determining some parameters so 
as to assure continuity of the wave function, the remainder 
are evaluated by use of the variational principle. 
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where ¢; and ¢2 both satisfy the condition 
¢/dr = — (H,)(¢), (12) 
one has the system 
dt/dr ¢2°n) (Hy —H.)/W(¢) =0, (13) 
given) (Hy —H,)/W(e) =0, (14) 
where W(¢) indicates the Wronskian of ¢, ¢o, 
W(¢) = ¢1(0¢2/0r) — (15) 


The Wronskian, under the condition (12), is a 
constant for all values of r°. 

It is a simple matter to confirm that (11) 
does actually satisfy (2) provided that also (12) 
and (13) are satisfied. 

Equations (13) and (14) may, for convenience, 
be written in matrix form. 


where &, and 7, are arbitrary values assigned to 
and » for r=c. 


3. Solution of the integral equation 

Equation (16) is a Volterra integral equation, 
of the second kind, for the unknowns é and 7. 
If one denotes the unknown matrix by y, the 
constant matrix by g(c), the solution is given by 
the series :° 


y= 8 +2 Kjg(c), (17) 
where K is so defined that 


f K(r)ar; f 
xf (18) 


and K denotes the four-component matrix 
of (16). 


8 The derivative vanishes identically. 

® Gerhard Kowalewski, Integral Gleichungen (Walter de 
Gruyter and Co., Berlin, 1930), first edition, pp. 49-90. See 
also E. G. C. Poole, Theory of Linear Differential Equations 
(Oxford Clarendon Press, 1936), Chapter I, for discussions 
on the convergence and uniqueness of the expansion. 
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One may easily demonstrate the convergence 
of (17) as follows: If ¢: and ¢g» are finite in the 
interval r to c, let us replace each element of the 
matrix by where exceeds both g; and 
over the entire interval. Also, let us take the 
absolute value of /7,—//7, at every point in the 
interval. One then obtains a series, term by 
term larger than (17), of the form, 


— 
W(¢) | 
1\ f/é-+n- 
)( ). (19) 
1 2 


The series (17) thus is absolutely conver- 
gent and summable in any order," provided 
|\(H,)—(H,)| is integrable in the interval c, r 
The series likewise will still converge as r 
approaches infinite values, provided the integral 
converges. 


4. Differential equation for 5, 


In order to evaluate the series (17) it is con- 
venient to use a vector representation instead of 
a matrix representation. One then considers the 
pair as a two-component vector and the matrix of 
four elements as a product (dyadic) of two two- 
component vectors, which are 


¢=(¢1, ¢2); S=(¢2, (20), (21) 


g and ¢ are thus mutually perpendicular. 
If one denotes the vector &, 7 by P, then 
corresponding to (13) and (14), one has, 


dP (22) 


where P denotes the scalar product ¢27. 
One may set P equal to the product of a 
scalar and a unit vector of components cos y 


5.. Explicit series for 5, 


and sin y. The resulting equation for ¥ is, 


dy + 1, — H,) 
— sin? (y+a)=0, (23) 
dr W(¢) 


where tan a= ¢)/¢2. 

y may be identified with the phase angle 6, 
depending on asymptotic behavior of ¢; and ¢g» 
atinfinity. Asa particular example, if ¢;=sin (kr), 
ge=cos (kr), then = 

From (23), one may derive the perturbation 
formula 


Y=Yo—exp fexp 


where 


f —H,) sin 2(yo+a)dr; 


G2= sin? (yo+a) 
W(¢) 
by setting y equal to yo+(y—vo), where yo 
satisfies (23) for the potential V(r). The Born 
formula (4) is a less accurate expression, ob- 
tained by equating yo to zero, and retaining only 
term linear in V. Another perturbation formula 
could be obtained by using the Picard second 
approximation" on (23). 

One may also formulate the eigenvalue problem 
from (23), for 6; by requiring the solution of the 
following equation to approach zero as r ap- 
proaches infinity : 
dy 1/8x°*m l(l+1) 
V(r)+ ) 
dr I? 


l 
Xsin? =0. (25) 


It is not the only solution as the addition of 
any multiple of z to the solution of (25) is also 
an acceptable solution. 


Using (22), one may write (17) in the following manner, 


"E.G. C. Poole, reference 9. 


“E. T. Whittaker and G. N. Watson, Modern Analysis (Cambridge University Press, 1935), fourth edition, pp. 18-28. 
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where I'(s) =[H,(r.) —H7,(r.) ]/W(¢) and the following inequalities hold, 
(27a) 


The ratio of — to 7 then determines the phase 6;. In the event & and , are not specified, and ¢; is 
bounded at the origin, & and , are replaced by 1 and 0, and ¢ is taken to be zero. This procedure is 
permissible, however, only when the series still converges. In general this will hold provided 17, — HZ, 
is everywhere finite, and decreases faster than 1/r as r approaches infinity.” 

If one employs the functions r!J (kr) and riJ (kr) for ¢: and ¢» in (26) and (27), the Wronskian 


of ¢ becomes (—1)!'*'(2/m), and 7,—H, becomes —8x?mI'(r)/h?. To terms of the second order, 
then, one has for 6,, with c=0, &)=1, no =0, 


16 


J(kr)r V(r)dr+(—1)! f f J(kr)J(kr)rV(r)dr 
h? Jo ht o Jo 4h 


1672m?2 x r2 
dr (J (kre) J (Rri) — (kre) J (Rr) (kre) (rari V (re) Viraddri}. (28) 


In general, for the choice of functions given above, one may write for tan 6), the following quotient, 


+4 
tan 6,=(—1)! (29) 


1+(-1)! J(kr)J(kr)r V(r)dr 
h? o j=2 
where 
72 jo2 
—|(-1) II {J (ars) J (rss) — (kr) 
h? 0 0 +) —I-} 


TL (reVir) (30) 


and 
72 

(—1) | IIL {J (ers) J — (kr) (Rrs—1) 
h? 0 0 —I-} —I-} I+} 


I+} 


If the range of the potential V(r) is not too great, one may use an approximation for ¢1(7;) ¢2(7s—1) 
— ¢2(r.)¢1(rs—-1), obtained by expanding the functions ¢; and ¢2 of 7,1, about the point r=r,. Re- 
taining only the first term, one has, 


where the Wronskian W(¢) is a constant. 

Three factors limit the magnitude of (30) and (31): the inequalities (27a), the range of V(r), 
and the form of the integrand, (32). The range is a constant, and one may easily see that the approxi- 
mation (32) is less serious the higher the order of a; or };. 


1 The asymptotic form of the wave function is periodic with respect to (kr—const.Xlog r). Reference 2, Chapter 
III; reference 10, Chapter XVI, for V~1/r. 
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6. Application to various models 


When the interaction potential is such that &, 
and 7. of (26) and (27) can be determined (c >0), 
it is possible to use simpler functions than those 
employed in (29). A more convenient set is 7/*! 
and r~'. With this choice one finds that 7,—TH, 
Now if Vir)=x, O0<r<e. 
Viry=> ur", c<r<d then one may set 
ne=—c' and calculate d&/dr, and 
én/dr, in known functions for r=d. At this 
point the solution may be joined to that given 
by the functions kr) and (kr) and again 


one may determine from the asymptotic behavior 
at infinity, the phase 6:. 

Instead of taking V(r) to be infinite for 
0<r<c, one may also assume V(r)=G;r’, 
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O0<r<c. Again the functions r'*', may be 
employed to determine & and »., repeating the 
previous analysis. 

It is possible to show that the series (26) and 
(27), with the above choice of functions, is 
absolutely convergent, being less than the corre- 
sponding series, near the origin 


r | ? 
k?—- 
0 h | 


—, 33 


One may, therefore, treat potentials, with these 
functions, which are singular at the origin, pro- 
vided such singularities are of an order less than 
the inverse square. 
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Letters to the Editor 


PHYSICAL REVIEW 


ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. Closing dates for this department are, for 
the first issue of the month, the eighteenth of the preceding 
month, for the second issue, the third of the month. Because 
of the late closing dates for the section no proof can be shown 
to authors. The Board of Editors does not hold itself responsible 
for the opinions expressed by the correspondents. 
Communications should not in general exceed 600 words 
in length. 


Radioactive Isotopes of Osmium 


G. T. SEABORG AND G. FRIEDLANDER 


Department of Chemistry, Radiation Laboratory, Department of Physics, 
University of California, Berkeley, California 
January 31, 1941. 


— has the following stable isotopes:! 


Mass 184-186 «1187's «188 189-190) 192 


Abundance 
in % 0.018 1.59 1.64 13.3 16.1 264 41.0 


Slow neutron bombardment of osmium should, therefore, 
be expected to yield at least two radioactive osmium iso- 
topes, Os"! and Os'. When the present work was started, 
only one osmium activity had been reported,? its half-life 
being given as 40 hours. Meanwhile Zingg* has reported 
two periods obtained by slow neutron bombardment of 
osmium, one of 29.8+1 hours and one of about 10 days. 

We bombarded osmium metal, surrounded by paraffin, 
with neutrons from the 60-inch cyclotron for several days. 
The chemically separated osmium fraction showed two 
periods, one of 32+2 hours and one of 17+1 days. From 
absorption measurements, the upper energy limits of the 
beta-particles from these two activities were found to be 
about 1.5 Mev and 0.35 Mev, respectively. Both activities 
have gamma-rays associated with them. 

How these two periods are to be distributed between 
Os"! and Os" cannot be decided by bombardment of neigh- 
boring elements with neutrons, protons, deuterons, or 
alpha-particles, assuming the ordinary nuclear reactions. 
A clue to the assignment might be expected from fast 
neutron bombardment of osmium. This should produce 
Os"™! by the reaction Os!(n,2m)Os'!, while Os" cannot be 
formed by a similar reaction. In addition, n,y reactions 
leading to both isotopes will occur, even with moderately 
fast neutrons. 

We surrounded one sample of osmium with boron carbide 
and cadmium, another one with paraffin, and bombarded 
them with beryllium+16-Mev deuteron neutrons under 
otherwise identical conditions. Both the 32-hour and the 
17-day periods were stronger in the slow neutron sample 
and the ratio of 32-hour to 17-day activity was twice as 
large with the fast as with the slow neutrons. The latter 
result may indicate that the 17-day period is due to Os" 
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and the 32-hour period to Os™, although this is not the 
only possible interpretation and the possibility of isomerism 
cannot be excluded. 

It is a pleasure to thank Professor E. O. Lawrence for his 
continued interest and to express our gratitude to the 
Rockefeller Foundation and the Research Corporation 
whose financial support made this work possible. 

1A. O. Nier, Phys. Rev. 52, 885 (1937). 

21. V. Kurtchatow, C. D. Latyschew, L. M. Nemenow and I. P. 


Selinow, Physik. Zeits. Sowjetunion 8, 589 (1935). 
3 E. Zingg, Helv. Phys. Acta 13, 219 (1940). 


Coulomb Exchange Energy in Light Nuclei 
M. PHILLIPS 
Brooklyn College, Brooklyn, New York 
AND 
E. FEENBERG 
New York University, New York, New York 
January 22, 1941 


NUMBER of isobaric mass differences of the type 

(Z+1,Z)—(Z,Z+1) are known with considerable 
accuracy. The recent determinations of Haxby, Shoupp, 
Stephens and Wells! yield 


(7) Be?7—Li?= 9.4+0.3 in 10-‘mu, 
(9) B°—Be® =11.6+0.2 “ 
(11) C"—B" =21.340.2 “ 
(13) N8—C8 =23.8+0.4 “ 


In addition? 
(15) O85 =30.8+2.0 


The large differences (11)—(9) and (15)—(13) and the 
small differences (9)—(7) and (13)—(11) find a simple 
qualitative explanation in the periodic variation of the 
Coulomb exchange energy as calculated in the Hartree 
approximation. * In the notation of references 3 and 4 the 
theoretical results are 

{Be?—Li?} — {Li'— = L.4+2K./3, 

{B°—Be®} — {Be?7—Li’} = L.—6K,., 

{Cu—B"} — {B*—Be®} =L.4+2K,/3, 

{N8—C8} —{C"—B"} =L.—6K,, 

{O'—N®} — =L.42K./3. 


Here L, and K, are the direct and exchange Coulomb 
integrals, respectively, for protons in the p shell. With the 
exception of (1) and (2) these equations are listed in refer- 
ence 4. The theory yields numerical values L.~4.9X 10 
mu, K.~0.30X10~* mu, but no quantitative significance 
is claimed for these numbers. The experimental mass 
differences are fitted reasonably well by L.=7.0 10™ mu 
and K,.=0.7X10-4 mu. In view of the crudeness of the 
Hartree approximation we prefer not to stress the possi- 
bility of numerical agreement, especially since supple- 
mentary effects, such as the one studied by Bethe,’ may 
account for part of the observed periodicity. 

1 Haxby, Shoupp, Stephens and Wells, Phys. Rev. 58, 1035 (1940). 

2 Fowler, Delsasso and Lauritsen, Phys. Rev. 49, 561 (1936). The 
upper limit of the position spectrum from O" is 1.7 Mev by inspection 


and 2.0 Mev by the K-U plot. We use the value 1.85 +0.20 Mev to com- 


pute the mass difference (15). 
3 E. Feenberg and E. Wigner, Phys. Rev. 51, 95 (1937). 
4E. Feenberg and M. Phillips, Phys. Rev. 51, 597 (1937). 
5H. A. Bethe, Phys. Rev. 54, 436 (1938). 
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Cosmic Rays at a Depth Equivalent to 1400 
Meters of Water 


Y. Nisuina, Y. Sekipo, Y. Miyazaki AND T. Masupa 


Institute of Physical and Chemical Research, Komagome, Hongo, 
Tokyo, Japan 


December 31, 1940. 


ARIOUS authors! have studied the nature of cosmic 
rays under very thick layers of matter. For the same 
purpose we have been performing experiments in Simidu 
tunnel (600 m above sea level, geographic ¢=36.8°N, 
\=138.9°E, geomagnetic ¢=26.6°N) of the government 
railways. An observation hut was built in one of the man- 
holes of the tunnel under a vertical thickness of about 
500 m of rocks, which consists mainly of diorite, the average 
density being probably 2.8, thus making the depth equiv- 
alent to about 1400 m of water. 

The experimental arrangements are shown in Fig. 1. 
We shall denote the coincidence discharge of any one 
Geiger-Miiller counter in each of the trays A, B and C, 
respectively, as being caused by a “single ray” and that 
of any one counter in each of the trays A and B and any 
two counters in the tray C, respectively, as being caused 
by a “shower.” It is thus possible that “single rays’ are 
sometimes showers in reality. These two types of coin- 
cidences are recorded simultaneously through two separate 
sets of amplifiers on the same paper tape, which is moved 
continuously by means of a clockwork. Lead absorbers of 
various thickness were placed either at the position I or 
II on a wood support 8 cm in thickness. Measurements 
were carried out for 4600 hours of observation between 
September, 1939, and July, 1940. 

The results of the experiments are shown in Fig. 2, from 
which we can draw the following conclusions. 

(1) There exist both “single rays’’ and “showers” which 
are hard and absorbed only gradually by 30 cm Pb. The 
relative abundance of “showers” to “singles” is higher than 
that at sea level, as was observed by various authors.! 
From this abundance, geometrical conditions and the above 
meaning of the word “‘singles”’ as well as from the simi- 
larity in the form of the absorption curves I for “singles” 
and “showers” in Fig. 2, we think that a greater part of 
“singles” are fragments of hard showers. We come thus to 
the conclusion that there exists a very high proportion of 


Fic. 1. Arrangement 
of counters. 


I Pb 


(Ocm 


hard showers at a depth of about 1400 m water equivalent. 
Barnéthy and Forré observed a sharp minimum of the 
absorption curve at 10 cm Pb, while our minimum lies 
within the range of statistical fluctuations. 

(2) As seen from Fig. 2, the position of lead either at I or 
at II makes no statistically significant difference in the 
number of coincidences. This fact seems to show the absence 
of non-ionizing primary rays for the hard ionizing showers 
which we observe, as mentioned above, in our apparatus. It 
is, however, also probable that the non-ionizing rays, even if 
they exist, do not produce the hard secondaries within 
such a thin layer of matter as 20 cm Pb, because the hard 
showers mentioned have an average range much larger 
than 30 cm Pb, and these must therefore have mostly 
been produced within a very thick layer of rock. It is thus 
difficult to decide from our results whether the primaries 


Single lay 


Shower 


= 
= 
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Fic. 2. Number of single rays and showers as a function of 
thickness of lead. 


for the hard showers observed are ionizing particles such 
as mesons or protons,? or non-ionizing particles—say neu- 
trinos—as was proposed by Barnéthy and Forré.’ 

This work was carried out as a part of the program of 
the Cosmic-Ray Sub-Committee of the Japan Society for 
the Promotion of Scientific Research. Our measurements 
were made possible by the kind assistance rendered by the 
Department of Railways and its various offices at Ueno, 
Minakami and Yusawa stations and we wish to express 
our cordial thanks to all officials concerned. We acknow!l- 
edge the kind help given by our colleagues Messrs. M. lio 
and T. Kameta for the preparations as well as for the 
observations throughout this work. 


1J. Barnéthy and M. Forré, Phys Rev. 55, 870 (1939); 58, 844 
(1940); V. C. Wilson, Phys. Rev. 55, 6 (1939); J. Clay, Rev. Mod. 
Phys. 11, 128 (1939). 

2J. Clay, reference 1. 

3J. Barnéthy and M. Forré, Nature 138, 325 (1936); Zeits. f. 
Physik 104, 744 (1937). 
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Change in Height of a Mesotron-Producing 
Layer of Air 
NIEL F. BEARDSLEY 


University of Chicago, Chicago, Illinois 
January 30, 1941 


N connection with the correlation studies on the effect 
of air conditions aloft on cosmic-ray intensity at the sur- 
face! it was convenient to apply the theory of Blackett? to a 
particular set of observations. Blackett assumed a meso- 
tron-producing layer of air at about 16 km which had a 
mean change of height between summer and winter of 500 
m and a mean temperature change of 10°C. From our data 
the mean density at this height is 0.1773 g-cm~. 
By linear interpolation of density and pressure between 
15-16 km and 16-17 km we computed the height Z, of this 
density for each of the 123 days in 1939 for which the radio- 
sonde flight at the Anacostia Naval Air Station reached 
the necessary height and also the Compton-Bennett 
shielded ionization meters were in operation at both Chelt- 
enham and at Huancayo, Peru. 

We made a least-squares determination of the constants 
in an equation connecting the variation of cosmic-ray 
intensity with (a) the mass of air above this mesotron- 
producing layer represented by P., (b) the height of the 
layer Z, and (c) the mass, (Po=P.), of air below this layer. 
We found 


= —0.935(Po— Pz) — 33.585Z, —1.156P, 
(pressures in mb) 


6J’/I’ is the variation from balance in units of 0.1 percent 
of the intensity at Cheltenham corrected for world-wide 
changes as measured by the meter at Huancayo. Let 6/’’/J’ 
be the variation in intensity because of the change in height 
of the mesotron-producing layer only. Then from 61’’/I’ 
= —6Z/L we get the mean path before disintegration of 
the mesotrons as L,=29.8 km. Blackett assumed 32 km. 
A two-constant equation for these same data gives 


6I'/I' = —0.926P)—32.355Z, 


with L,=30.9 km. In both these cases approximately 60 
percent of the variance of J’ is associated with the variables 
chosen. 

From our present knowledge of cosmic rays, it seems 
probable that the mean mesotron-producing layer ought 
to be associated with a given pressure rather than with a 
given density. However, using a constant pressure of 107 
mb we find 

= —0.946P)— 40.7952, 


with L,=24.52 km and only 50 percent of the variance of 
I’ associated with P» and Z,. 

These results indicate that Blackett’s assumption of a 
decaying mesotron formed at a height which changes from 
day to day is a possible explanation for some of the varia- 
tions observed in cosmic-ray intensity at the surface. 
Some method for using the probable fact that mesotrons 
are produced at varying rates at different levels will be 
needed to make this correlation more exact. 


1N. F. Beardsley, Phys. Rev. 59, 233 (1941). 
?P. M.S. Blackett, Phys. Rev. 54, 973 (1938). 


THE EDITOR 


Capture Cross Sections for Slow Neutrons 


K. SInMA AND F. YAMASAKI 


Nuclear Research Laboratory, Institute of Physical and Chemical 
Research, Tokyo, Japan 


January 8, 1941 


ECENTLY Rasetti! and Goldhaber and O’ Neal? have 
determined the capture cross sections of various 
isotopes for C neutrons by comparing the intensities of 
8-rays from radioactive isotopes produced by capture with 
that of an isotope Mn**, whose absolute cross section of 
formation was measured by Lapointe and Rasetti* and was 
found to be ¢=9.4X 

Since the summer of 1939, we have been measuring the 
capture cross sections for slow neutrons with a similar 
method by using a strong neutron source (Be+D) obtained 
with our cyclotron. We shall also give here our results. 
A paraffin cylinder 14 cm in diameter and 14 cm in height 
was placed about 30 cm away from the source, and sub- 
substances to be measured were placed at its center, to- 
gether with a standard sample. The activity of the sub- 
stance in question was then compared with that of the 
standard and thus the relative capture cross sections were 
determined for about fifty isotopes. The measurements of 
activity were made with a Geiger-Miiller tube counter or 
a Lauritsen electroscope, both with a window of aluminum 
foil 0.013 mm thick. 

Absolute capture cross sections were then calculated in 
the following manner. Dunning, Pegram, Fink and Mitchell‘ 
obtained large absorption cross sections for slow neutrons 
in the elements Rh, Ag, Re and Au. From the table of 
isotopes, we can conclude that these absorption processes 
are almost entirely due to the capture of slow neutrons. 
By comparing our results for these elements with the ab- 
solute cross sections of the above-mentioned authors, we 
can thus convert our relative cross sections into the ab- 
solute values. The results are given in Table I. 

From this we can see that Rasetti’s values are system- 
atically larger than ours. The origin of the discrepancies 


TABLE I. Conversion of our relative capture cross sections 
into absolute values. 


Iso- DEcAY IN 107*4| Iso- Decay IN 
TOPE PERIOD cm? TOPE PERIOD cm? 
14.8 hr. 0.38 Ru! 20 hr. 0.48 
Mg 10.2 min. 0.028 | Rhi 4.2 min. 8.8 
Al?? 2.4 min. 0.15 Rhtos 44 sec. 140 
Sise 170 min. 0.063 | Pdies 13 hr. 20 
pai 14.3 days 0.15 Pd'0 17 min. 0.55 
Cis? 37 min. 0.38 Agi? 2.3 min. 26.7 
K4 12.4 hr. 0.7 Agios 22 sec. 70 
Sc 85 days 2.8 In'3 48 days 24 
Visi 3.9 min. 3.5 In'5 54 min. 125 
Mn 2.59 hr. 6.0 Sb!2i 2.8 days 8.0 
Co%? 11 min. 48 Sb!3 60 days 1.6 
Co'?® 7 yr. 3.8 [27 25 min. 6.5 
Ni® 2.6 hr. 0.35 Balss 86 min. 3.0 
Cus 12.8 hr. 1.3 Lals9 31 hr. 5.5 
Cu 5 min. 1.5 Eu! 9.2 hr. 530 
Zn‘ 57 min 0.48 Dy64 2.5 hr. 1600 
Ga*® 20 min 0.68 Tals 97 days 4.5 
Ga?! 14 hr 1.3 W186 24 hr. 33 
26.8 hr 1.5 90 hr. 85 
Br? 4.4 hr 1.4 Re!87 16 hr. 63 
Br79 18 min 4.3 Iris 60 days 190 
Br®! 34 hr 1.1 Ir193 19 hr. 63 
ys 60 hr 0.73 Pti98 31 min. 9.5 
Nb* 6.6 min 0.005 | Au’? 2.7 days 130 
Rule 4hr 0.19 T1203 4 min. 0.19 
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probably lies in the difference in the absolute values taken 
for reference in converting the relative cross sections into 
absolute ones. 

A more detailed report is now in press and will soon 
appear in the Scientific Papers of the Institute of Physical 
and Chemical Research (Tokyo). 

We wish to express our best thanks to Dr. Y. Nishina 
for his valuable help and encouragement. 

1F. Rasetti, Phys. Rev. 58, 869 (1940). 

2M. Goldhaber and R. D. O'Neal, Phys. Rev. $9, 109A (1941). 

3C. Lapointe and F. Rasetti, ngs —— 58, 554 (1940 


). 
4jJ. R. Dunning, G. B. Pegram, G. A. Fink and D. P. Mitchell, 
Phys. Rev. 48, 265 (1935). 


Radioactive »,Sc"', ,,A® and 


D. R. anp L. D. P. KincG 
Purdue University, West Lafayette, Indiana 
January 31, 1941 


HE production of these Wigner-type nuclei has 

already been reported.' Energy measurements of 
the emitted positrons have now been completed by means 
of a cloud chamber with air and alcohol vapor at atmos- 
pheric pressure. 

Targets were held in the beam and then placed in a 
holder which could be swung up to a thin window on the 
side of the cloud chamber. A picture was taken within 
two seconds after the sample was removed from the beam. 
The cyclotron was on only during the actual bombard- 
ments which varied from 0.5 to 5 seconds, depending on 
the sample used. Only one photograph was made for each 
bombardment, and four bombardments were made per 
minute. To eliminate the influence of longer periods in- 
duced in the samples, six samples were used successively. 

The half-life for »Sc* was reported as 0.87+0.03 sec- 
ond.'! Figure 1 shows the spectrum for the positrons ob- 
served from the reaction 2Ca*(d,n)2Sc. The magnetic 
field used was 1276 oersteds; the observed p maximum was 
142 mm. The upper energy limit obtained for the »:;Sc* 
positrons is 4.94+0.07 Mev. 

The maximum energy of the positrons and the half-life 
of A® and S* have been recently also reported by the 
Princeton group.? Their value of 3.2+0.02 seconds for the 
half-life and 3.85+0.07 Mev for the upper limit of the 


THE EDITOR 403 


Positrons from 
w 
/ 


150} / 


/ 
2 
/ 
‘= / 
/ 
t 
5 / 
/ 
| 
= 
(4) i i i i 
20 50 


Radius eof Curvature in Millimeters 


Fic. 1. Positrons from Sc*!. Maximum energy 4.94+0.07 Mev, 
half-life 0.87 +0.03 second. 


positrons for S* are in good agreement with our values, 
3.18+0.04 seconds! * and 3.87+0.15 Mev, respectively. 

For A® the value of 2.2+0.2 seconds for the half-life 
and 4.38+0.07 Mev for the maximum positron energy, 
obtained at Princeton are to be compared with our values 
of 1.91 seconds and 4.41+0.09 Mev." * Since the half-life 
for A® did not agree with the Princeton results, this period 
was remeasured with greater accuracy giving 1.88+0.04 
seconds. Because of the high energy of the positrons from 
all of these isotopes, it has been possible to place suitable 
filters between the samples and the counter so that prac- 
tically pure periods were obtained. 

All of these results are in good agreement with calculated 
values.? 

'L. D. P. King and D. R. Elliott, Phys. Rev. 59, 108A (1941). 


* White, Cruetz, ry and Wilson, Phys. Rev. 59, 63 (1941). 
3L. D. P. King and D. R. Elliott, Phys. Rev. 58, 846 (1940). 
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